Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:15:42

0001 C
0002 C
0003 C
0004         FUNCTION GMRE(X)
0005         IMPLICIT REAL*8(A-H,O-Z)
0006         Z=X
0007         IF(X.GT.3.0D0) GO TO 10
0008         Z=X+3.D0
0009 10      GMRE=0.5D0*DLOG(2.D0*3.14159265D0/Z)+Z*DLOG(Z)-Z+DLOG(1.D0
0010      1  +1.D0/12.D0/Z+1.D0/288.D0/Z**2-139.D0/51840.D0/Z**3
0011      1  -571.D0/2488320.D0/Z**4)
0012         IF(Z.EQ.X) GO TO 20
0013         GMRE=GMRE-DLOG(Z-1.D0)-DLOG(Z-2.D0)-DLOG(Z-3.D0)
0014 20      CONTINUE
0015         RETURN
0016         END