function H=htest1(x)
// Hessian of Rosenbrock, 1960
// x0 = [-1.2, 1]'
// xo = [1, 1]'
// S(xo) = 0

H(1,1)=100*(12*x(1)^2-4*x(2))+2;
H(1,2)=-400*x(1);
H(2,1)=H(1,2);
H(2,2)=200;

endfunction
