start = np.array([0.0, 0.0])
result = minimize(
    objective,
    start,
    method="BFGS",
    jac=gradient,
    options={"gtol": 1e-10},
)