C_a = 1e-2
C_b = 2.5e-2
v_a = 50.0
v = linspace (0,60,num=61,endpoint=True)
nH3Op,nHOm = ReactionTitrante (C_a, C_b, v_a, v)

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