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;

endfunction
