function [g,h]=restr(x)
// Edgar & Himmelblau, 1988
// x0 = [-1, -2]'
// xo = [0.8229, 0.9114]'
// S(xo) = 1.3935

g =[];

h(1,1) = x(1)-2*x(2)+1;
h(2,1) = (x(1).^2)/4+x(2).^2-1 + x(3);

endfunction
