## Monoacide de pKa = 4,8
import numpy as np
import matplotlib.pyplot as plt

pH = np.linspace(0,14,100)
h = 10 ** (-pH)

pKa = 4.8
Ka = 10 ** (-pKa)

alpha0 = 1 / (1 + Ka / h)
alpha1 = 1 - alpha0

plt.figure()
plt.plot(pH, alpha0, label=r"$\alpha_0$")
plt.plot(pH, alpha1, label=r"$\alpha_1$")
plt.legend()
plt.grid()
plt.show()

##
# Tétraacide
import numpy as np
import matplotlib.pyplot as plt

pH = np.linspace(0,14,100)
h = 10 ** (-pH)

pKa1 = 1.5
Ka1 = 10 ** (-pKa1)

pKa2 = 2.3
Ka2 = 10 ** (-pKa2)

pKa3 = 6.6
Ka3 = 10 ** (-pKa3)

pKa4 = 9.3
Ka4 = 10 ** (-pKa4)

alpha0 = 1 / (1 + Ka1 / h + Ka1 * Ka2 / h ** 2 + Ka1 * Ka2 * Ka3 / h ** 3 + Ka1 * Ka2 * Ka3 * Ka4 / h ** 4)

alpha1 = Ka1 * alpha0 / h
alpha2 = Ka2 * alpha1 / h
alpha3 = Ka3 * alpha2 / h
alpha4 = Ka4 * alpha3 / h


plt.figure()
plt.plot(pH, alpha0, label=r"5")
plt.plot(pH, alpha1, '.-', label=r"4")
plt.plot(pH, alpha2, '--', label=r"3")
plt.plot(pH, alpha3, '.', label=r"2")
plt.plot(pH, alpha4, '-.', label=r"1")

plt.legend()
plt.grid()
plt.show()


## Triacide Acide citrique
# Tétraacide
import numpy as np
import matplotlib.pyplot as plt

pH = np.linspace(0,14,100)
h = 10 ** (-pH)

pKa1 = 3.14
Ka1 = 10 ** (-pKa1)

pKa2 = 4.76
Ka2 = 10 ** (-pKa2)

pKa3 = 6,40
Ka3 = 10 ** (-pKa3)

alpha0 = 1 / (1 + Ka1 / h + Ka1 * Ka2 / h ** 2 + Ka1 * Ka2 * Ka3 / h ** 3)

alpha1 = Ka1 * alpha0 / h
alpha2 = Ka2 * alpha1 / h
alpha3 = Ka3 * alpha2 / h

plt.figure()
plt.plot(pH, alpha0, label=r"4")
plt.plot(pH, alpha1, '.-', label=r"3")
plt.plot(pH, alpha2, '--', label=r"2")
plt.plot(pH, alpha3, '.', label=r"1")

plt.legend()
plt.grid()
plt.show()