def myHorner (dd, xp, x):

    n= np.size(dd)
    p = dd[n-1]
    for k in range (1, n, 1):
        p = dd[n-k]+(x-xp[n-k])*p
    return p