est1.m

motecalo仿真实现源代码

m=5;
n=500;
u=rand(n,m)
uu=1./(1+u.*u);
data=zeros(n,m);
data(1,:)=4*uu(1,:);
for j=2:n
    data(j,:)=4*sum(uu(1:j,:))/j;
end

est=data(n,:);
est1=sum(est)/m;

plot(data,'k')