Re: Maxima - precision
- From: rjf <fateman@xxxxxxxxx>
- Date: 20 May 2007 17:11:00 -0700
On May 20, 12:59 am, adam majewski <adamm...@xxxxx> wrote:
Hi
I'm computing complex roots of polynomial:
P(n):=if n=0 then 0 else P(n-1)^2+c;
bigfloat(allroots(%i*P(3)));
How can I set precision ( number of significant digits) ?
Changing fpprec seems to have no effect .
You are computing the roots of a polynomial in machine precision.
Allroots is just a transcription of the Jenkins-Traub floating point
polynomial zero-finder, and does not use bigfloats at all.
given the machine-float roots, you then convert them to bigfloat
precision. This does not make them any more accurate. At least you
would do this if you were to spell bfloat(...) correctly.
What you want is a program different from allroots that computes the
roots to a precision that depends on fpprec. This is not, so far as I
know, part of the Maxima system, though a version could easily be
produced.
It would probably have some glitches because zero-finding programs
have to deal with nasty problems.
But maybe your polynomials are not nasty.
What you can do is choose a zero computed from allroots, and then
refine it.
Usually you can use newton iteration, and run that until you are happy
with the answer.
.
- References:
- Maxima - precision
- From: adam majewski
- Maxima - precision
- Prev by Date: Jordan form of the BLOCK Frobenius matrix and its inverse
- Next by Date: Re: Help to prove a linear algebra equation
- Previous by thread: Maxima - precision
- Next by thread: Jordan form of the BLOCK Frobenius matrix and its inverse
- Index(es):
Relevant Pages
|
|
