Re: [Bug-gsl] Gauss's hypergeometric function

[Brian Gough]

At Sat, 20 Aug 2011 04:58:42 -0700 (PDT),
Hatef Monajemi wrote:
> 
> More specifically, I am trying the following: 
> ====================== 
> for (size_t i=1; i<5 ; ++i){ 
> double a = gsl_sf_hyperg_2F1 (i-0.5, -0.5, i+1, 0.9); 
> std::cout<< "a:"<< a <<std::endl; 
> } 
> 
> ====================== 
> Note that when 'z < .5', I get the right results. for .5<z<1.0 which is the 
> case for me, the results are different than those of Matlab. I am assuming 
> Matlab implementation is correct. I do not have access to Mathematica though. 
> It seems for 'c-a-b = m' when m is an integer and 0.5<z<1.0 the code is not 
> working properly. Any comment? 
> 

Thanks for the bug report.  I agree there is a problem with those
values.  I have added a test case for the problem, but I don't have a
fix.

-- 
Brian Gough