Re: [Bug-gsl] possible bug in gsl_sf_hyperg_U

[Brian Gough]

At Thu, 3 Apr 2008 11:22:01 +0200,
Stefan Gerlach wrote:
> I have a problem using gsl_sf_hyperg_U with certain arguments.
> running this example
> 
>        for(z=3.0;z>0.0;z-=0.1)
>                 printf("%g %g\n",z,gsl_sf_hyperg_U(-2.0,4.0,z));
> 
> i get 
> 
>       ...
>       2.2 2.84
>       2.1 3.41
>       2 4
>       gsl: psi.c:381: ERROR: domain error
>       Default GSL error handler invoked.
>       Aborted
> 

Thanks for the bug report.  I have logged it in the bug database and
added a test case for it to the regression tests.

-- 
Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/