Methode de newton pour resolution d'equation

Bonjour ,

je fais appel à vous car j'ai un petit problème et je ne vois vraiment pas d'où vient le problème.

Je cherche à calculer les solutions d'une d'une fonction de classe C2 en python grace a la methode de newton en utilisant la formule

f'(x)=f(x0+h)-f(x0-h))/2*h quand h tends vers 0 mais ça bug et j'arrive pas a detecter le probleme .

Merci d'avance pour votre aide

voila ce que j'ai essayé

def appro_derivee(f,h,x0):


return (f(x0+h)-f(x0-h))/2*h

def f(x) :

return x**2-2

def resolution_Newton(f,u0,epsilon,h):

u=u0

v=u-f(u)/appro_derivee(f,h,u)

while abs(v-u)>epsilon :

u,v =v,v-f(v)/appro_derivee(f,h,v)

return v