Re: [Bug-gsl] bug in hypergeometric U

[Krzysztof Pachucki]

Dear Brian,
the example follows, gsl 1.3:

main()
{
  double t,gsl_sf_hyperg_U_int();
  t =gsl_sf_hyperg_U_int(1,4,0.5); 
  printf("16.8f \n",t);
           }

gcc bug.c -lgsl -lgslcblas
It is the case when U functiopn is a finite series expansion.

Krzysztof

On Mon, 14 Apr 2003, Brian Gough wrote:

> Krzysztof Pachucki writes:
>  > Dear Developers, gsl_sf_hyperg_U(1.0_16, 4.0_16, double x)
>  > crashes (integer version as well) with error message
>  > gsl: gamma.c:774: ERROR: error
>  > any hints to overcome it are appreciated
>  > 
> 
> Please send an example program that reproduces the bug.
> Thanks.
> 
> 

-- 
Krzysztof Pachucki

Institute of Theoretical Physics
Warsaw University
Hoza 69, 00-681 Warsaw, Poland
Tel: + (48-22) 5532246
Fax: + (48-22) 6219475