answered 1 year ago
If you need exact results, you can use:
x = PolynomialRing(QQbar, 'x').gen() f = x^4 + 1 rs = f.roots(ring=QQbar,multiplicities=False) [f(u).radical_expression() for u in rs] [0, 0, 0, 0]