MATLAB中索引超过数组元素个数怎么修改，直接告诉我修改代码吧，这个程序应该如何修改，谢谢您

这个程序是蒙特卡洛模拟光子的米氏散射，不知道该如何调试和应用，感谢各位

第23行的z(1)应该改为z(i)，否则z(i)就没有定义了。

第23行 z(i)=z(i-1)+L*cos(b); %把在z(1)改成在z(i),i变量。
不再报错，但运行时间有点长，不知何时结束。crtl+c中断后，显示计算了不少值。