Solve crash

[parisse]

This works without user intervention with Xcas on the desktop:

Code:


s:=proot(p:=48*x*(1+x)^60 -(1+x)^60 +1);
(48*pcoeff(s)-symb2poly(p))./symb2poly(p);
point(s);

I'm afraid that multi-precision floats are required to get the roots with sufficient accuracy, compare with normal precision:

Code:


m:=companion(p);
point(diag(egvl(evalf(m))));