[解決済み] Pythonを使って2つの非線形方程式を解くには？

質問

を解決するための（最良の）方法は何でしょうか？ Pythonを使用して非一次方程式のペア。 (Numpy、ScipyまたはSympy)

例えば

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

上記のペアを解決するコードスニペットは素晴らしいでしょう。

どのように解決するのですか？

数値計算を行うには、fsolveを使用します。

http://docs.scipy.org/doc/scipy/reference/generated/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))