Mercurial > hg > octave-nkf

diff scripts/deprecated/lognormal_rnd.m @ 5692:fe226f54d259

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2006-03-17 17:02:32 by jwe]
author jwe
date Fri, 17 Mar 2006 17:02:32 +0000
parents 2a16423e4aa0
children 93c65f2a5668
line wrap: on
line diff
--- a/scripts/deprecated/lognormal_rnd.m
+++ b/scripts/deprecated/lognormal_rnd.m
@@ -34,6 +34,22 @@
 
 function rnd = lognormal_rnd (varargin)
 
- rnd =  lognrnd (varargin{:});
+  if (nargin > 1)
+    a = varargin{2};
+    idx = a >= 0;
+    a(idx) = log (a(idx));
+    a(!idx) = NaN;
+    varargin{2} = a;
+  endif
+
+  if (nargin > 2)
+    v = varargin{3};
+    idx = v >= 0;
+    v(idx) = sqrt (v(idx));
+    v(!idx) = NaN;
+    varargin{3} = v;
+  endif
+
+ rnd = lognrnd (varargin{:});
 
 endfunction