Mercurial > octave-nkf

comparison scripts/statistics/distributions/fpdf.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 a1dbe9d80eee
children 1740012184f9
comparison
equal deleted inserted replaced
10548:479536c5bb10 10549:95c3e38098bf
50 k = find ((x > 0) & (x < Inf) & (m > 0) & (n > 0)); 50 k = find ((x > 0) & (x < Inf) & (m > 0) & (n > 0));
51 if (any (k)) 51 if (any (k))
52 if (isscalar (m) && isscalar (n)) 52 if (isscalar (m) && isscalar (n))
53 tmp = m / n * x(k); 53 tmp = m / n * x(k);
54 pdf(k) = (exp ((m / 2 - 1) .* log (tmp) 54 pdf(k) = (exp ((m / 2 - 1) .* log (tmp)
55 - ((m + n) / 2) .* log (1 + tmp)) 55 - ((m + n) / 2) .* log (1 + tmp))
56 .* (m / n) ./ beta (m / 2, n / 2)); 56 .* (m / n) ./ beta (m / 2, n / 2));
57 else 57 else
58 tmp = m(k) .* x(k) ./ n(k); 58 tmp = m(k) .* x(k) ./ n(k);
59 pdf(k) = (exp ((m(k) / 2 - 1) .* log (tmp) 59 pdf(k) = (exp ((m(k) / 2 - 1) .* log (tmp)
60 - ((m(k) + n(k)) / 2) .* log (1 + tmp)) 60 - ((m(k) + n(k)) / 2) .* log (1 + tmp))
61 .* (m(k) ./ n(k)) ./ beta (m(k) / 2, n(k) / 2)); 61 .* (m(k) ./ n(k)) ./ beta (m(k) / 2, n(k) / 2));
62 endif 62 endif
63 endif 63 endif
64 64
65 endfunction 65 endfunction