ASA245
Logarithm of Gamma Function

ASA245 is a FORTRAN77 library which evaluates the logarithm of the gamma function, by Allan Mcleod;

ASA245 is Applied Statistics Algorithm 245. Source code for many Applied Statistics Algorithms is available through STATLIB.

Languages:

ASA245 is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

FN, a FORTRAN77 library which approximates elementary and special functions using Chebyshev polynomials, by Wayne Fullerton.

GSL, a C++ library which includes many routines for evaluating special functions.

SLATEC, a FORTRAN90 library which evaluates many special functions.

SPECFUN, a FORTRAN77 library which evaluates certain special functions using fitted data.

TEST_VALUES, a FORTRAN77 library which contains sample values for a number of distributions.

TOMS291, a FORTRAN77 library which evaluates the logarithm of the Gamma function.

TOMS708, a FORTRAN77 library which evaluates several special functions, including the logarithm of the Gamma function.

Author:

Original FORTRAN77 version by Allan Mcleod; This FORTRAN77 version by John Burkardt.

Reference:

1. Cornelius Lanczos,
   A precision approximation of the gamma function,
   SIAM Journal on Numerical Analysis, B,
   Volume 1, 1964, pages 86-96.
2. Allan Mcleod,
   Algorithm AS 245: A Robust and Reliable Algorithm for the Logarithm of the Gamma Function,
   Applied Statistics,
   Volume 38, Number 2, 1989, pages 397-402.

Source Code:

• asa245.f, the source code.
• asa245.sh, commands to compile the source code.

Examples and Tests:

• asa245_prb.f, a sample calling program.
• asa245_prb.sh, commands to compile and run the sample program.
• asa245_prb_output.txt, the output file.

List of Routines:

• ALNGAM computes the logarithm of the gamma function.
• GAMMA_LOG_VALUES returns some values of the Log Gamma function.
• LNGAMMA computes Log(Gamma(X)) using a Lanczos approximation.
• TIMESTAMP prints out the current YMDHMS date as a timestamp.

You can go up one level to the FORTRAN77 source codes.