ybjzf_zuoye2.m

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

function result=ybjzf_zuoye2(a1,b1,a2,b2,mm)
%a1是外重积分的下限
%b1是外重积分的上限
%a2是内重积分的下限
%b2是内重积分的上限
%积分函数cos(x)+sin(y)+3
%mm 是随机实验次数
xrandnum = unifrnd(a2,b2,1,mm);
yrandnum = unifrnd(a1,b1,1,mm);
fsum=sum(cos(xrandnum)+sin(yrandnum)+3);
result=fsum*(b1-a1)*(b2-a2)/mm

return
 ybjzf_zuoye2(0,4,1,3,10000)
result =  24.4723