global e

xo=[];
So=[];
F1=[];
F2=[];
x0=-2;
for e=0:0.1:3
   [x,S,ns,lambda]=sqp('test25e','restr25e',x0);
   x0=x;
   xo=[xo x];
   So=[So S];
   F1=[F1 atan(x+2).^2];
   F2=[F2 atan(x-2).^2];
end
e=[0:0.1:3];

plot(e,xo,'o')
figure(2);plot(xo,F1,'bo',xo,F2,'ro')
figure(3);plot(F1,F2)
