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

