2 DOUBLE PRECISION FUNCTION d9gmit (A, X, ALGAP1, SGNGAM, ALX)
28 DOUBLE PRECISION a, x, algap1, sgngam, alx, ae, aeps, algs, alg2,
29 1 bot,
eps, fk, s, sgng2, t, te,
d1mach, dlngam
40 IF (x .LE. 0.d0)
CALL xermsg (
'SLATEC',
'D9GMIT',
41 +
'X SHOULD BE GT 0', 1, 2)
44 IF (a.LT.0.d0) ma = a - 0.5d0
48 IF (a.LT.(-0.5d0)) ae = aeps
60 CALL xermsg (
'SLATEC',
'D9GMIT',
61 +
'NO CONVERGENCE IN 200 TERMS OF TAYLOR-S SERIES', 2, 2)
63 30
IF (a.GE.(-0.5d0)) algs = -algap1 +
log(s)
64 IF (a.GE.(-0.5d0))
GO TO 60
66 algs = -dlngam(1.d0+aeps) +
log(s)
72 t = x*t/(aeps-(m+1-k))
78 algs = -ma*
log(x) + algs
79 IF (s.EQ.0.d0 .OR. aeps.EQ.0.d0)
GO TO 60
81 sgng2 = sgngam * sign(1.0d0, s)
82 alg2 = -x - algap1 +
log(
abs(s))
84 IF (alg2.GT.bot) d9gmit = sgng2 * exp(alg2)
85 IF (algs.GT.bot) d9gmit = d9gmit + exp(algs)
OCTAVE_EXPORT octave_value_list or N dimensional array whose elements are all equal to the base of natural logarithms The constant ex $e satisfies the equation log(e)
double precision function d1mach(i)