function res = vanvusse_EE(x,rpar)
%definições das variáveis
Ca=  x(1);
Cb = x(2);

% dCadt = dxdt(1);
% dCbdt = dxdt(2);

%parâmetros externos
Ca0 = rpar(1); %Cain
T  = rpar(2) + 273.15;
F  = rpar(3);

%parâmetros comuns

%parâmetros cinéticos da eq Arrhenius
k0 = [1.287e12 1.287e12 9.043e9];
E_R = [9758.3 9758.3 8560.0];

k = k0.*exp(-E_R/T);

V  = 10;


res(1) = (F/V)*(Ca0 - Ca) - k(1)*Ca - k(3)*(Ca^2);
res(2) = -(F/V)*Cb +  k(1)*Ca - k(2)*Cb ;
