Re: problem with consistency of root finding and sorting

[Henry F. Mollet]

on 10/10/07 10:31 AM, Ben Abbott at address@hidden wrote:

> Henry, 
> 
> Since there has been some problem with getting my version of residue to work
> on all platforms, I'm still on the look out for the unfortunate possibility of
> more bugs ... arrgghhh :-(
> 
> If by chance you are using cvs to keep everthing as bleeding-edge as possible,
> please do an update and build.
> 
> If not, please grab copies of residue.m and mpoles.m from the web interface to
> Octave's cvs server.
> 
> http://velveeta.che.wisc.edu/cgi-bin/cvsweb.cgi/octave/scripts/polynomial/
> 
> At present, I believe the versions checked in are functioning correctly.
> 
> If you place these to files into a folder (by themselves and without renaming
> them) and run octave from there, you'll be able to verify all is working by
> issuing the command "test residue"
> 
> I'd love to hear that everything works on your end, so that I can move on to
> other improvements/modifications.
> 
> Ben

Looks OK now. I've grabbed copies of residue.m and mpoles.m from the web
interface to Octave's cvs server. Downloaded Rev 1.30 of residue.m and Rev
1.2 of poles.m and saved to documents/BenTest
Henry
Details:
octave-2.9.14:1> cd documents/BenTest

octave-2.9.14:2> ls
mpoles.m        residue.m
 
octave-2.9.14:3> test residue
PASSES 2 out of 2 tests

octave-2.9.14:4> b = [1 0 1];
octave-2.9.14:5> a = [1 0 18 0 81];
octave-2.9.14:6> [r, p, k, e] = residue(b, a)
r =

   0.00000 - 0.09259i
   0.22222 + 0.00000i
   0.00000 + 0.09259i
   0.22222 - 0.00000i

p =

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

k = [](0x0)
e =

   1
   2
   1
   2

octave-2.9.14:7> [e, indx] = mpoles (p)
e =

   1
   2
   1
   2

indx =

   1
   2
   3
   4

octave-2.9.14:8>