Say I have a multivariate polynomial ring:
R.<w,x,y,z> = PolynomialRing(ZZ, 4, order='lex')
and a polynomial contains only one variable, for example:
f = w^4 - 1
How can I find the roots of f. Thx.
sage: f.univariate_polynomial().roots() [(1, 1), (-1, 1)]