# Fixpunktiteration zur Berechnung
# der Loesung von x**2=a, x>0, a>0

def g(x,a):
    y=(x**3+3*a*x)/(3*x**2+a)
    return y

a=100
#Loesung ist 10
n_iter=28
x=1e4
fehler_alt=0
for k in range(n_iter):
    x=g(x,a)
    #res=abs(x**2-a)
    if k>0:
        fehler_alt=fehler
        fehler=abs(x-10)
        if abs(fehler_alt)<1e-15:
            print("konvergiert")
            break
        quotient=fehler/fehler_alt**2    
        print(x,fehler,quotient)
    fehler=abs(x-10)