Figure 9.9:

Several replicate data sets, \ln k versus 1/T.

Code for Figure 9.9

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
# Converted from severaldata.m - 10 Arrhenius datasets
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
np.random.seed(0)
nexpts = 500
lnkmeas = lnk[:, None] + np.sqrt(measvar) * np.random.randn(ndata, nexpts)
tmp = np.column_stack([1.0/Tmeas, lnkmeas[:, :10]])
save_ascii('severaldata.dat', tmp)