[Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps

[Marco Caliari]

Follow-up Comment #7, bug #47800 (project octave):

I had almost finished to implement the options scaledupper and scaledlower
when I realized that Matlab supports negative x, while slatec (and thus
Octave) does not. On the other hand, gammainc(-10,10) gives a complex result,
while wolframalpha a real result (Matlab result's real part) and the usual
series definition agrees with wolframalpha. gammainc(-10,10.1) gives the
correct result, in Matlab and as the usual series. Now I think that a plain m
file should be used (like expint, for instance).

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?47800>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/