[zhyanchao]

Z(j)=Z(j)-M(j)+M(j+1); 在j=100时M(j+1)超出了矩阵的大小。
另外你这个循环有问题，这个循环只能看到i=1000是的结果。应该在最后一个end前加个你要观看的矩阵，不要加分号。