]>
Commit | Line | Data |
---|---|---|
e74335a4 | 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 |