shape_simple_support.m

function w =shape_simple_support(x,q)
%The function above calculated the displacement w(x) from an input of q and
%x where q is the transverse distributed load and x is the position on the
%beam.
%System Parameters 
E=70E9;%Young's modulus (70 GPa)
L=1; %Beam length (meters)
b=0.1; %Beam width (meters)
h=0.01; %Beam height (meters)

I=(b*(h^3))/12;%Second moment of inertia
%Function in terms of x 
w =(((q*L*(x^3))/(12*E*I)))-(((q*(x^4))/(24*E*I)))-((q*(L^3)*x)/(24*E*I));
end