Comment calculer les courbes de répartition des espèces ?

Comment calculer les courbes de répartition des espèces ?
Explication en moins de 10 minutes de l'établissement des formules pour tracer les courbes de répartition des espèces acide/bases. Le code python associé est donné...ci-après :
import numpy as np
import matplotlib.pyplot as plt


# Définir la valeur du pKa

m1 = 9.2 # valeur du pka

# Créer un tableau de valeurs de pH=x
x = np.linspace(0, 14, 1000)

# Calculer les valeurs de l'abondance en %

y1 = 100*(10**(m1))*10**(-x)/(1 + (10**(m1))*10**(-x))
y2 = 100*(10**(-m1))*10**(x)/(1 + (10**(-m1))*10**(x))

# Tracer les courbes
plt.figure(dpi=300)

plt.plot(x, y1, color='grey')
plt.plot(x, y2, color='blue')

plt.xlim(0, 14)
plt.ylim(-5, 105)

# Ajouter des éléments au graphique

plt.title('Diagramme E=f(pH) du chlore')
plt.xlabel('pH', fontsize = 16, color='black')
plt.ylabel('Abondance relative en %', fontsize = 16, color='black')
plt.axhline(0, color='black',linewidth=0.1, ls='--')
plt.axvline(0, color='black',linewidth=0.1, ls='--')

# ymin et ymax sont en coordonnées fractionnaires y
plt.axvline(x=9.2,ymin=0, ymax=1,color='black',linewidth=2, ls=':')
plt.axhline(y=50,xmin=0, xmax=1,color='black',linewidth=2, ls=':')

plt.grid(color = 'gray', linestyle = '--', linewidth = 0.2)

# Si on veut voir une légende
# plt.legend()

plt.title('Répartition des espèces acide-base pour $NH_{4}^{+}$/$NH_{3}$ dans $H_{2}O$')

plt.text(12, 90, '$NH_{3}$', fontsize = 20, color='blue')
plt.text(3, 90, '$NH_{4}^{+}$', fontsize = 20, color='grey')
plt.text(8.7, -13, '9,2', fontsize = 14, color='black')

# Pour afficher le graphique à l'ecran : plt.show() après savefig

plt.savefig('myimage.jpg')