//w=0.5;
global w

xo=[];
So=[];
F1=[];
F2=[];
for w=0:0.01:1
   [x,S,ns]=bfgs('test25',2*(2*(w<0.5)-1));
   xo=[xo x];
   So=[So S];

   F1=[F1 atan(x+2).^2];
   F2=[F2 atan(x-2).^2];
end
w=[0:0.01:1];
