my_pi_montocalo.m

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

function y=my_pi_montocalo(n)
m=0;
for k=1:n
    if rand(1)^2+rand(1)^2<=1
        m=m+1;
    end
end
m=4.*m;
vpa(m/n)