function [dg,dh]=restr(x)
// x0 = [1, 1]
// xo = [0, 0]
// S(xo) = 0

dg=[];
dh=[1;4*x(2).^3];

endfunction
