对于数值解，可以使用fsolve：
http://docs.scipy.org/doc/scipy/reference/generation/scipy.optimize.fsolve.html#scipy.optimize.fsolve
from scipy.optimize import fsolveimport mathdef 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))
