function S=test(x)
% lucro de um processo catalitico
% x(1) -->  catalisador (fracao molar)
% x(2) -->  vapor (razao vapor/produto)
% x0 = [0.05 0.5]
% xo = [0.0741 1.0184]
% S(xo) = -2.5926

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;

%[xo,So,nS] = complex('test16',[0.05, 400],1,[0 100],[0.3 500],[],[],[],0.01)  % sem normalização
%[xo,So,nS] = complex('test16',[0.05, 0.5],1,[0 0.1],[0.3 2],[],[],[],0.01)
