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

[Hatef Monajemi]

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? 

Bests, 
Hatef 



----- Original Message -----
From: "Hatef Monajemi" <address@hidden> 
To: address@hidden 
Sent: Thursday, August 18, 2011 4:01:29 AM 
Subject: [Bug-gsl] Gauss's hypergeometric function 

Hello, 

I am trying to use the " gsl_sf_hyperg_2F1 ( double a, double b, double c, 
double x ) " function. I noticed that the results are wrong for positive x. Is 
this a known bug? If yes, do you plan to fix it? 

Bests, 
Hatef 

_______________________________________________ 
Bug-gsl mailing list 
address@hidden 
https://lists.gnu.org/mailman/listinfo/bug-gsl