share/colnew/ex2/solutn.f

   1 C................................................................       SOL   10
   2       SUBROUTINE SOLUTN(X, Z, DMVAL)                                    SOL   20
   3       implicit real*8 (a-h, o-z)
   4       COMMON EPS, DMU, EPS4MU, GAMMA, XT
   5       DIMENSION Z(4), DMVAL(2)
   6       CONS = GAMMA*X*(1.-.5*X*X)
   7       DCONS = GAMMA*(1.-1.5*X*X)
   8       D2CONS = -3.*GAMMA*X
   9       IF (X.GT.XT) GO TO 10
  10       Z(1) = 2.*X
  11       Z(2) = 2.
  12       Z(3) = -2.*X + CONS
  13       Z(4) = -2. + DCONS
  14       DMVAL(2) = D2CONS
  15       GO TO 20
  16    10 Z(1) = 0.
  17       Z(2) = 0.
  18       Z(3) = -CONS
  19       Z(4) = -DCONS
  20       DMVAL(2) = -D2CONS
  21    20 DMVAL(1) = 0.
  22       RETURN
  23       END