anchor.diff:   50   IF (X .GT. XMAX) CALL XERMSG ('SLATEC', 'DGAMMA',
        +   'X SO BIG GAMMA OVERFLOWS', 3, 2)
   C
  -      DGAMMA = 0.D0
  +      UTIL_DGAMMA = 0.D0
         IF (X .LT. XMIN) CALL XERMSG ('SLATEC', 'DGAMMA',
        +   'X SO SMALL GAMMA UNDERFLOWS', 2, 1)
         IF (X.LT.XMIN) RETURN
   C
  -      DGAMMA = EXP ((Y-0.5D0)*LOG(Y) - Y + SQ2PIL + D9LGMC(Y) )
  +      UTIL_DGAMMA = EXP ((Y-0.5D0)*LOG(Y) - Y + SQ2PIL + D9LGMC(Y) )
         IF (X.GT.0.D0) RETURN
   C
         IF (ABS((X-AINT(X-0.5D0))/X) .LT. DXREL) CALL XERMSG ('SLATEC',
  
positive.code: C   891214  Prologue converted to Version 4.0 format.  (BAB)
  C   900315  CALLs to XERROR changed to CALLs to XERMSG.  (THJ)
  C   920618  Removed space from variable name.  (RWC, WRB)
  C***END PROLOGUE  DGAMMA
        DOUBLE PRECISION X, GAMCS(42), DXREL, PI, SINPIY, SQ2PIL, XMAX,
       1  XMIN, Y, D9LGMC, DCSEVL, D1MACH
        LOGICAL FIRST
  
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
  