Mercurial > hg > octave-nkf

diff scripts/control/base/lsim.m @ 4174:aabb8307021f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2002-11-13 00:45:22 by jwe]
author jwe
date Wed, 13 Nov 2002 00:45:23 +0000
parents 0d411821682c
children bdbee5282954
line wrap: on
line diff
--- a/scripts/control/base/lsim.m
+++ b/scripts/control/base/lsim.m
@@ -73,10 +73,13 @@
   t(2)-t(1);
   u=u';
   n = max(size(t));
+
   for ii = 1:(n-1)
 
     ## check if step size changed
-    if (t(ii+1) - t(ii) != Ts)
+    ## XXX FIXME XXX -- this is probably not the best test, but it is
+    ## better than a test for exact equality.
+    if (abs (t(ii+1) - t(ii) - Ts) > 10 * eps)
       Ts = t(ii+1) - t(ii);
       ## [F,G] = c2d(a,b,Ts);
       dsys = c2d(sys, Ts);