[Octave-bug-tracker] [bug #62329] Incorrect value (regression) for betainc

[Rik]

Follow-up Comment #2, bug #62329 (project octave):

This is certainly caused by the change introduced for bug #60682.  Adding
Michele Ginesi to the CC list.

The code in question is the fast path code when a==1.


I(a_one) = 1 - (1 - x(a_one)) .^ b(a_one);


The specific issue is 


(1 - x(a_one))


x is 2e-20 which is less than eps(1) so in the subtraction the result is just
a perfect 1.


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?62329>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/