[Python] 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.)

Continue reading...