• Couple acide‑base : ${\text{C}}_3{\text{H}}_7\text{COOH(aq)/}$${\text{C}}_3{\text{H}}_7{\text{COO}}^{-}\text{(aq),}$ ${\text{H}}_3{\text{O}}^{+}{\text{(aq)/H}}_2\text{O(l)}$ et ${\text{H}}_2{\text{O(l)/HO}}^{-}\text{(aq)}$

import numpy as np
import matplotlib.pyplot as plt

# Conditions experimentales
c_A =
c_B = 
V_A =

# Coefficients stoechiometriques
a = 
b = 
c = 
d = 

V_B = np.arange(0, 25, 0.01)
x = c_B*V_B/b
n_A = c_A*V_A - a*x
n_B = np.zeros(len(V_B))
n_C = c*x
n_D = d*x

for i in range(len(V_B)) :
	if c_B*V_B[i]/b > c_A*V_A/a :
		print(c_B*V_B[i]/b, c_A*V_A/a, i)
		x[i] = c_A*V_A/a
		n_A[i] = 0
		print(x[i])
		n_B[i] = c_B*V_B[i] - b*x[i]
		n_C[i] = c*x[i]
		n_D[i] = d*x[i]

plt.xlabel('Volume verse de solution titrante (mL)')
plt.ylabel('Quantites de matiere (mol)')
plt.title('Evolution des quantites de matiere \n en fonction du volume de solution titrante versee')
plt.grid()
plt.plot(V_B, n_A,color = 'red', label='Quantite n_A')
plt.plot(V_B, n_B,color = 'blue', label='Quantite n_B')
plt.plot(V_B, n_C,color = 'yellow', label='Quantite n_C')
plt.plot(V_B, n_D,"--", color = 'orange', label='Quantite n_D')
plt.legend()
plt.show()