c9_example5.asv

来自「motecalo仿真实现源代码」· ASV 代码 · 共 14 行

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()