Mercurial > octave-nkf

comparison scripts/miscellaneous/bincoeff.m @ 10549:95c3e38098bf

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Untabify .m scripts
author Rik <code@nomad.inbox5.com>
date Fri, 23 Apr 2010 11:28:50 -0700
parents f0c3d3fc4903
children 693e22af08ae
comparison
equal deleted inserted replaced
10548:479536c5bb10 10549:95c3e38098bf
78 ind = (k == 0); 78 ind = (k == 0);
79 b(ind) = 1; 79 b(ind) = 1;
80 80
81 ind = ((k > 0) & ((n == real (round (n))) & (n < 0))); 81 ind = ((k > 0) & ((n == real (round (n))) & (n < 0)));
82 b(ind) = (-1) .^ k(ind) .* exp (gammaln (abs (n(ind)) + k(ind)) 82 b(ind) = (-1) .^ k(ind) .* exp (gammaln (abs (n(ind)) + k(ind))
83 - gammaln (k(ind) + 1) 83 - gammaln (k(ind) + 1)
84 - gammaln (abs (n(ind)))); 84 - gammaln (abs (n(ind))));
85 85
86 ind = ((k > 0) & (n >= k)); 86 ind = ((k > 0) & (n >= k));
87 b(ind) = exp (gammaln (n(ind) + 1) 87 b(ind) = exp (gammaln (n(ind) + 1)
88 - gammaln (k(ind) + 1) 88 - gammaln (k(ind) + 1)
89 - gammaln (n(ind) - k(ind) + 1)); 89 - gammaln (n(ind) - k(ind) + 1));
90 90
91 ind = ((k > 0) & ((n != real (round (n))) & (n < k))); 91 ind = ((k > 0) & ((n != real (round (n))) & (n < k)));
92 b(ind) = (1/pi) * exp (gammaln (n(ind) + 1) 92 b(ind) = (1/pi) * exp (gammaln (n(ind) + 1)
93 - gammaln (k(ind) + 1) 93 - gammaln (k(ind) + 1)
94 + gammaln (k(ind) - n(ind)) 94 + gammaln (k(ind) - n(ind))
95 + log (sin (pi * (n(ind) - k(ind) + 1)))); 95 + log (sin (pi * (n(ind) - k(ind) + 1))));
96 96
97 ## Clean up rounding errors. 97 ## Clean up rounding errors.
98 ind = (n == round (n)); 98 ind = (n == round (n));
99 b(ind) = round (b(ind)); 99 b(ind) = round (b(ind));
100 100