anchor.diff:  C***AUTHOR  Fullerton, W., (LANL)
   C***DESCRIPTION
   C
  -C DGAMMA(X) calculates the double precision complete Gamma function
  +C UTIL_DGAMMA(X) calculates the double precision complete Gamma function
   C for double precision argument X.
   C
   C Series for GAM        on the interval  0.          to  1.00000E+00
  
positive.code:       DATA PI / 3.1415926535 8979323846 2643383279 50 D0 /
        DATA SQ2PIL / 0.9189385332 0467274178 0329736405 62 D0 /
        DATA FIRST /.TRUE./
  C***FIRST EXECUTABLE STATEMENT  DGAMMA
        IF (FIRST) THEN
           NGAM = INITDS (GAMCS, 42, 0.1*REAL(D1MACH(3)) )
  C
  
negative.code: C   891214  Prologue converted to Version 4.0 format.  (BAB)
  C   900727  Added EXTERNAL statement.  (WRB)
  C***END PROLOGUE  DGAMR
        DOUBLE PRECISION X, ALNGX, SGNGX, UTIL_DGAMMA
        EXTERNAL UTIL_DGAMMA
  C***FIRST EXECUTABLE STATEMENT  DGAMR
        DGAMR = 0.0D0
        IF (X.LE.0.0D0 .AND. AINT(X).EQ.X) RETURN
  