2 FUNCTION r9gmit (A, X, ALGAP1, SGNGAM, ALX)
28 DATA eps, bot / 2*0.0 /
33 IF (x .LE. 0.0)
CALL xermsg (
'SLATEC',
'R9GMIT',
34 +
'X SHOULD BE GT 0', 1, 2)
37 IF (a.LT.0.0) ma = a - 0.5
41 IF (a.LT.(-0.5)) ae = aeps
53 CALL xermsg (
'SLATEC',
'R9GMIT',
54 +
'NO CONVERGENCE IN 200 TERMS OF TAYLOR-S SERIES', 2, 2)
56 30
IF (a.GE.(-0.5)) algs = -algap1 +
log(s)
57 IF (a.GE.(-0.5))
GO TO 60
59 algs = -alngam(1.0+aeps) +
log(s)
71 algs = -ma*
log(x) + algs
72 IF (s.EQ.0.0 .OR. aeps.EQ.0.0)
GO TO 60
74 sgng2 = sgngam*sign(1.0,s)
75 alg2 = -x - algap1 +
log(
abs(s))
77 IF (alg2.GT.bot) r9gmit = sgng2*exp(alg2)
78 IF (algs.GT.bot) r9gmit = r9gmit + 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)