function S = test(x)

%w=0.5;
global w
S=w*atan(x+2).^2+(1-w)*atan(x-2).^2;

% [x,S,ns]=bfgs('test25',2)
