Add a quadratic funtion to test the analytical Jacobian