Mercurial > octave-nkf

diff scripts/statistics/distributions/frnd.m @ 6339:20c48710b2c7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-02-22 18:22:24 by dbateman] ChangeLog
author dbateman
date Thu, 22 Feb 2007 18:22:24 +0000
parents 34f96dd5441b
children 6304d9ea0a30
line wrap: on
line diff
--- a/scripts/statistics/distributions/frnd.m	Thu Feb 22 07:07:20 2007 +0000
+++ b/scripts/statistics/distributions/frnd.m	Thu Feb 22 18:22:24 2007 +0000
@@ -80,7 +80,7 @@
 
   if (isscalar (m) && isscalar (n))
     if ((m > 0) && (m < Inf) && (n > 0) && (n < Inf))
-      rnd =  finv (rand (sz), m, n);
+      rnd = n ./ m .* randg(m/2,sz) ./ randg(n/2,sz);
     else
       rnd = NaN * ones (sz);
     endif
@@ -96,7 +96,7 @@
     k = find ((m > 0) & (m < Inf) &
               (n > 0) & (n < Inf));
     if (any (k))
-      rnd(k) = finv (rand (size (k)), m(k), n(k));
+      rnd(k) = n(k) ./ m(k) .* randg(m(k)./2,size(k)) ./ randg(n(k)./2,size(k));
     endif
   endif