from z3 import *
x, y = Reals("x y")
solve(x * x + 3 * x == 3)
[x = -3.7912878474?]
solve(y ** 5 + 2 * y ** 4 + 3 * y ** 3 + y ** 2 + y == 3)
[y = 0.7121129580?]