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

[Marco Caliari]

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

I release my contribution under the GPL, too.

The only functionality missing with respect to the current Octave
implementation (and Matlab as well) is the possibility to work with a and x
vectors. I do not see an easy way to vectorize: in expint.m, for instance,
there is a simple 


for t = 1:numel (x)
...
end


    _______________________________________________________

Reply to this item at:

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

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