|
| From: | Nir Krakauer |
| Subject: | [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps |
| Date: | Mon, 9 May 2016 16:18:20 +0000 (UTC) |
| User-agent: | Mozilla/5.0 (Macintosh; Intel Mac OS X 10_11_3) AppleWebKit/601.4.4 (KHTML, like Gecko) Version/9.0.3 Safari/601.4.4 |
Follow-up Comment #15, bug #47800 (project octave):
Thanks, Marco.
I added the exponential sum algorithm for small integer a and small x, which
gives near-exact results.
Also, I found a more recent (2012) discussion involving Temme at
http://oai.cwi.nl/oai/asset/20433/20433B.pdf Their Fortran code is posted at
http://personales.unican.es/gila/incgam.zip , unfortunately with no license.
(file #37114)
_______________________________________________________
Additional Item Attachment:
File name: ngammainc.m Size:5 KB
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?47800>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
| [Prev in Thread] | Current Thread | [Next in Thread] |