La fonction solve(f,a,b,epsilon) recherche une solution de l’équation f(x)=0 entre a et b avec une précision de epsilon par dichotomie. Les images de a et b par f, f(a) et f(b), doivent être de signes différents.
from math import * def solve(f,a,b,epsilon): if f(a)*f(b)>0: return None else: while (b-a)>=epsilon: c=(a+b)/2 if f(b)*f(c)<=0: a=c else: b=c return c