Figure 3.4:

Partial pressures of components A and C versus liquid-phase composition in a nonideal solution.

Figure 3.4

Code for Figure 3.4

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
19
20
21
# Converted from margules.m
import numpy as np
from misc import save_ascii

ac = 1.4
ca = 2.0
pa0 = 0.65
pc0 = 0.5

npts = 100
x = np.linspace(0, 1, npts)

pa = pa0 * x * np.exp((1 - x)**2 * (ac + 2*(ca - ac)*x))
pc = pc0 * (1 - x) * np.exp(x**2 * (ca + 2*(ac - ca)*(1 - x)))
ptot = pa + pc
paid = pa0 * x
pcid = pc0 * (1 - x)
ptotid = paid + pcid

table = np.column_stack([x, pa, pc, ptot, paid, pcid, ptotid])
save_ascii('margules.dat', table)