sjtdf_zuoye2.asv

来自「此文件包括蒙特卡罗方法的基本思想及应用」· ASV 代码 · 共 14 行

function result=sjtdf(a,b,m,mm)
%a1是外重积分的下限
%b1是积分的上限
%m是函数的上界%积分函数cos(x)+2
%mm 是随机实验次数
frq=0;
xrandnum = unifrnd(a,b,1,mm);
yrandnum = unifrnd(0,m,1,mm);
for ii=1:mm
    if (cos(xrandnum(1,ii))+2>=yrandnum(1,ii))
        frq=frq+1;
    end    
end    
result=frq*m*(b-a)/mm