import numpy as np # Numpy permet de faire des calculs simples sur des arrays
import matplotlib.pyplot as plt # Matplotlib pour tracer des graphes

# Solution commerciale utilisée
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
t = 0.23    #titre massique de la solution commerciale
d = 1.18    #densité de la solution commerciale
r = 5       #rapport de dilution 
N = 18      #nombre de mesures réalisées
c0 = 1000 * d * t /36.46    #concentration en mol/L de la solution commerciale

# Concentrations des solutions 
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
c = []
for i in range (18):
    c.append(c0 / r**i)     # calcul des concentrations des solutions diluées

# Résultats des mesures de pH
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~
pHmes = []

# Valeurs calculées du pH
# ~~~~~~~~~~~~~~~~~~~~~~~
# pHcalc1 =         # valeurs calculées avec pH=-log(c)
# pHcalc2 =         # valeurs calculées en tenant compte de l'autoprotolyse

# Tracé du graphique
# ~~~~~~~~~~~~~~~~~~
fig=plt.figure(figsize=(16,9))
fig.patch.set_facecolor([0.8,0.8,0.8])

plt.plot(-np.log10(c), pHmes, 'o', color = 'navy', label = 'valeurs mesurées') # représentation des valeurs mesurées
#plt.plot(-np.log10(c), pHcalc1, '-', color = 'steelblue', label = 'valeurs calculées avec pH = -log(c)') # représentation des valeurs calculées
#plt.plot(-np.log10(c), pHcalc2, '-', color = 'plum', label = 'valeurs calculées en tenant compte de l\'autoprotolyse') # représentation des valeurs calculées

plt.xlabel('-log(c)')
plt.ylabel('pH')
plt.grid()
plt.legend()

plt.show() # affichage de la figure