function [gr,H]=gradhess(X)
x=X(1);y=X(2);
gr=[4*(x-2)^3+2*(x-2*y); 
    -4*(x-2*y)];

H=[12*(x-2)^2+2 , -4 ;
   -4 ,         8];
