%a)
T = 300
RT = 8.314*T;
deltaH = [-279 -200].*1e3;

k11 = 0.7;
k12 = k11.*exp(deltaH(2)./RT);
k21 = 0.6;
k22 = k21.*exp(deltaH(2)./RT);

%plot(T,k11,T,k12,T,k21,T,k22)
%legend('k_1_1','k_1_2','k_2_1','k_2_2')
%xlabel('Temperatur (K)')
%title('Temperaturberoendet hos hastighetskonstanterna')
%grid

%d)
T=100:600;
k11 = 20.4055*exp(-1.0111e3./T);
k12 = k11.*exp(deltaH(2)./RT);
k21 = 126.3822*exp(-1.6099e3./T);
k22 = k21.*exp(deltaH(2)./RT);
plot(T,k11,T,k12,T,k21,T,k22)
xlabel('Temperatur (K)')
title('Temperaturberoendet hos hastighetskonstanterna')
legend('k_1_1 = 20.4055exp(-1.0111e3/T)','k_1_2','k_2_1 = 126.3822exp(-1.6099e3/T)','k_2_2')
