clc
clear all
% beem
tic
l=4.5;
sx=240;
Fy=2400;
Q=1.55;
sigmaQ=.155;
%sigmaQ=.2*1.55;
%sigmaQ=.4*1.55;
sigmaFy=240;
%sigmaFy=.2*240;

    
n=100;
fy=random('normal',[2344.2 2455.8],n,1);
%fy1=2344.2+(2455.8-2344.2)*fy;
    
    q=random('normal',[1.5139 1.5861],n,1);
   %q1=1.5139+(1.5861-1.5139)*q;
sum=0;
for i=1:n
    

    
    g=(sx*fy)-(q*(l^2/8)*10^5);
    if g<0
        s=1;
    else
        s=0;
    end
    sum=sum+s;
end
pf=sum/n

toc