Example 5

%matplotlib inline
import sympy as sm
from sympy.physics.continuum_mechanics.beam import Beam

E = sm.Symbol('E')
I = sm.Symbol('I')
b = Beam(10, E, I)
b.apply_support(0, type="pin")
b.apply_support(10, type="pin")
R_0, R_10 = sm.symbols('R_0, R_10')
b.solve_for_ild_reactions(-1, R_0, R_10)
b.ild_reactions

{R_0: 1 - x/10, R_10: x/10}

b.plot_ild_reactions();