pufeng.m

来自「本文件包括葡丰概率实验测定pi的方法介绍和matlab实现方法」· M 代码 · 共 18 行

function [pai,M]=pufeng(L,a,N)
if(L>=a)
    fprintf('参数设置错误，无法计算！\n');
    pai=0;M=0;
    return;
end
M=0; 
for k=1:N 
    x=a/2*rand(1); 
    p=pi*rand(1); 
    if x<0.5*L*sin(p) 
        M=M+1; 
    end 
end 
p=M/N; 
pai=(2*L)/(p*a);
fprintf('细针长度为%g,平行线间宽度为%g，\n共试验%d次，有%d次细针与平行线相交，\n算得的圆周率的近似值为%g\n'...
    ,L,a,N,M,pai);