Newton's method for root of 33

newton.py

def newton_method(x0, fn, dfn):
    xn = x0
    while True:
        xn -= fn(xn) / dfn(xn)
        yield xn

root_33 = newton_method(1, lambda x: 33 - x*x, lambda x: -2*x)

for x in range(6):
    print(next(root_33))