Mercurial > octave
diff scripts/specfun/betainc.m @ 30400:168da23530b4 stable
maint: strip trailing spaces from code base.
author | Rik <rik@octave.org> |
---|---|
date | Mon, 29 Nov 2021 16:43:46 -0800 |
parents | 01de0045b2e3 |
children | 796f54d4ddbf |
line wrap: on
line diff
--- a/scripts/specfun/betainc.m Mon Nov 29 15:58:51 2021 -0800 +++ b/scripts/specfun/betainc.m Mon Nov 29 16:43:46 2021 -0800 @@ -136,7 +136,7 @@ non_trivial = a_not_one & b_not_one; a_one &= b_not_one; b_one &= a_not_one; - + if (strcmpi (tail, "lower")) y(a_b_one) = x(a_b_one); y(a_one) = 1 - (1 - x(a_one)) .^ b(a_one); @@ -146,13 +146,13 @@ y(a_one) = (1 - x(a_one)) .^ b(a_one); y(b_one) = 1 - x(b_one) .^ a(b_one); endif - + ## Non-Trivial cases ## In the following, we use the fact that the continued fraction Octave uses ## is more efficient when x <= a / (a + b). Moreover, to compute the upper ## version, which is defined as I_x(a,b,"upper") = 1 - I_x(a,b) we use the ## property I_x(a,b) + I_(1-x) (b,a) = 1. - + x = x(non_trivial); a = a(non_trivial); b = b(non_trivial); @@ -182,7 +182,7 @@ + (gammaln (a + b) - gammaln (a) - gammaln (b)) + log (f); y_nt = real (exp (y_nt)); y_nt(fflag) = 1 - y_nt(fflag); - + y(non_trivial) = y_nt; ## Restore original shape