JC法可靠度计算

%问题1：
%某受弯构件极限状态方程为：Z = g（R SG SQ） = R－SG－SQ，
%R服从LN（μσ），μR = 155.7，δR = 0.15
%SG服从N（μSGσSG），μSG = 30，δSG = 0.11；
%SQ服从极值Ⅰ型分布，μSQ = 40，δSQ = 0.15；
%用monte carlo法求受弯构件可靠性指标β。
% R=[3039.9	3018.8	3001.2	2939.6	2912.5 2872.1	2851.7	2848.6	2805.5	2637.4];
% P1=[335.2 371.8	268.5 296.8 315.6 301.5 335.5 33.6 352.4 344.1];
% QC=1058.67;
% P2=91.23;
% ZZ=959.75;
% P3=41.26;
function [bb]=kkdmc
mm=1000;
nn=mm^2;