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

[Marco Caliari]

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

First of all, the trick with


jj = abs ((Deltaj - 1) ./ y) > eps;


is missing in ngammainc_l. Now included.

Second, ngammainc_l is used for x>0, therefore always real results. What about
to write an oct file just for ngammainc_l? it should not be difficult.


(file #37335)
    _______________________________________________________

Additional Item Attachment:

File name: gammainc.m                     Size:10 KB


    _______________________________________________________

Reply to this item at:

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

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