["python","numpy","scipy","sympy"]
How can I solve a pair of nonlinear equations using Python?

What's the (best) way to solve a pair of nonlinear equations using Python? (NumPy, SciPy, or SymPy) For example: x+y^2 = 4 e^x+ xy = 3 A code snippet which solves the above pair will be great. (These are not the equations that I want to solve, but the method will be applicable to my case.)
For a numerical solution, you can use fsolve: from scipy.optimize import fsolve import math def equations(p): x, y = p return (x+y**2-4, math.exp(x) + x*y - 3) x, y = fsolve(equations, (1, 1)) print equations((x, y))