# HG changeset patch
# User jwe
# Date 1200020338 0
# Node ID 204d16b1d17f62b9bc9ad39146bb081074f6e09e
# Parent  807cf4d677055e5174d9d8c80ee43ef6355d29b7
[3-0-0-branch @ 2008-01-11 02:58:58 by jwe]

diff --git a/scripts/ChangeLog b/scripts/ChangeLog
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,9 +1,13 @@
+2008-01-10  Ben Abbott  <bpabbott@mac.com>
+
+	* polynomial/mpoles.m: Avoid cases where poles could be assigned
+	to more than one multiplicity group.
+
 2008-01-09  John W. Eaton  <jwe@octave.org>
 
 	* plot/drawnow.m: Fail if filename includes a directory part that
 	does not exist.
 
-
 2008-01-07  John W. Eaton  <jwe@octave.org>
 
 	* miscellaneous/copyfile.m, miscellaneous/movefile.m:
diff --git a/scripts/polynomial/mpoles.m b/scripts/polynomial/mpoles.m
--- 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];