dgami.f

Go to the documentation of this file.

*DECK DGAMI
      DOUBLE PRECISION FUNCTION dgami (A, X)
C***BEGIN PROLOGUE  DGAMI
C***PURPOSE  Evaluate the incomplete Gamma function.
C***LIBRARY   SLATEC (FNLIB)
C***CATEGORY  C7E
C***TYPE      DOUBLE PRECISION (GAMI-S, DGAMI-D)
C***KEYWORDS  FNLIB, INCOMPLETE GAMMA FUNCTION, SPECIAL FUNCTIONS
C***AUTHOR  Fullerton, W., (LANL)
C***DESCRIPTION
C
C Evaluate the incomplete gamma function defined by
C
C DGAMI = integral from T = 0 to X of EXP(-T) * T**(A-1.0) .
C
C DGAMI is evaluated for positive values of A and non-negative values
C of X.  A slight deterioration of 2 or 3 digits accuracy will occur
C when DGAMI is very large or very small, because logarithmic variables
C are used.  The function and both arguments are double precision.
C
C***REFERENCES  (NONE)
C***ROUTINES CALLED  DGAMIT, DLNGAM, XERMSG
C***REVISION HISTORY  (YYMMDD)
C   770701  DATE WRITTEN
C   890531  Changed all specific intrinsics to generic.  (WRB)
C   890531  REVISION DATE from Version 3.2
C   891214  Prologue converted to Version 4.0 format.  (BAB)
C   900315  CALLs to XERROR changed to CALLs to XERMSG.  (THJ)
C***END PROLOGUE  DGAMI
      DOUBLE PRECISION a, x, factor, dlngam, dgamit
C***FIRST EXECUTABLE STATEMENT  DGAMI
      IF (a .LE. 0.d0) CALL xermsg('SLATEC', 'DGAMI',
     +   'A MUST BE GT ZERO', 1, 2)
      IF (x .LT. 0.d0) CALL xermsg('SLATEC', 'DGAMI',
     +   'X MUST BE GE ZERO', 2, 2)
C
      dgami = 0.d0
      IF (x.EQ.0.0d0) RETURN
C
C THE ONLY ERROR POSSIBLE IN THE EXPRESSION BELOW IS A FATAL OVERFLOW.
      factor = exp(dlngam(a) + a*log(x))
C
      dgami = factor * dgamit(a, x)
C
      RETURN
      END