Mercurial > octave-nkf

comparison scripts/statistics/distributions/frnd.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 3306cfcb856e
children be55736a0783
comparison
equal deleted inserted replaced
10548:479536c5bb10 10549:95c3e38098bf
36 36
37 if (nargin > 1) 37 if (nargin > 1)
38 if (!isscalar(m) || !isscalar(n)) 38 if (!isscalar(m) || !isscalar(n))
39 [retval, m, n] = common_size (m, n); 39 [retval, m, n] = common_size (m, n);
40 if (retval > 0) 40 if (retval > 0)
41 error ("frnd: m and n must be of common size or scalar"); 41 error ("frnd: m and n must be of common size or scalar");
42 endif 42 endif
43 endif 43 endif
44 endif 44 endif
45 45
46 46
52 error ("frnd: c must be a positive integer"); 52 error ("frnd: c must be a positive integer");
53 endif 53 endif
54 sz = [r, c]; 54 sz = [r, c];
55 55
56 if (any (size (m) != 1) && 56 if (any (size (m) != 1) &&
57 ((length (size (m)) != length (sz)) || any (size (m) != sz))) 57 ((length (size (m)) != length (sz)) || any (size (m) != sz)))
58 error ("frnd: m and n must be scalar or of size [r,c]"); 58 error ("frnd: m and n must be scalar or of size [r,c]");
59 endif 59 endif
60 elseif (nargin == 3) 60 elseif (nargin == 3)
61 if (isscalar (r) && (r > 0)) 61 if (isscalar (r) && (r > 0))
62 sz = [r, r]; 62 sz = [r, r];
65 else 65 else
66 error ("frnd: r must be a positive integer or vector"); 66 error ("frnd: r must be a positive integer or vector");
67 endif 67 endif
68 68
69 if (any (size (m) != 1) && 69 if (any (size (m) != 1) &&
70 ((length (size (m)) != length (sz)) || any (size (m) != sz))) 70 ((length (size (m)) != length (sz)) || any (size (m) != sz)))
71 error ("frnd: m and n must be scalar or of size sz"); 71 error ("frnd: m and n must be scalar or of size sz");
72 endif 72 endif
73 elseif (nargin == 2) 73 elseif (nargin == 2)
74 sz = size(a); 74 sz = size(a);
75 else 75 else
78 78
79 79
80 if (isscalar (m) && isscalar (n)) 80 if (isscalar (m) && isscalar (n))
81 if (isinf (m) || isinf (n)) 81 if (isinf (m) || isinf (n))
82 if (isinf (m)) 82 if (isinf (m))
83 rnd = ones (sz); 83 rnd = ones (sz);
84 else 84 else
85 rnd = 2 ./ m .* randg(m / 2, sz); 85 rnd = 2 ./ m .* randg(m / 2, sz);
86 endif 86 endif
87 if (! isinf (n)) 87 if (! isinf (n))
88 rnd = 0.5 .* n .* rnd ./ randg (n / 2, sz); 88 rnd = 0.5 .* n .* rnd ./ randg (n / 2, sz);
89 endif 89 endif
90 elseif ((m > 0) && (m < Inf) && (n > 0) && (n < Inf)) 90 elseif ((m > 0) && (m < Inf) && (n > 0) && (n < Inf))
91 rnd = n ./ m .* randg (m / 2, sz) ./ randg (n / 2, sz); 91 rnd = n ./ m .* randg (m / 2, sz) ./ randg (n / 2, sz);
92 else 92 else
93 rnd = NaN (sz); 93 rnd = NaN (sz);