function S=test13(x)
% Himmelblau, 1972
% x0 = [3, 1]'
% xo = [2.5, 0]'
% S(xo) = -3.25

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