GAMMA(III) 5/15/74 GAMMA(III)
NAME
gamma - log gamma function
SYNOPSIS
jsr pc,gamma
double gamma(x)
double x;
DESCRIPTION
If x is passed (in fr0) gamma returns ln |G(x)| (in fr0).
The sign of G(x) is returned in the external integer
signgam. The following C program might be used to calculate
G:
y = gamma(x);
if (y > 88.)
error( );
y = exp(y);
if(signgam)
y = -y;
DIAGNOSTICS
The c-bit is set on negative integral arguments and the max-
imum value is returned. There is no error return for C pro-
grams.
BUGS
No error return from C.
- 1 -