mybuffon.m

来自「algotithm to calculate pi」· M 代码 · 共 12 行

function y=mybuffon(n)
count=0;
for k=1:n
    a=2*rand(1);
    b=2*rand(1);
    if a>=1&&b<=1 count=count+1;
    elseif a>=1&&b>1;
    elseif a<1&&b>=1 count=count+1;
    else;
    end
end
vpa(2*n/count)