Figure 8.29:

Reaction rate versus concentration of limiting reactant; rate expression is neither convex nor concave.

Code for Figure 8.29

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
# Converted from leven_rate.m - Levenspiel rate expression
import numpy as np
from misc import save_ascii

a = 5.
b = 0.05

def rate(c, a, b):
    return c/(1. + a*c*c) + b*c

c     = np.linspace(0.01, 10., 300)
y     = rate(c, a, b)
table = np.column_stack([c, y])
save_ascii('leven_rate.dat', table)