[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
I think I have found another problem in residue.
From: |
Doug Stewart |
Subject: |
I think I have found another problem in residue. |
Date: |
Sun, 09 Dec 2007 12:23:34 -0500 |
User-agent: |
Thunderbird 1.5.0.13 (Windows/20070809) |
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
- I think I have found another problem in residue.,
Doug Stewart <=