[Bug-gsl] bug on gsl_sf_hyperg_1F1

[Weibin Li]

There seems a bug in the function “gsl_sf_hyperg_1F1" in version 1.16.

See the attached test program, in which the generated value is correct when 
ri=1000 while is wrong when ri is large (ri=4000 in the program).
The result was compared with Mathematica. 

When I tried version 1.9, which somehow gives the correct result.
However I am not sure how this happens in
This message and any attachment are intended solely for the addressee and may 
contain confidential information. If you have received this message in error, 
please send it back to me, and immediately delete it.   Please do not use, copy 
or disclose the information contained in this message or in any attachment.  
Any views or opinions expressed by the author of this email do not necessarily 
reflect the views of the University of Nottingham.



This message has been checked for viruses but the contents of an attachment

may still contain software viruses which could damage your computer system, you 
are advised to perform your own checks. Email communications with the 
University of Nottingham may be monitored as permitted by UK legislation.

gsl_hyperg.c
Description: Binary data

 other versions.

I used a Macbook pro with Yosemite.

best

Weibin