sjtdf_pi2.m

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

function piguji=sjtdf_pi2(mm)
%mm 是随机实验次数
frq=0;
xrandnum = unifrnd(0,1,1,mm);
yrandnum = unifrnd(0,1,1,mm);
for ii=1:mm
    if (sqrt(1-xrandnum(1,ii)^2)>=yrandnum(1,ii))
        frq=frq+1;
    end    
end    
piguji=4*frq/mm