Figure 9.7:

Measured rate constant at several temperatures.

Figure 9.7

Code for Figure 9.7

Text of the GNU GPL.

main.py



 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# Converted from untransform.m - Arrhenius data in original (untransformed) scale
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(382)
lnkmeas = lnk + np.sqrt(measvar) * np.random.randn(ndata)
data    = np.column_stack([Tmeas, np.exp(lnkmeas)])
save_ascii('untransform.dat', data)