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