    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 -