1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
//calculation of activation energy from rate constant and Temperature data
clear;
clc;
printf("\t Example 13.8\n");
R=8.314;//gas constant, kJ/mol
T=[700,730,760,790,810];//temperature(data given), K
x=T^-1;//1/T values corresponding to Temp values above, K-1
k=[0.011,0.035,0.105,0.343,0.789];//rate constant(data given) in 1/M^1/2 s corresponding to temperature values
lnk=log(k);//lnk values corresponding to k
[slope]=reglin(x,lnk);//lnk=-Ea/(R*T)+lnA, slope of the line between lnk and 1/T gives -Ea/R value
Ea=-slope*R;//activation energy, kJ/mol
printf("\t the activation energy for the decomposition is : %4.2f*10^2 kJ/mol\n",Ea*10^-5);
//End
|
