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