slvsml.for

Numerical Recipes一书中例子的源码所用到的函数集,William H. Press 和 Saul A. Teukolsky 所著

      SUBROUTINE slvsml(u,rhs)
      DOUBLE PRECISION rhs(3,3),u(3,3)
CU    USES fill0
      DOUBLE PRECISION h
      call fill0(u,3)
      h=.5d0
      u(2,2)=-h*h*rhs(2,2)/4.d0
      return
      END