GNOME Bugzilla – Bug 96739
GAMMAINV
Last modified: 2004-12-22 21:47:04 UTC
GAMMAINV for p close to 1 returns an error. A1=gammainv(b1,c1,d1) B1=1-1e-13 C1=1 D1=1 (B2=1-B1 just to prove that B1 has not become 1.) A1=#NUM! for me. We can do better.
Hmm... This comes from R code (qgamma in mathfunc.c) which specifically cuts off at 1e-12 from the endpoint.
[From Ethan Glasser-Camp <glasse@rpi.edu>] =gammainv(.3,.6,64) --> #NUM! =gammainv(.3,64,.6) --> #NUM! What gives?
The latter stupidity has been fixed in cvs. (1.1 and 1.0.)
Fixed in 1.1 a while ago.