function [dg,dh]=grestr21(x)
% Papageorgiou, 1996
%
% x0 = [2 0]

% weights: w = 0.5

% xo = [0.142857 1.02041]
% S(xo) = 1.02041

F1s=0.5;
F2s=0.5;
global w
w1=w;
w2=1-w;

dg(1,:)= [x(1) -w1];
dg(2,:)= [4*(2*x(1)-1) -w2];
dg(3,:)= [-1     0];

dh=[];

% [x,S,ns,lambda]=sqp('test21','restr21',[2 0])
