Mercurial > hg > octave-nkf

diff scripts/polynomial/residue.m @ 5463:ab7973a599cf

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2005-09-22 21:25:59 by jwe]
author jwe
date Thu, 22 Sep 2005 21:25:59 +0000
parents 4c8a2e4e0717
children 9d0aab486882
line wrap: on
line diff
--- a/scripts/polynomial/residue.m
+++ b/scripts/polynomial/residue.m
@@ -205,9 +205,16 @@
   p = roots (a);
   lp = length (p);
 
+
+  ## Determine if the poles are (effectively) zero.
+  index = find (abs (p) < toler);
+  if (length (index) != 0)
+    p (index) = 0;
+  endif
+
   ## Determine if the poles are (effectively) real.
 
-  index = find (abs (imag (p) ./ real (p)) < toler);
+  index = find (abs(p)>=toler && ( abs(imag(p)) ./ abs(p) < toler ));
   if (length (index) != 0)
     p (index) = real (p (index));
   endif