bfloat roots of polynomials



On 9/4/07, Barton Willis <willisb at unk.edu> wrote:
>   (%i4) list_of_poly_coeffs(p,x) := (p : rat(p,x),
> makelist(ratcoef(p,x,k),k,0,hipow(p,x)))$

Alas, this won't work for bfloats.  For floats, you can set
keepfloat:true.  If I remember correctly, implementing keepfloat:all
might be non-trivial because the rat package assumes that coefficients
are atoms.  Or maybe it's enough to change its test (don't have source
here at work).

         -s