function fx=D2Rosenbrock(x)
fx=[ 800*x(1)^2-400*(x(2)-x(1)^2)+2,  -400 *x(1);    -400*x(1), 200];
end
