%a)
T = 400
RT = 8.314*T;
deltaH = [-1875 -1319 -2216.8].*1e3;

k11 = 2e-3;
k21 = 4e-4;
k31 = 6e-4;
k12 = k11.*exp(deltaH(2)./RT);
k22 = k21.*exp(deltaH(2)./RT);
k32 = k31.*exp(deltaH(3)./RT);
