Re: gammai in 2.0.8

[Ted Harding]

> > I have found that Kurt Hornik's gammai(a,x) function seg-faults
> > when run in my octave-2.0.8.
> 
> It seems to work well in 2.0.9 (Linux  1.2.13  i586):
> 
> > Octave, version 2.0.9 (i586-pc-linux-gnu).

Thanks, Mario. If it's any clue, I have since found that bessel
gives a FORTRAN exception:

    octave-2:1> bessel(1.3,5)
     ***MESSAGE FROM ROUTINE DGAMMA IN LIBRARY SLATEC.
     ***FATAL ERROR, PROGRAM ABORTED.
     *  X IS 0
     *  ERROR NUMBER = 4
     *
     ***END OF MESSAGE

     ***JOB ABORT DUE TO FATAL ERROR.
    error: exception encountered in Fortran subroutine xdgamma

However, my mileage varies:

    octave-2:1> bessel(.5,5)
    ans = -0.34217
    octave-2:2> bessel(3.5,2.5)
    error: Segmentation fault -- stopping myself...
    attempting to save variables to `octave-core'...
    save to `octave-core' complete

Ted.                                    address@hidden