MICROCERN/gamma.F

   1 *
   2 * $Id$
   3 *
   4 #include "gen/pilot.h"
   5 #if (defined(CERNLIB_DOUBLE))&&(!defined(CERNLIB_IBM))
   6       FUNCTION GAMMA(X)
   7       DOUBLE PRECISION DGAMMA,D
   8
   9       SROUND(D)=D+(D-SNGL(D))
  10       GAMMA=SROUND(DGAMMA(DBLE(X)))
  11       RETURN
  12       END
  13 #endif