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

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