Simultit

from scipy.optimize import fsolve #fonction pour résoudre une équation
from math import log10
equation = 'x - (10 ** -14) / x' # électroneutralité
def f(x): #fonction à résoudre pour trouver la concentration en ion hydrogène
return eval(equation)
c = input('Concentration acide fort ?') # on ajoute un acide fort de concentration c
equation = equation + ' - ' + c #equation après ajout d'un acide fort
print ('equation à résoudre : ' + equation)
print(round(- log10(fsolve(f,0.5)),3))
#calcul du pH à partir de la concentration en ion hydogène
#3 chiffres après la virgule