%a)
T = 400
RT = 8.314*T;
deltaH = [-131.5 -124.4].*1e3;

k11 = 0.412;
k12 = k11.*exp(deltaH(2)./RT);
k21 = 0.055;
k22 = k21.*exp(deltaH(2)./RT);

