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