buffon.m

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

function piguji=buffon(llength,mm)
%llength 是针的长度
%mm 是随机实验次数
frq=0;
xrandnum = unifrnd(0,0.5,1,mm);
phi= unifrnd(0,pi,1,mm);
for ii=1:mm
    if (xrandnum(1,ii)<=(llength*sin(phi(1,ii))/2))
        frq=frq+1;
    end    
end    
piguji=2*llength/(frq/mm)