Example 6

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

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

{R_0: 3/2 - 3*x/20, M_0: x/2 - 5, R_10: 3*x/20 - 1/2}

b.plot_ild_reactions();