# -*- coding: cp1250 -*-
from numpy import prod
wektor_xi=[1,2,4,5]
wektor_fxi=[1,4,16,25]
x=input("podaj wartosc x dla której szukamy wartosci: ")

def Interpolacja_lagrangea(vektor_miejsc,vektor_wartosci,poszukiwany_x):
    suma=0
    for i in range(len(vektor_wartosci)):
        lxi=1
        for j in range(len(vektor_miejsc)):
            if i!=j:
                #print x,wektor_xi[i],wektor_xi[i],wektor_xi[j]
                lxi*=((((poszukiwany_x-vektor_miejsc[j])*1.0)/(vektor_miejsc[i]-vektor_miejsc[j])))
        suma+=vektor_wartosci[i]*lxi
    return suma
print Interpolacja_lagrangea(wektor_xi,wektor_fxi,x)