function  [jacob,ires] = vanvusse_jac_2eq(t,x,dxdt,cj,par)
%definições das variáveis
Ca=  x(1);
%Cb = x(2);
% T  = x(3);
% Tk = x(4);
%parâmetros externos
FV  = par(1);
% Qk = par(2);

%parâmetros comuns


T = 135 + 273.15;
%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);



A(1,1) = - FV - k(1) - 2*Ca*k(3);
A(1,2) = 0;
A(2,1) = k(1);
A(2,2) = - FV - k(2);


jacob = A - cj*eye(2,2);

jacob = jacob';

ires = 0;
 

