Question

Fortran：program exception--array bounds exceeded

Program linear_equation
Implicit none
Integer::i,j,k,imax,t
Real::max,n
Real,dimension(:.:),allocatable::a,m
Real,dimension(:),allocatable::x
Print*,’输入方程组维数n:’
Read*,t
Allocate(a(t,t+1),m(t,t+1),x(t))
Print*’按行输入方程组的增广矩阵a:’
Read*,((a(i,j),j=1,t+1),i=1,t)
Do k=1,t-1
Max=abs(a(k,k))
Imax=k
Do i=k+1,t+1
If (abs(a(i,k))>max) then
Max=abs(a(i,k))
Imax=i
End if
End do
Do j=k,t+1
M(k,j)=a(k,j)
A(k,j)=a(imax,j)
A(imax,j)=m(k,j)
End do
Do i=k+1,t
M(i,k)=a(i,k)/a(k,k)
Do j=k+1,t+1
a(i,j)=a(i,j)-m(i,k)*a(k,j)
End do
End do
End do
X(t)=a(t,t+1)/a(t,t)
Do k=t-1,1,-1
N=0
Do j=t,k+1,-1
N=n+a(k,j)*x(j)
End do
X(k)=(a(k,t+1)-n)/a(k,k)
End do
Print*,’方程组的解为：’
Do i=1,t
Print*,x(i)
End do
End program linear_equation
编译之后都没有问题，但是运行之后就会显示
program exception--array bounds exceeded

，求大神帮我看一看，具体能怎么解决，

Answer 1

Intel(R) Visual Fortran Compiler XE 12.0.5.221 [IA-32] 编译后，

输入方程组维数n:
2
按行输入方程组的增广矩阵a:
3 4 5
2 3 5
方程组的解为：
2.4305608E+25
-1.6203738E+25
请按任意键继续. . .

运行正常

追问

这。。运行环境不一样。请问能解决下我程序的问题吗？有问题存在还是得搞清楚才好啊，谢谢你

追答

如果你的编译器上会出错，那么自己头从下断点，单步Debug调试。
我的编译器是没问题的，所以无法替你完成这工作。

Answer 2

Do i=k+1,t+1
If (abs(a(i,k))>max) then
Max=abs(a(i,k))
Imax=i
End if
End do

a矩阵是t*（t+1）的，而abs(a(i,k))>max这一段里，当i=t+1时，出现了a(t+1,k)，所以数组越界