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];
