On 2013-09-30, Wilhelm Haager <wilhelm.haager at htlstp.ac.at> wrote:
> I think, a new bug has been introduced in Maxima 5.31.1-1:
> realroots hangs when trying to calculate the real roots
> even of a simple polynomial:
>
> (%i1) p:x^3+2*x^2-2*x-1;
> (%o1) x^3+2*x^2-2*x-1
> (%i2) realroots(p); .... HANGUP!
The problem is specific to GCL (works as expected with other Lisps).
Problem seems to be that the global variable 'rootsepsilon' is zero
while it should be nonzero. Try this: realroots(p, 1e-7); or
realroots(p), rootsepsilon=1e-7; Does it work OK?
I don't know what changed between 5.30 and 5.31 to affect it.
best
Robert Dodier