function [g,h]=restr(x)
// Papageorgiou, 1996
//
// x0 = 2
// goals: S* = [0.5 0.5]'
// weights: w = [0.5 0.5]'
// xo = 0.1429
// S(xo) = [1.0102 1.0102]'

F1s=1;
F2s=0.5;
global w
w1=w;
w2=1-w;
g(1,:)=1+0.5*x(1).^2-w1*x(2)-F1s;
g(2,:)=4*(x(1)-0.5).^2+0.5-w2*x(2)-F2s;
g(3,:)=-x(1);

h=[];

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

endfunction
