function H = hessian(x)
H = [      22  -40*x(2)
     -40*x(2)  120*x(2)^2-40*x(1)];
