* difficult region for gsl_sf_hyperg_1F1
@ 2009-12-15 18:36 Brian Gough
0 siblings, 0 replies; only message in thread
From: Brian Gough @ 2009-12-15 18:36 UTC (permalink / raw)
To: GSL Discuss Mailing List
We have gsl_sf_hyperg1F1(-37.8, 2.01, 103.58) => -3.16661e+14 +/-
3.31+19 from the series. Correct answer is -6.21927211009e17(*).
Does anyone know a better algorithm for this region?
http://savannah.gnu.org/bugs/?28267
(*) In pari:
poch(a,x) = { gamma(a+x)/gamma(a) }
t(a,b,x,k) = { (poch(a,k)/poch(b,k)) * (x^k)/(k!) }
suminf(k=0,t(-37.8,2.01,103.58,k))
^ permalink raw reply [flat|nested] only message in thread
only message in thread, other threads:[~2009-12-15 18:36 UTC | newest]
Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2009-12-15 18:36 difficult region for gsl_sf_hyperg_1F1 Brian Gough
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).