# HG changeset patch
# User jwe
# Date 1200020342 0
# Node ID 1e716cb4b2a3b507a561a40ab24f4592bf519ee2
# Parent  8fb8d6985395b42de1a96e9250901c7684aa9609
[project @ 2008-01-11 02:59:02 by jwe]

diff --git a/scripts/ChangeLog b/scripts/ChangeLog
--- a/scripts/ChangeLog
+++ b/scripts/ChangeLog
@@ -1,3 +1,8 @@
+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-10  John W. Eaton  <jwe@octave.org>
 
 	* plot/gnuplot_binary.m: New file.
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];