solve(x ** 3 + 2 * x + 3, x)
p = Symbol('p')
q = Symbol('q')
sorted(solve(x ** 2 + p * x + q, x))
solve_poly_system([y - x, x - 5], x, y)
solve_poly_system([y ** 2 - x ** 3 + 1, y * x], x, y)
Solving polynomial equations in Python