shape_simple_support.m

function w =shape_simple_support(x, q)
%The Differential equation for the shape of the beam was integrated to
%solve for the equation of the beam. The input is the postition on the beam
%and the constant distributed load and the output is the shape of the beam.
%The plot shows the maximum deflection that occurs at loads ranging from 1 to 1000 
E=70e9;
b=0.1;
h=0.01;
I=(b*h^3)/12;
w=q*((x^4)/(24*E*I))-((q*(x^3))/(12*E*I))+((q*x)/(24*E*I));
q1=1:1:1000;
for n= 1:1:1000
d= (n/(384*E*I))-(n/(96*E*I))+(n/(48*E*I));
D(n)=d;
end
plot(q1,D,'.');
xlabel 'Distributed Load'
ylabel 'Maximum Deflection'
title('q vs Maximum Deflection')