Figure 9.8:

Transformed data set, \ln k versus 1/T.

Code for Figure 9.8

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
# Converted from onedata.m - single Arrhenius data set
import numpy as np
from misc import save_ascii

lnk0  = 1.0
E     = 100.0
ndata = 10
Tmin  = 300.0
Tmax  = 500.0
Tmeas = np.linspace(Tmin, Tmax, ndata)
X     = np.column_stack([np.ones(ndata), -1.0/Tmeas])
lnk   = X @ np.array([lnk0, E])
measvar    = 1e-3
measstddev = np.sqrt(measvar)
np.random.seed(382)
lnkmeas = lnk + measstddev * np.random.randn(ndata)

tmp = np.column_stack([1.0/Tmeas, lnkmeas])
save_ascii('onedata.dat', tmp)