montecarlo_support.m

function [mean1, stand_dev ]=montecarlo_support(N)
%Uses monte carlo model to input number of randomn varaibles to generate varying variables b and
%h so that they have a standard deviation
%The maximum deflection is then solved for and the mean and standard deviation of that is returned.
E=70e9;
q=50;
for i=1:N
b=normrnd(0.1,0.0001);
h=normrnd(0.01,0.00001);
I=(b*h^3)/12;
D= (q/(384*E*I))-(q/(96*E*I))+(q/(48*E*I));
d(i)=D;
end
mean1=mean(d);
stand_dev=std(d);
end