function [S,g] = opt_res(x)

S=x(1)^2+x(2)^2-3*x(1)+15*x(2);

g(1)=(x(1)+x(2))^2-4*(x(1)-x(2));

%S=S+10000*g(1)*g(1);

% x1=[1.5 0.5], S=5.5
% x2=[4 0],     S=4

