function [dg,dh]=grestr(x)
// Edgar & Himmelblau, 1988
// x0 = [-1, -2]'
// xo = [0.8229, 0.9114]'
// S(xo) = 1.3935

dg = [];
dh(1,:) = [ 1            -2         0];
dh(2,:) = [(1/2 * x(1))  (2 * x(2)) 1];

endfunction
