function [L,g] = lucro(x)

a=129.98;
b=4000;
c=0.5;
d=25;
x0=0.02;

L=-(a-b*x(2)-c./x(1)-d*x(1)./x(2));

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

endfunction
