скорость сходимости метод Ньютона

Question

Используя метод Ньютона, я могу вычислить корень нелинейного уравнения. Но теперь я хотел бы рассчитать скорость сходимости с помощью этой функции:

 def rate(x, x_exact):
     e = [abs(x_ - x_exact) for x_ in x]
     print (e)
     q = [(log(e[n+1]/e[n]))/(log(e[n]/e[n-1])) for n in range(4)]
     return q

где x - список сходимости решений x. метод Ньютона постепенно добавляет x в список до тех пор, пока он не превратится в x_exact. (так что здесь последний элемент списка должен быть x_exact)

Но я не понимаю формулу, потому что если я вставлю x_exact в список, тогда x=x_exact в конце и ошибка объединится (деление на 0)

Как мне это сделать? Если я нажимаю каждый x в списке, но не последний, то скорость сходимости не будет точной.

Любые идеи? Благодарность