function res=Keq(T)

R  = 8.314;                              % Ideal gas constant  [kJ/(kmol*K)]
res=zeros(2,1);

res(1)=1.0404e6*exp(19-50784/(R*T));       % Bommannan 1989, Lee 1997 Canged 29.09.05 in order to limit reaction 1 
res(2)=9.869e-3*exp(-10.12+ 8823/(R*T));   % Bommannan 1989, Lee 1997 Canged 29.09.05 in order to increase reaction 1 
