montecarlo_model.m

function [Mean, standard_deviation ]=montecarlo_model(N)
%determines the mean and standard deviation of the 
%maximum deflection of the beam. b and h are random variables 
%with 0.1% standard deviation and q=50 N/m
E=70e9; %Youngs Modulus 70 GPA
q=50; %Input of loading 50 N/m
for i=1:N
b=normrnd(0.1,0.01);
h=normrnd(0.01,0.001);
I=(b*h^3)/12;
L=1;
x=0.5;
W =-(((q*L*(x^3))/(12*E*I)))-(((q*(x^4))/(24*E*I)))-((q*(L^3)*x)/(24*E*I));
w(i)=W;
end
Mean=mean(w);
standard_deviation=std(w)
end