function jac = K(x)
% x = vector of x's
% jac = Jacobian matrix of partial derivatives
% Used for Example 6.2

jac = [2*x(1)+x(2), x(1); 3*x(2)^2, 1+6*x(1)*x(2)]; 

