c9_estimatepi.m

motecalo仿真实现源代码

% File c9_estimatepi.m

m=input('Enter M,number of experiments>');% m个pi的估计值
n=input('Enter N,number of trials / experiments> ');% 得到每个估计值所作重复随机试验的次数

z=zeros(1,m);
data=zeros(n,m);
for j=1:m
    x=rand(1,n);
    y=rand(1,n);
    k=0;
    for i=1:n
        if x(i)^2+y(i)^2<=1  %  点落在四分之一圆内
            k=k+1;
        end
        data(i,j)=4*(k/i);
    end
    z(j)=data(n,j)
end
    plot(data,'k')
    xlabel('number of trials')
    ylaber('estimate of pi')
    % End of script file.