Unfortunately I only have maxima 5.13.0 and I can't use bfallroots. Since my polynomials are of degree 3 or 4, we have formulas for the solutions and I think I can use bfloat(solve(...)); Thanks