slvsm2.for

来自「Numerical Recipes一书中例子的源码所用到的函数集,William」· FOR 代码 · 共 12 行

      SUBROUTINE slvsm2(u,rhs)
      DOUBLE PRECISION rhs(3,3),u(3,3)
CU    USES fill0
      DOUBLE PRECISION disc,fact,h
      call fill0(u,3)
      h=.5d0
      fact=2./h**2
      disc=sqrt(fact**2+rhs(2,2))
      u(2,2)=-rhs(2,2)/(fact+disc)
      return
      END