import matplotlib.pyplot as plt

C_a = 1e-2
C_b = 2.5e-2
V_a = 50.0
V = np.linspace(0, 60, num=61, endpoint=True)
nH3Op, nHOm = ReactionTitrante(C_a, C_b, V_a, V)

plt.plot(V, nH3Op, label=r"$n_{H_3O^+}$")
plt.plot(V, nHOm, label=r"$n_{HO^-}$")
plt.legend()
plt.xlabel("Volume de solution titrante (mL)")
plt.ylabel("Quantités de matière (mol)")
plt.show()