def interpol (xp,x,y):
    yp=0
    for i in range (len(x)):
        g=1
        for j in range (len(x)):
            if i!=j:
                g=g*((xp-x[j])/(x[i]-x[j]))
        yp+=g*y[i]    
    return yp