python system of nonlinear equations

import scipy.optimize as opt<br /> from numpy import exp</p> <p>def f(variables) :<br /> (x,y) = variables<br /> first_eq = x + y**2 – 4<br /> second_eq = exp(x) + x*y – 3<br /> return [first_eq, second_eq]</p> <p>solution = opt.fsolve(f, (0.1, 1)) # fsolve(equations, X_0)<br /> print(solution)</p> <p># [ 0.62034452 1.83838393]