如何使用Python解决一对非线性方程式？

Question

使用Python 解决一 对非线性方程式 的（最佳）方法是什么？（Numpy，Scipy或Sympy）

例如：

• x + y ^ 2 = 4
• e ^ x + xy = 3
  
  

解决以上问题的代码片段会很棒

Answer 1

对于数值解，可以使用fsolve：

http://docs.scipy.org/doc/scipy/reference/generation/scipy.optimize.fsolve.html#scipy.optimize.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))