%問題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;???