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

endfunction
