On 29-Feb-2008, Rafael Laboissiere wrote: | * Rafael Laboissiere <[EMAIL PROTECTED]> [2008-02-24 18:23]: | | > * Ólafur Jens Sigurðsson <[EMAIL PROTECTED]> [2008-02-24 12:34]: | > | > > I just tried some values and the limit seems to be poisscdf(604,604), | > > after that the function fails, that is poisscdf(605,605) fails and | > > poisscdf(604,604) succeeds. | > | > It fails for smaller values than that, like poisscdf (599, 600). The | > culprit is the gammainc function, which ultimately calls the function | > defined in d9lgit.f. | | Just for the record, the gamma_inc_P function in the gsl octave-forge | package works fine in the range for which gammainc from d9lgit.f fails. For | instance: | | octave:1> gamma_inc_P (1e5,1e5) | ans = 0.50042 | | John, is there any chance of integrating GSL directly in Octave? At least | making the compilation in files like gammainc.cc conditional on the presence | of libgsl and, if not, falling back to Slatec.
Yes, but maybe not until after 3.1 is released unless someone shows up with a patch. jwe
