Mercurial > hg > octave-lyh

diff scripts/polynomial/mpoles.m @ 7360:1e716cb4b2a3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2008-01-11 02:59:02 by jwe]
author jwe
date Fri, 11 Jan 2008 02:59:02 +0000
parents aeeb646f6538
children eb63fbe60fab
line wrap: on
line diff
--- a/scripts/polynomial/mpoles.m
+++ b/scripts/polynomial/mpoles.m
@@ -100,6 +100,10 @@
       p0 = abs (p(n));
     endif
     k = find (dp < tol * p0);
+    ## Poles can only be members of one multiplicity group.
+    if (numel (indx))
+      k = k(! ismember (k, indx));
+    endif
     m = 1:numel (k);
     multp(k) = m;
     indx = [indx; k];