function [v, g, H] = f(x) 

v = 2*x(1)^4 + 3*x(2)^4 ...
    + 2*x(1)^2 + 4*x(2)^2 + x(1)*x(2) ...
    - 3*x(1) - 2*x(2)  ;

g = [ 8*x(1)^3 + 4*x(1) + x(2) - 3
      12*x(2)^3 + 8*x(2) + x(1) - 2];

H = [ 24*x(1)^2+4   1
          1        36*x(2)^2+8];
end

