Re: [Help-gsl] Gauss hypergeometric

[Brian Gough]

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