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