ctr=0
def computeRoot(poly, x_0, epsilon):
    ctr=0
    y=42
    while not(abs(y)<=abs(epsilon)):
        ctr=ctr+1
        y=0
        dy=0
        for thing in range(len(poly)):
            y+=poly[thing]*(x_0**thing)
            dy+=(thing-1)*poly[thing]*(x_0**(thing-1))
        x_0=x_0-(y/dy)
    return [x_0,ctr]