* 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. -- Rafael