[Bug-gsl] Bug in gsl_sf_hyperg_2F1

[Andrew Benson]

Under GSL 1.14 (compiled with either gcc version 4.4.4 or gcc version 4.1.2) 
the gsl_sf_hyperg_2F1 function crashes with the following arguments:

gsl_sf_hyperg_2F1(-10.34,2.05,3.05,0.1725)

The expected answer (computed using Maple for example) is 0.3104735522

I've attached an example case:

$ g++ test.cpp -o test.x -lgsl -lgslcblas
$ test.x
gsl: hyperg_2F1.c:755: ERROR: error
Default GSL error handler invoked.
Abort

The crash seems to occur only when the first argument is less than or equal to 
-10. For example, 

gsl_sf_hyperg_2F1(-10.0,2.05,3.05,0.1725)

crashes, but

gsl_sf_hyperg_2F1(-9.99999999999,2.05,3.05,0.1725)

evaluates correctly to 0.321419346301974773.

test.cpp
Description: Text Data