import sympy as sy
x, y = sy.symbols("x y")
# nsolve needs the (in this case: two) equations, the names of the variables
# (x,y) we try to evaluate solutions for, and an initial guess (1,1) for the
# solution
print sy.nsolve((x**3+sy.exp(y)-4,x+3*y),(x,y),(1,1))
The result s...