La question : Compléter la ligne en gras [pH.append()] qui permet de calculer le pH de la solution après équivalence.
Le code :
<code>from math import log10
Bibliothèques importées
import numpy as np
import matplotlib.pyplot as plt
Définition des variables
CAHi : concentration initiale en acide éthanoïque; CB concentration 'de la solution contenant le réactif titrant VA : volume de solution d'acide dosé VBE : volume équivalent VEAU : Volume d eau
CAi,CB,VA,VBE,VEAU = 0.013, 0.025, 20, 10.5, 20
Initialisation de 6 listes (concentrations en AH, A-, HO-, H3O+ pH et volume
CAH, CA, CHO, CH, pH, V =[], [], [], [], [], []
Remplissage des 6 listes suivant le volume VB versé
for VB in np.arange(1,26,0.5) :
-----if VB < VBE :
---------- V.append(VB)
---------- CAH.append((CAi*VA-CB*VB)/(VA+VB+VEAU))
---------- CA.append((CB*VB)/(VA+VB+VEAU))
---------- CHO.append(0)
La valeur 999999 indiquée pour la concentration en ion oxonium est simplement là pour indiquer que ce programme ne cherche pas à déterminer la concentration de cet ion avant l équivalence
----------CH.append(999999)
----------pH.append(4.8 + log10((CB*VB)/(CAi*VA-CB*VB)))
-----elif VB > VBE:
----------V.append(VB)
----------CAH.append(0)
----------CA.append((CB*VBE)/(VA+VB+VEAU))
----------CHO.append((CB*(VB-VBE))/(VA+VB+VEAU))
----------CH.append(10**-14/((CB*(VB-VBE))/(VA+VB+VEAU)))
----------pH.append()
Affichage des graphiques
plt.figure(figsize=(16,10))
plt.subplot(121)
plt.plot(V, CAH, 'bx-', label="CAH")
plt.plot(V, CA, 'rx-', label="CA-")
plt.legend()
plt.subplot(122)
plt.ylim(0,14)
plt.plot(V, pH, 'rx-', label="pH")
plt.legend()
plt.show()</code>
