import matplotlib.pyplot as plt

import numpy as np



c_B=0.1  #concentration (mol/L) de la solution titrante 

V_A=20   #volume initial (mL) de solution titrée

V_E=20   #volume équivalent (en mL) 

c_A=c_B*V_E/V_A



V_B1=np.linspace(0,V_E,V_E*10+1)

n_a1=(c_A*V_A-c_B*V_B1)*1e-3

n_b1=c_B*V_B1*1e-3

n_B1=0*V_B1



V_B2=np.linspace(V_E,V_E*1.5,V_E*10+1)

n_a2=0*V_B2

n_b2=(0*V_B2+c_B*V_E)*1e-3

n_B2=c_B*(V_B2-V_E)*1e-3



plt.plot(V_B1,n_a1,color="blue")

plt.plot(V_B1,n_b1,color="green")



plt.plot(V_B2,n_a2,color="blue")

plt.plot(V_B2,n_b2,color="green")



plt.show()