可靠性方法：蒙特卡洛法的Matlab源代码

function MCS（NEN）
%蒙特卡罗法
%适用范围：随机变量为任意分布，变量间独立
%适用范围：可进行串联多失效模式计算
%输入参数：NE - 极限状态函数索引，整数
%输入参数：N  - 抽样数量，整数
    global Prob
    clc
    tic
    if nargin<1 NE = 0; end                                 %如果未定义方程索引，默认使用索引0
    if nargin<2 N = 10000000; end                           %如果未定义抽样数量，默认使用10000000
    Prob = FunDist（ NE ）;                                    %获得方程相关公式、变量数目、变量参分布、变量参数
    x = XSamp（ProbN）;                                       %根据变量数目、分布、参数，以及抽样数量，生成变量随机值
   
    %计算系统可靠性
    for i=1:Prob.Ny
        gValue（:i）=eval（Prob.Fung{i}）;                      %根据方程公式计算函数值
    end
    gmin=min（gValue[]2）;                                   %假设所有失效模式为串联关系，获取所有函数值小于零索引
    XI=find（gmin<=0）;                                        %搜索所有gmin值小于零的索引
    
    %