At Thu, 24 Jun 2010 21:11:53 +0200,
Heiko Bauke wrote:
> I would like to point out that the GSL functions gsl_sf_hyperg_2F1_e and
> gsl_sf_hyperg_2F1 which compute the Gauss hypergeometric function for
> -1<=x<1 may be extended with little effort to arguments in the interval
> -oo<x<1 by using relations (32) and/or (33) in [1] by something like
> this:
>
> double hyperg_2F1(double a, double b, double c, double x) {
> ..
> }
>
> I cannot say much about the accuracy and the stability of this approach.
> I think, however, it is reasonable to extend gsl_sf_hyperg_2F1_e and
> gsl_sf_hyperg_2F1 in this way. It would make these functions more
> general and useful.
Thanks for the suggestion. I have logged it in the bug tracker at
https://savannah.gnu.org/bugs/index.php?30324 - if you try out the
method you can post any comments about the accuracy there.
--
Brian Gough
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
