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" <[email protected]>
To: [email protected]
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
[email protected]
https://lists.gnu.org/mailman/listinfo/bug-gsl
_______________________________________________
Bug-gsl mailing list
[email protected]
https://lists.gnu.org/mailman/listinfo/bug-gsl
