from sympy.physics.continuum_mechanics.beam import Beam
from sympy import symbols
E, I = symbols('E, I')
b = Beam(6, E, I)
b.apply_support(0, "pin")
b.apply_support(3, "roller")
b.apply_load(4, 3, 0, end=6)
b.apply_load(4, 0, 1, end=3)
R_0, R_3 = symbols('R_0, R_3')
b.solve_for_reaction_loads(R_0, R_3)
p = b.draw()
p.show()
b.plot_shear_force()