Bug & Patch for mpoles.m

[Ben Abbott]

There is a bug in mpoles.m which occurs when a pole is associated with  
two, or more, groups of multiplicity. For example,

octave:2> p = [.9 1 1.1];
octave:2> toler = 0.15;
octave:3> [m, indx] = mpoles (p, toler, 1)
m =

    1
    1
    1
    2

indx =

    3
    2
    2
    1

There are only three poles, but the result is indicative of four. This  
is because "1" qualifies as a multiple of "0.9" and "1.1" and  causes  
residue.m to fail (terminate with an error).

I've attached a simple patch and ChangeLog.

John, this patch to mpoles.m is independent of the patch submitted for  
residue.m and is required for the current cvs version of residue as  
well as the patched version.

Ben

-------------- next part --------------
A non-text attachment was scrubbed...
Name: mpoles.patch
Type: application/octet-stream
Size: 541 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080102/590c1649/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ChangeLog
Type: application/octet-stream
Size: 176 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080102/590c1649/attachment-0001.obj