Taux d'avancement final d'une réaction entre un acide et l'eau

• Un exemple de code qui fonctionne

# Importation librairie utile
import numpy as np

# Renseigner ici les données chimiques
pKa=4.19
concentration=2.5E-3

# Calcul de la constante d'acidité
Ka=10**(-pKa)

################################################
### Résolution du polynôme du deuxieme degré ###
################################################
# Coefficients du polynôme
a=1
b=Ka
c=-Ka*concentration
# Calcul du discriminant
Delta=b*b-4*a*c
# Solution positive
h=(np.sqrt(Delta)-b)/(2*a)

##############################
### Conséquences chimiques ###
##############################
# Taux d'avancement final
tau=h/concentration
# pH final
pH=-np.log10(h)

### Ecriture des resultats
print("Taux d'avancement final :",round(tau*100,2),"%")
print("pH final :",round(pH,2))

# Résultats complémentaires
print("Concentration finale de l'acide :",format(concentration-h,"3.2E"),"mol/L")
print("Concentration finale de la base :",format(h,"3.2E"),"mol/L")