[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: problem with consistency of root finding and sorting
From: |
Henry F. Mollet |
Subject: |
Re: problem with consistency of root finding and sorting |
Date: |
Wed, 10 Oct 2007 10:58:51 -0700 |
User-agent: |
Microsoft-Entourage/11.1.0.040913 |
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>
- Re: problem with consistency of root finding and sorting, (continued)
- Re: problem with consistency of root finding and sorting, David Bateman, 2007/10/10
- Re: problem with consistency of root finding and sorting, Ben Abbott, 2007/10/10
- Re: problem with consistency of root finding and sorting, David Bateman, 2007/10/10
- Re: problem with consistency of root finding and sorting, Ben Abbott, 2007/10/10
- Re: problem with consistency of root finding and sorting, Ben Abbott, 2007/10/10
Re: problem with consistency of root finding and sorting, Henry F. Mollet, 2007/10/10