residue() confusion

[Ben Abbott]

As a result of reading through Hodel's 
http://www.nabble.com/bug-in-residue.m-tf4475396.html post  I decided to
check to see how my Octave installation and my Matlab installation responded
to the example

Using Matlab v7.3
--------------------------
 num = [1 0 1];
 den = [1 0 18 0 81];
 [a,p,k] = residue(num,den)

a =

        0 - 0.0926i
   0.2222 - 0.0000i
        0 + 0.0926i
   0.2222 + 0.0000i


p =

   0.0000 + 3.0000i
   0.0000 + 3.0000i
   0.0000 - 3.0000i
   0.0000 - 3.0000i


k =

     []
--------------------------

Using Octave 2.9.13 (via Fink) on Mac OSX
--------------------------
 num = [1 0 1];
 den = [1 0 18 0 81];
 [a,p,k] = residue(num,den)

a =

  -3.0108e+06 - 1.9734e+06i
  -3.0108e+06 + 1.9734e+06i
  3.0108e+06 + 1.9734e+06i
  3.0108e+06 - 1.9734e+06i

p =

  -0.0000 + 3.0000i
  -0.0000 - 3.0000i
   0.0000 + 3.0000i
   0.0000 - 3.0000i

k = [](0x0)
e =

   1
   1
   1
   1
--------------------------

These are different from both the result that 
http://www.nabble.com/bug-in-residue.m-tf4475396.html Hodel obtained , as
well as different from 
http://www.nabble.com/bug-in-residue.m-tf4475396.html Mollet's 

Thoughts anyone?


-- 
View this message in context: 
http://www.nabble.com/residue%28%29-confusion-tf4502015.html#a12839586
Sent from the Octave - General mailing list archive at Nabble.com.