|Author(s): K.S. Kölbig||Library: MATHLIB|
|Language: Fortran||Revised: 15.03.1993|
Function subprogram GAMDIS calculates the gamma distribution function (incomplete gamma function)
for real arguments and a > 0.
User Entry Name: GAMDIS
Files Referenced: Unit 6
External References: GAMMA, ALGAMA, MTLMTR, ABEND
In any arithmetic expression,
GAMDIS(X,A) has the value ,
The method is described in Ref. 1.
Approximately six digits are correct.
Error G106.1: or .
Error G106.2: Difficulties of convergence (unlikely).
The function value is set equal to zero, and a message is written on Unit 6, unless subroutine MTLSET (N002) has been called.
This subprogram is based on a Fortran program for the incomplete gamma functions published in Ref. 2.