Hi, 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
AFAIK the hypergeometricU is well defined for z> 0.0.
I get the same error with gsl-1.11 on openSUSE 10.X i586 and x86_64.
--
Stefan Gerlach
pgp1tGixUuGvr.pgp
Description: PGP signature
_______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
