shape_simple_support_plot.m

%The objective of this code is to plot the load distributed on the beam
%versus the maximum deflection  
max_load=100;
q=[1:max_load];
w=zeros(1,100);
x=0.5; %Location where max deflection occurs 

for n=1:max_load
    w(n) = shape_simple_support(x,n);
end

%Plot load vs max deflection
plot(q,w)
setdefaults 
xlabel('Applied Load (N/m)') %x-axis label
ylabel('Max Deflection (m)') %y-axis label
title('Applied Load versus Maximum Deflection')