Mercurial > hg > octave-nkf

diff scripts/polynomial/mpoles.m @ 11612:204d16b1d17f release-3-0-x

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[3-0-0-branch @ 2008-01-11 02:58:58 by jwe]
author jwe
date Fri, 11 Jan 2008 02:58:58 +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];