import math
from sympy import symbols
#Variable declaration
w=symbols('w')
L=symbols('L')
E=symbols('E')
I=symbols('I')
Rb=symbols('Rb')
n=-1
#Calculation
yBW=(n*w*(L**4))/(8*E*I)
yBR=(Rb*(L**3))/(3*E*I)
Rb=(3/8.0)*(w)*(L)
Ra=w*L - Rb
Ma=(1/2.0)*(w)*(L**2) - Rb*L
# Result
print ('Reactions at the supports for the prismatic beam :-')
print(Ra)
print ('Reactions at the supports for the loading :-')
print(Rb)