How to do failure probability estimation using the direct Monte Carlo simulation. my code does not run it gives me "untitled'' as ans

%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE

%The given limit state function for the given beam is

%g = -3.689 + 0.0184*D^2 + 0.2938*V ;

%Given :

D = 90.0219343 ;

V = 6.672408424 ;

mean_D = 90 ;

mean_V = 6.6724 ;

sigma_D = 10.29023 ;

sigma_V = 0.67708 ;

%Now use normrnd function

D = normrnd(mean_D, sigma_D, [n,1]);

V = normrnd(mean_V, sigma_V, [n,1]);

%

g = zeros(n,1);

K = 0;

%Create for loop

for i=1:n

g(i,1) = -3.689 + 0.0184*(D(i,1)^2) + 0.2938*(V(i,1)) ;

if g(i,1)<0

K = K + 1;

end

end

P_fail = K/n;

Answer 1

[Torsten]

%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE
%The given limit state function for the given beam is 
%g = -3.689 + 0.0184*D^2 + 0.2938*V ;
%Given :
D = 90.0219343 ;
V = 6.672408424 ;
mean_D = 90 ;
mean_V = 6.6724 ;
sigma_D = 10.29023 ;
sigma_V = 0.67708 ;
%Now use normrnd function
n = 1e6;
D = normrnd(mean_D, sigma_D, [n,1]);
V = normrnd(mean_V, sigma_V, [n,1]);
%
%Create for loop
K = nnz(-3.689 + 0.0184*D.^2 + 0.2938*V < 0);
P_fail = K/n

P_fail = 0

[Ashwini Subhash]

Thanks it worked. I also wanted to know how can I calculate beta function and show its histogram

[Ashwini Subhash]

Actually I ran first order reliability method first so I need to compare beta function and pf of FORM and Monte Carlo Simulation

[Torsten]

first order reliability method ?

compare beta function and pf of FORM ?

[Ashwini Subhash]

Forget about FORM. I was asking how can I find beta distribution for the given input parameters like in the attached file