33 alnsml =
log(r1mach(1))
38 xmin = xmin - xmin*((xmin+0.5)*xln - xmin - 0.2258 + alnsml)
40 IF (
abs(xmin-xold).LT.0.005) go to 20
42 CALL
xermsg(
'SLATEC',
'GAMLIM',
'UNABLE TO FIND XMIN', 1, 2)
44 20 xmin = -xmin + 0.01
46 alnbig =
log(r1mach(2))
51 xmax = xmax - xmax*((xmax-0.5)*xln - xmax + 0.9189 - alnbig)
53 IF (
abs(xmax-xold).LT.0.005) go to 40
55 CALL
xermsg(
'SLATEC',
'GAMLIM',
'UNABLE TO FIND XMAX', 2, 2)
58 xmin =
max(xmin, -xmax+1.)
octave_value log(void) const
charNDArray max(char d, const charNDArray &m)
subroutine xermsg(LIBRAR, SUBROU, MESSG, NERR, LEVEL)
subroutine gamlim(XMIN, XMAX)