python system of nonlinear equations

import scipy.optimize as opt
from numpy import exp

def f(variables) :
    (x,y) = variables
    first_eq = x + y**2 - 4
    second_eq = exp(x) + x*y - 3
    return [first_eq, second_eq]

solution = opt.fsolve(f, (0.1, 1)) # fsolve(equations, X_0)
print(solution)

# [ 0.62034452  1.83838393]