2 DOUBLE PRECISION FUNCTION d9lgit (A, X, ALGAP1)
29 DOUBLE PRECISION a, x, algap1, ax, a1x,
eps, fk, hstar, p, r, s,
32 SAVE eps, sqeps, first
41 IF (x .LE. 0.d0 .OR. a .LT. x)
CALL xermsg (
'SLATEC',
'D9LGIT',
42 +
'X SHOULD BE GT 0.0 AND LE A', 2, 2)
52 r = t/((ax+fk)*(a1x+fk)-t)
55 IF (
abs(p).LT.
eps*s)
GO TO 30
57 CALL xermsg (
'SLATEC',
'D9LGIT',
58 +
'NO CONVERGENCE IN 200 TERMS OF CONTINUED FRACTION', 3, 2)
60 30 hstar = 1.0d0 - x*s/a1x
61 IF (hstar .LT. sqeps)
CALL xermsg (
'SLATEC',
'D9LGIT',
62 +
'RESULT LESS THAN HALF PRECISION', 1, 1)
64 d9lgit = -x - algap1 -
log(hstar)
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)