]>
Commit | Line | Data |
---|---|---|
ef94df36 | 1 | *CMZ : 31/08/93 17.51.25 by Jonathan Butterworth |
2 | *-- Author : | |
3 | DOUBLE PRECISION FUNCTION JMFACT( N ) | |
4 | C -------------------------------------------------------------------------- | |
5 | C Factorial of N, given as double precision to allow for very large numbers. | |
6 | C -------------------------------------------------------------------------- | |
7 | IMPLICIT NONE | |
8 | INTEGER N, I | |
9 | DOUBLE PRECISION TOTAL | |
10 | ||
11 | TOTAL = 1.D0 | |
12 | DO I = 2, N | |
13 | TOTAL = TOTAL*DBLE(I) | |
14 | ENDDO | |
15 | JMFACT = TOTAL | |
16 | RETURN | |
17 | END |