---

def f(x):
    return np.array([
        1/(x[0] + 1.5) - x[1] - 1,
        (1.2 * x[0])**3 - x[1] - 1
    ])