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