function S = test25(x)

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

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

endfunction
