消主元解方程的程序问题

Fortran会长大

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

subroutine ch(iz)
  parameter nub0=2400
  common/a4/rm(nub0,nub0),x(nub0),u(nub0)
  write(*,*)rm
  do 20 ii=1,iz-1
   c=1.0/rm(ii,ii)
   do 15 ij=ii,iz
15        rm(ii,ij)=rm(ii,ij)*c
   u(ii)=u(ii)*c
   do 20 jj=ii+1,iz
    do 30 kk=ii+1,iz
    rm(jj,kk)=rm(jj,kk)-rm(ii,kk)*rm(jj,ii)
30        continue
   u(jj)=u(jj)-u(ii)*rm(jj,ii)
20       continue
  x(iz)=u(iz)/rm(iz,iz)
  do 10 ii=iz-1,1,-1
   as1=0.0
   do 40 jj=ii+1,iz
    as1=as1+x(jj)*rm(ii,jj)
40        continue
   x(ii)=u(ii)-as1
10       continue
      do 55 ii=1,iz-1
  write(*,*)u(ii)
55      continue
  return
end

这是一个消主元解微分方程的程序 主程序里rm已经设定为系数矩阵了，为什么我读u（ii）能读出  但是读x（ii）就结果就是NAN呢