def Distributions (pH, pKa): E=10**(pKa-pH) a = E/(E+1) b = 1/(E+1) return a,b pH = linspace (0,14) acide, base = Distributions (pH, 5.0) plot(pH,acide,'b') plot(pH,base ,'b--')