allocatemem(20*10^6);
do(P)=nfrootsof1(nfinit(P))[1];

do(polcyclo(23))
do(polresultant(y^3*x^3-y^2*x^2+1, polcyclo(23)))
do(x^54+9*x^51+18*x^48-81*x^45+387*x^42-729*x^39+1953*x^36-7560*x^33+14229*x^30-12393*x^27-270*x^24+6156*x^21+26136*x^18-77679*x^15+88452*x^12-49572*x^9+10287*x^6+972*x^3+27)
do(x^2+396735)
do(x^2+4372152)
do(x^2+x+99184)
do(x^16+2*x^15-x^14-4*x^13+x^12+4*x^11-2*x^9-3*x^8+7*x^6-9*x^4+4*x^3+4*x^2-4*x+1)
do(polcyclo(68))
do(polcyclo(85))
nfrootsof1(x^2023+x+1)
nfrootsof1(x^2024+x-1)

\\ errors
nfrootsof1(x^0)
