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


% dCadt = dxdt(1);
% dCbdt = dxdt(2);
% dTdt  = dxdt(3);
% dTkdt = dxdt(4);

%parâmetros externos
F  = rpar(1);
Qk = rpar(2);

%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];
DH = [-4.2 11 41.85]; %delta H

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



V  = 10;
rho = 0.9342;
Cp = 3.01;
Kw = 4032;
Ar = 0.215;
T0 = 130 + 273.15; %equivalente a Tin


Ca0 = 5.1; %equivalente a Cain

%parâmetros da camisa de resfriamento
mk = 5.0;   
cpk = 2.0; 

res(1) = (F/V)*(Ca0 - Ca) - k(1)*Ca - k(3)*(Ca^2);
res(2) = -(F/V)*Cb +  k(1)*Ca - k(2)*Cb ;
res(3) = (1/(rho*Cp))*(k(1)*Ca*DH(1) + k(2)*Cb*DH(2) + k(3)*Ca*Ca*DH(3)) +  (F/V)*(T0 - T) + ((Kw*Ar)/(rho*Cp*V))*(Tk - T);
res(4) = (1/(mk*cpk))*(Qk+Kw*Ar*(T-Tk)); %balanço na camisa


