* possible bug in gsl_sf_hyperg_U
@ 2008-04-03 12:45 Stefan Gerlach
0 siblings, 0 replies; only message in thread
From: Stefan Gerlach @ 2008-04-03 12:45 UTC (permalink / raw)
To: gsl-discuss
[-- Attachment #1: Type: text/plain, Size: 452 bytes --]
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. Any ideas?
--
Stefan Gerlach
[-- Attachment #2: Type: application/pgp-signature, Size: 189 bytes --]
^ permalink raw reply [flat|nested] only message in thread
only message in thread, other threads:[~2008-04-03 12:45 UTC | newest]
Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2008-04-03 12:45 possible bug in gsl_sf_hyperg_U Stefan Gerlach
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).