[Help-gsl] Computation of gsl_sf_hyperg_1F1 fails for some values

[Daniel Schury]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hi,

I've got some problems with the usage of the above mentioned function.

My usage is
"double part22 = gsl_sf_hyperg_1F1(-p/2.0, 1.0/2.0,
- -1.0/2.0*((E1-x)/dE)^2.0);"
This works for most values, but fails for some. With p = 2.09, E1 = 8, x
= 26.8 and dE = 0.5 my program crashes.

I'm using the library shipped with Origin 9.0 in C, Windows XP, 32 bit.

I'm not completely into the mathematics about the hypergeometric
function, but mathematica for example can compute this function with the
same values, so for me it seems not like a mathematical, but a
programming problem.

Can somebody confirm this problem maybe on a different platform or has
an idea how to trace the problem?

Many thanks in advance.

- -- 
Freundliche Grüße
Daniel Schury
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)
Comment: Using GnuPG with undefined - http://www.enigmail.net/

iEYEARECAAYFAlGMwVwACgkQk7XSfqVgB0oIsQCfU42V6cW45UbT+o20pdehy+Yc
70QAoLh3P/AHi3KAPXYlLM/PYU71FwAs
=8vxP
-----END PGP SIGNATURE-----