Mercurial > hg > octave-lyh

comparison scripts/statistics/distributions/nbinrnd.m @ 10525:3306cfcb856e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Replace constructs like "NaN * one()" with "NaN()" and "Inf * ones ()" with "Inf()"
author David Bateman <dbateman@free.fr>
date Fri, 16 Apr 2010 10:32:07 +0200
parents 16f53d29049f
children 95c3e38098bf
comparison
equal deleted inserted replaced
10524:1c6ff93c025a 10525:3306cfcb856e
75 print_usage (); 75 print_usage ();
76 endif 76 endif
77 77
78 if (isscalar (n) && isscalar (p)) 78 if (isscalar (n) && isscalar (p))
79 if ((n < 1) || (n == Inf) || (n != round (n)) || (p <= 0) || (p > 1)); 79 if ((n < 1) || (n == Inf) || (n != round (n)) || (p <= 0) || (p > 1));
80 rnd = NaN * ones (sz); 80 rnd = NaN (sz);
81 elseif ((n > 0) && (n < Inf) && (n == round (n)) && 81 elseif ((n > 0) && (n < Inf) && (n == round (n)) &&
82 (p > 0) && (p <= 1)) 82 (p > 0) && (p <= 1))
83 rnd = randp ((1 - p) ./ p .* randg (n, sz)); 83 rnd = randp ((1 - p) ./ p .* randg (n, sz));
84 else 84 else
85 rnd = zeros (sz); 85 rnd = zeros (sz);