# Newt. int. pol. interpolira f-ju f(x)=x*log(x) u čvorovima 2,4,5,7

from numpy import *

x = array([2,4,5,7])
def f(x):
    return x*log(x)
a = []
n = len(x)

for i in range(n):
    a.append(x[i])

for j in range(1,n):
    for i in range(n-1,j-1,-1):
        a[i]=float(a[i]-a[i-1])/float(x[i]-x[i-j])

print (x)
print (y)
print (a)

print('p3(x)=',a[0],'+',a[1],'(x-',x[0],')+',a[2],'(x-',x[0],')(x-',x[1],')+',a[3],'(x-',x[0],')(x-',x[1],')(x-',x[2],')')