newton's method

            print "Starting from",            for k in range(ITERATIONS):                print '\t', est                divisor = fprime(est)                assert divisor == fprime(est)                assert f(est) == f(est)                if divisor != 0:                    est = est - f(est)/divisor