I think the smart money on this computation is heading towards finding the
eigenvalues of the companion matrix.
So maybe what we need is a really neat bfloat eigenvalue program.
RJF
> -----Original Message-----
> From: maxima-bounces at math.utexas.edu
> [mailto:maxima-bounces at math.utexas.edu] On Behalf Of Raymond
> Toy (RT/EUS)
> Sent: Friday, August 31, 2007 2:17 PM
> To: maxima at math.utexas.edu
> Subject: Re: [Maxima] bfloat roots of polynomials
>
> Raymond Toy (RT/EUS) wrote:
> > Hope nobody minds, but I'm appending a small maxima program
> to compute
> > the roots of a polynomial with bfloat coefficients. It uses
> > Jenkins-Traub to compute the roots, much like allroots. But this
>
> Oops. I forgot the new_poly function.
>
> Ray
>
> new_poly(a) :=
> block([p : map(lambda([c], expand(rectform(abs(c)))), a)],
> p[length(p)] : -p[length(p)],
> p);
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>