share/colnew/ex2/dfsub.f

   1 C................................................................       DFS   10
   2       SUBROUTINE DFSUB(X, Z, DF)                                        DFS   20
   3       implicit real*8 (a-h, o-z)
   4       DIMENSION Z(4), DF(2,4)
   5       COMMON EPS, DMU, EPS4MU, GAMMA, XT
   6       DF(1,1) = 1./X/X + (1.+Z(3)/X)/EPS4MU
   7       DF(1,2) = -1./X
   8       DF(1,3) = -(1.-Z(1)/X)/EPS4MU
   9       DF(1,4) = 0.
  10       DF(2,1) = (1.-Z(1)/X)/DMU
  11       DF(2,2) = 0.
  12       DF(2,3) = 1./X/X
  13       DF(2,4) = -1./X
  14       RETURN
  15       END