如前贴，由于矩阵太大，出现内存分配不足问题，所以我暂时更改了矩阵大小。
当执行时，出现如果所示提示。是不是数据寻找出现了问题？
源代码如下： parameter(np=4100)
        integer n,np,indx(np)
        real*8 a(np,np),a2(np,np),y(np,np),prod(np,np),d
        character result*8
        open(unit=8,file='d:\matrixprogram\matrix.dat',status='unknown')
        open(unit=10,file='d:\matrixprogram\prepared.dat',status='old')

其中prepared.dat是已经存在且待变换的矩阵。

大家帮忙参考一下，谢谢！