请教高人，我的程序哪里出错了
 

该程序是为研究n个独立同正态分布的随机变量的积的方差与n的关系，我的程序考虑的积的方差课能也与原正态分布的方差有关，利用随机数产生器产生一定量的数据拟合
muo=7.4863;
sigmao=9.8520;
format long
muo=muo/1000+1;
sigmao=sigmao/1000;
j=0;
for i=0.00012:0.00002:0.01
a=normrnd(muo,i,500);
j=j+1;
z=eye(500);
for n=1:500
z=z.*a;
t=z(:);
[mu,sigma]=normfit(t);
s(j,n)=sigma;
end
end
sn=s(:);
j=1:500;
j=j';
jn=repmat(j,500,1);
n=1:500;
nn=repmat(n,500,1);
nn=nn(:);
surf(sn,jn,nn);
这里尚未完成的到的2500组数据拟合方法也请教一下，还有拟合出以后可不可以取定一个原正态分布的方差来确定n个随机变量的积方差和n的关系