[c,ceq]=alkylation_constraints(x)

% Title Solution of an Alkylation Process via Constrained Nonlinear Programming

% Reference: Edgar, Himmelblau and Lasdon in Optimization of
% Chemical Processes, McGraw-Hill International Edition, 2001

% Author of the matlab file: Housam Binous, PhD, Tunis, Tunisia
% Email: [e-mail]

% The same problem was solved using Mathematica and GAMS.
% GAMS ans Mathematica files are available from the author upon request.

function [c,ceq]=alkylation_constraints(x)

ceq(1)=x(5)-1.22*x(4)+x(1);
ceq(2)=x(6)-98000*x(3)/(x(4)*x(9)+1000*x(3));
ceq(3)=x(8)-(x(2)+x(5))/x(1);

c(1)=-(x(1)*(1.12+0.13167*x(8)-0.00667*(x(8)^2)))+99/100*x(4);
c(2)=(x(1)*(1.12+0.13167*x(8)-0.00667*(x(8)^2)))-100/99*x(4);

c(3)=-(86.35+1.098*x(8)-0.038*(x(8)^2)+0.325*(x(6)-89))+99/100*x(7);
c(4)=(86.35+1.098*x(8)-0.038*(x(8)^2)+0.325*(x(6)-89))-100/99*x(7);

c(5)=-(35.82-0.222*x(10))+9/10*x(9);
c(6)=(35.82-0.222*x(10))-10/9*x(9);

c(7)=-(-133+3*x(7))+99/100*x(10);
c(8)=(-133+3*x(7))-100/99*x(10);