newton rhapson

def f(x):
    return x**x-10.
def NR(x):
    e=10**-10#make this smaller to be more accurate
    return x-e*f(x)/(f(x+e)-f(x))
s=2
for i in range(10):
    print str(i)+"th interation, x="+str(s)
    s=NR(s)