Andrej Vodopivec wrote:
> The issue is how to force the complex algorithm from algsys (for the
> initial system allroots gets called with the above polynomial).
> Multiplying with %i does not work in this case. Maybe maxima should
> try with the complex algorithm if the special algorithm for real
> polynomials fails.
>
Out of curiosity, I tried the original Fortran code (TOMS algorithm 493)
on this polynomial. It computes the roots just fine, so I conclude that
the translation in maxima is wrong in some way.
I'll add this to the bug list. :-( But it shouldn't be too hard to fix
this by comparing our translation with the Fortran code.
Ray