Figure 5.18:

Fractional coverage versus adsorbate concentration for different values of the adsorption constant, K.

Figure 5.18

Code for Figure 5.18

Text of the GNU GPL.

main.py



 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
# Converted from isotherm_general.m - Langmuir adsorption isotherms
import numpy as np
from misc import save_ascii

P = np.arange(0, 251, dtype=float)
CO_s1 = (0.01 * P) / (1 + 0.01 * P)
answer = np.column_stack([P, CO_s1])

for K in range(1, 101):
    CO_s = (0.1 * K * P) / (1 + 0.1 * K * P)
    answer = np.column_stack([answer, CO_s])

save_ascii('isotherm_general.dat', answer)