function?MVFOSM（NE）
%均值一次二階矩法
%適用范圍：隨機變量為正態分布，變量間獨立
%特點：線性功能函數獲得失效概率精確解，非線性程度不高的功能函數獲得失效概率近似解
%缺點：對于物理意義相同而數學表達是不同的非線性問題可能得到完全不同的失效概率
%輸入參數：NE?-?極限狀態函數索引，整數
????global?Prob
????tic
????if?nargin<1?NE?=?0;?end????????????????????????????????%如果未定義方程索引，默認使用索引0
????
????%1.定義符號變量
????syms?x1??x2??x3??x4??x5??x6?x7??x8??x9??xa??xb??xc??xd??xe??xf??xg;
????X=[?x1??x2??x3??x4??x5??x6?x7??x8??x9??xa??xb??xc??xd??xe??xf??xg?];
????XX=[‘x1‘;??‘x2‘;??‘x3‘??;‘x4‘?;?‘x5‘;??‘x6‘?;‘x7‘;??‘x8‘;??‘x9‘;??‘xa‘;??‘xb‘;??‘xc‘;??‘xd‘;??‘xe‘;??‘xf‘;?‘xg‘];

????%2.獲得方程相關公式、變量數目、變量分布、變量參數
????Prob?=?FunDist（?NE?）;
????
????%3.判斷方法適用性
????for?i=1:Prob.Nx
????????if?strcmp（Prob.Dist{i}?‘norm‘）==0?&&?strcmp（Prob.Dist{i}?‘Normal‘）==0
????????	error（‘This?method?can?only?be?used?for?normal?dis