Figure 4.10:

Reaction rate versus concentration for nth-order kinetics, r=kc_A^n, n\leq 0, k=1 for all orders.

Code for Figure 4.10

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
# Converted from rateneg.m
import numpy as np
from misc import save_ascii

npts = 100
cfin = 3.0
cmin = 0.001
k = 1.0
orders = np.array([-2, -1, -0.5, 0], dtype=float)
norders = len(orders)
store = np.empty((npts, 0))

for n in orders:
    c = np.linspace(cmin, cfin, npts)
    rate = k * c**n
    store = np.column_stack([store, c, rate]) if store.size > 0 else np.column_stack([c, rate])

save_ascii('rateneg.dat', store)