function S=test4(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;

endfunction
