Mercurial > octave
changeset 25022:fdff879d8ae8
maint: merge stable to default.
author | Rik <rik@octave.org> |
---|---|
date | Tue, 27 Mar 2018 11:38:05 -0700 |
parents | d0f9826f677d (current diff) 592db8745f26 (diff) |
children | f2f813a55eec |
files | |
diffstat | 1 files changed, 2 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/scripts/specfun/betainc.m Mon Mar 26 21:17:17 2018 -0700 +++ b/scripts/specfun/betainc.m Tue Mar 27 11:38:05 2018 -0700 @@ -147,8 +147,8 @@ f = __betainc__ (x, a, b); ## Divide continued fraction by B(a,b) / (x^a * (1-x)^b) to obtain I_x(a,b). - y = a .* log (x) + b .* log1p (-x) + gammaln (a + b) - ... - gammaln (a) - gammaln (b) + log (f); + y = a .* log (x) + b .* log1p (-x) ... + + (gammaln (a + b) - gammaln (a) - gammaln (b)) + log (f); y = real (exp (y)); y(fflag) = 1 - y(fflag);