ybjzf1.m

此文件包括蒙特卡罗方法的基本思想及应用

function result=ybjzf1(a,b,m,mm)
%a是积分的下限
%b是积分的上限
%积分函数
%mm 是随机实验次数
sum=0;
xrandnum = unifrnd(a,b,1,mm);
for ii=1:mm
    sum=sum+sqrt(1-xrandnum(1,ii)^2);
end    
result=sum*(b-a)/mm;
result=result*4