# Berechnung von log(2) mit dem Newton-Verfahren

import numpy as np

def f(x):
    y=np.exp(x)-2
    return y

def df(x):
    y=np.exp(x)
    return y

x=7
n_iter=15

xlim=np.log(2)
print("Vergleichswert log(2)=",np.log(2))

fehler_alt=0


for k in range(n_iter):
    x=x - f(x)/df(x)
#    print(x)
    if k>0:
        fehler_alt=fehler
        fehler=abs(x-xlim)
        if abs(fehler_alt)<1e-15:
            print("konvergiert")
            break
        quotient=fehler/fehler_alt**2
        print(x,fehler,quotient)
    fehler=abs(x-xlim)