[Maxima] bfloat roots of polynomials

[Raymond Toy (RT/EUS)]

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);