gmre.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file gmre.f

C
C
C
        FUNCTION gmre(X)
        IMPLICIT REAL*8(a-h,o-z)
        z=x
        IF(x.GT.3.0d0) go to 10
        z=x+3.d0
10      gmre=0.5d0*dlog(2.d0*3.14159265d0/z)+z*dlog(z)-z+dlog(1.d0
     1  +1.d0/12.d0/z+1.d0/288.d0/z**2-139.d0/51840.d0/z**3
     1  -571.d0/2488320.d0/z**4)
        IF(z.EQ.x) go to 20
        gmre=gmre-dlog(z-1.d0)-dlog(z-2.d0)-dlog(z-3.d0)
20      CONTINUE
        RETURN
        END