function [g,h]=restr(x)
% lucro de um processo catalitico
% x(1) -->  catalisador (fracao molar)
% x(2) -->  vapor (razao vapor/produto)

x2=x(2)*400;

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

S=-lucro;
g(1)=S;
g(2)=-x(1);
g(3)=-x(2);
g(4)=x(1)-1;
h=[];
