HIJING/hijing1_36/gmre.F

   1 * $Id$
   2 C
   3 C
   4 C
   5         FUNCTION GMRE(X)
   6         IMPLICIT REAL*8(A-H,O-Z)
   7         Z=X
   8         IF(X.GT.3.0D0) GO TO 10
   9         Z=X+3.D0
  10 10      GMRE=0.5D0*DLOG(2.D0*3.14159265D0/Z)+Z*DLOG(Z)-Z+DLOG(1.D0
  11      1  +1.D0/12.D0/Z+1.D0/288.D0/Z**2-139.D0/51840.D0/Z**3
  12      1  -571.D0/2488320.D0/Z**4)
  13         IF(Z.EQ.X) GO TO 20
  14         GMRE=GMRE-DLOG(Z-1.D0)-DLOG(Z-2.D0)-DLOG(Z-3.D0)
  15 20      CONTINUE
  16         RETURN
  17         END