def zpn(f,fprime,a):
    for k in range(10):
        a = a-f(a)/fprime(a)
    return a