Mercurial > octave-nkf

comparison scripts/statistics/distributions/gampdf.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 50
51 k = find ((x > 0) & (a > 0) & (a <= 1) & (b > 0)); 51 k = find ((x > 0) & (a > 0) & (a <= 1) & (b > 0));
52 if (any (k)) 52 if (any (k))
53 if (isscalar(a) && isscalar(b)) 53 if (isscalar(a) && isscalar(b))
54 pdf(k) = (x(k) .^ (a - 1)) ... 54 pdf(k) = (x(k) .^ (a - 1)) ...
55 .* exp(- x(k) ./ b) ./ gamma (a) ./ (b .^ a); 55 .* exp(- x(k) ./ b) ./ gamma (a) ./ (b .^ a);
56 else 56 else
57 pdf(k) = (x(k) .^ (a(k) - 1)) ... 57 pdf(k) = (x(k) .^ (a(k) - 1)) ...
58 .* exp(- x(k) ./ b(k)) ./ gamma (a(k)) ./ (b(k) .^ a(k)); 58 .* exp(- x(k) ./ b(k)) ./ gamma (a(k)) ./ (b(k) .^ a(k));
59 endif 59 endif
60 endif 60 endif
61 61
62 k = find ((x > 0) & (a > 1) & (b > 0)); 62 k = find ((x > 0) & (a > 1) & (b > 0));
63 if (any (k)) 63 if (any (k))
64 if (isscalar(a) && isscalar(b)) 64 if (isscalar(a) && isscalar(b))
65 pdf(k) = exp (- a .* log (b) + (a-1) .* log (x(k)) 65 pdf(k) = exp (- a .* log (b) + (a-1) .* log (x(k))
66 - x(k) ./ b - gammaln (a)); 66 - x(k) ./ b - gammaln (a));
67 else 67 else
68 pdf(k) = exp (- a(k) .* log (b(k)) + (a(k)-1) .* log (x(k)) 68 pdf(k) = exp (- a(k) .* log (b(k)) + (a(k)-1) .* log (x(k))
69 - x(k) ./ b(k) - gammaln (a(k))); 69 - x(k) ./ b(k) - gammaln (a(k)));
70 endif 70 endif
71 endif 71 endif
72 72
73 endfunction 73 endfunction