function H=htest4(x)
// Whitte & Holst, 1964
// x0 = [-1.2, 1]'
// xo = [1, 1]'
// S(xo) = 0

//S=100*(x(2)-x(1).^3).^2+100*(1-x(1)).^2;

H(1,1)=1800*x(1)^4-1200*(x(2)-x(1).^3)*x(1)+200;
H(1,2)=-600*x(1);
H(2,1)=H(1,2);
H(2,2)=200;

endfunction
