怎么就是得不出正确结果呢？

fcode

a(imax,j)=m(k,i)
改为
a(imax,j)=m(k,j)

你的代码可以进行很大的优化：

[Fortran] 纯文本查看 复制代码

program S_equation  implicit none
  integer::i,k,imax,t
  real::max,n
  !用矩阵（实型动态数组）将线性方程组表示出来!
  real,dimension(:,:),allocatable::a,b
  real,dimension(:),allocatable::x,m
  write(*,*)'请输入线性方程组的元数n:'
  read(*,*)t
  !给动态数组分配内存!
  allocate(a(t,t+1),m(t+1),x(t),b(t,t+1))
  write(*,*)'请依次按行输入线性方程组的增广矩阵a：'
  read(*,*)( b(i,:),i=1,t)
  a=b
  do k=1,t-1
    !选取每列的最大列元素!
    max=maxval(a(k:,k))!abs(a(k,k))
    imax=maxloc(a(k:,k),1) + k - 1    
    !将最大列元素所在的行与第K行进行交换!
    m(k:)    = a(k,k:)
    a(k,k:)    = a(imax,k:)
    a(imax,k:) = m(k:)
    !对方程组按X1,X2，...，Xt的顺序进行依次消元，将增广矩阵化为上三角矩阵!
    do i=k+1,t
      a(i,k+1:)=a(i,k+1:) - a(i,k)/a(k,k)*a(k,k+1:)
    end do
  end do
  !先算出Xt!
  x(t)=a(t,t+1)/a(t,t)
  !回带!
  do k=t-1,1,-1
    n= sum( a(k,k+1:t)*x(k+1:) )
    x(k)=(a(k,t+1)-n)/a(k,k)
  end do
  write(*,*)'线性方程组的解为：'
  write(*,*) x
  write(*,*) matmul( b(1:t,1:t) , reshape(x,[t,1]) )!//验证
  read(*,*)
end program S_equation