I think I have found another problem in residue.

[Doug Stewart]

I would like Ben Abbott to try the following code

n=poly([ -11 -10 ])*1000

d=poly([ -3-3j -3+3j -3-3j -3+3j -3-3j -3+3j -5.6 -1 ])

roots(d)

residue(n,d)


I get

p =

-5.60000 + 0.00000i

-3.00010 + 3.00003i

-3.00010 - 3.00003i

-2.99992 + 3.00007i

-2.99992 - 3.00007i

-2.99998 + 2.99990i

-2.99998 - 2.99990i

-1.00000 + 0.00000i

ans =

-1.3195 - 0.0000i

-3.7935 + 8.3944i

14.5612 + 10.6427i

16.3548 - 13.2580i

-3.7935 - 8.3944i

14.5612 - 10.6427i

16.3548 + 13.2580i

8.9065 - 0.0000i


but it should be:

p =

-5.60000 + 0.00000i

-3.00010 + 3.00003i

-3.00010 - 3.00003i

-2.99992 + 3.00007i

-2.99992 - 3.00007i

-2.99998 + 2.99990i

-2.99998 - 2.99990i

-1.00000 + 0.00000i

ans =

-1.3195 - 0.0000i

-3.7935 + 8.3944i
-3.7935 - 8.3944i

14.5612 + 10.6427i

14.5612 - 10.6427i

16.3548 - 13.2580i
16.3548 + 13.2580i

8.9065 - 0.0000i


Could you verify what I am getting please.?
Doug Stewart