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