% Smith parameters by Bommannan (partial pressure in kPa) (E taken from case 1)
%           1             2               3              4               Reaction number
k0e=[ 1.147865E+26   7.981719E-06    0.0             0.0         ]';
dEe=[-5.078400E+04   8.823000E+03    0.0             0.0         ]';
k0f=[ 1.186446E+08   4.170505E+13    4.588110E+18    4.588110E+18]';
dEf=[-1.511570E+05   2.635470E+05   -2.830580E+05   -2.830580E+05]'; % Case 1 Bommannan

kf=k0f.*exp(dEf./(R*T));                                     % Forward reaction
Kp=k0e.*exp(dEe./(R*T));                                     % Equilibrium reaction

