function [g,h]=restr(x)

global e
g=atan(x+2).^2-e;
h=[];

% [x,S,ns,lambda]=sqp('test25e','restr25e',2)
