MATLAB蒙特卡洛模拟代码求助

b1=randn(361,5000);
a1=randn(361,5000);
h(1,1)=0;
for i=1:361;
    j=1:5000;
    h(i+1,j)=h(i,j)+exp(-h(i,j))*0.0002*(-0.0015972-0.0002*exp(h(i,j)))-0.5*1.9933*1.9933*exp(-h(i,j))/360+exp(-0.5*h(i,j))*0.105053841*sqrt(1/360)*a1;
    y(i,j)=-0.3493*exp(h(i,j))+0.3184*(-0.3175)*((h(i+1,j)-h(i,j))*exp(h(i,j))-0.0002*(-7.986-exp(h(i,j)))+0.5*1.9933*1.9933/360)/1.9933+0.3184*exp(0.5*h(i,j))*sqrt(1-0.3175*0.3175)*sqrt(1/360)*b1;
end
end
plot(1:length(y(i,j)),y(i,j))
??? Index exceeds matrix dimensions.
求助
这段代码为啥没法输出图形，我知道是矩阵溢出，但是具体是哪里的错误，能告诉我吗

你是直接打在命令窗口吧！建议你还是在M文件中编写吧！首先程序多了一个end。然后是h(i+1,j)=h(i,j)+exp(-h(i,j))*0.0002*(-0.0015972-0.0002*exp(h(i,j)))-0.5*1.9933*1.9933*exp(-h(i,j))/360+exp(-0.5*h(i,j))*0.105053841*sqrt(1/360)*a1;这个公式有问题！指数超过矩阵尺寸！你自己再检查检查！