Re: [Bug-gsl] Bug in normalized incomplete Gamma function/upper tail Gamma dist'n probs.

[Brian Gough]

Brian  Williams writes:
 > Version number of GSL:  1.6
 > Hardware & OS:  Dell PowerEdge 6850, Intel Xeon processor at 3.0GHz/8MB
 > Cache, Linux 2.4.21-27.EL
 > Compiler:  gcc 3.2.3 20030502 (Red Hat Linux 3.2.3-52)
 >    gcc -o test.out -Wall -lgsl -lgslcblas -lm test.c
 > Description:  For some inputs, answers are returned that are incorrect. 
 > Here is a test program (test.c):

Thanks for the bug report.  This will be fixed in the next release,
which should be available in the next week or so.  The cause was an
instability in a continued fraction, you can find the corrected code 
in cvs at http://sources.redhat.com/gsl/devel.html (specfunc/gamma_inc.c).

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/