// 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]'

global e

xo=[];
So=[];
F1=[];
F2=[];
for e=1:0.01:1.125
   [x,S,ns,lambda]=sqp('test20','restr20',2);
   xo=[xo x];
   So=[So S];
   F1=[F1 1+0.5*x.^2];
   F2=[F2 4*(x-0.5).^2+0.5];
end
e=[1:0.01:1.125];
