function [L,g] = catalis(x)

L=(2+10*x(1)+20*x(1)^2)+(1+0.003*x(2)+2e-6*x(2)^2)-...
  50*(0.1+0.3*x(1)+1e-4*x(2)-1e-4*x(1)*x(2));

g(1)=L;
g(2)=-x(1);
g(3)=-x(2);
g(4)=x(1)-1;

endfunction
