# Lagrange koji apr. f-ju f(x)=(1+x**3)*cos(pi*x)
# čvorovi su vrijed. 0,0.7... itd. iz čvorovi.txt
# izračunajte apsolutnu grešku aprox. u čvoru x=1

from numpy import *
import numpy as np
from scipy import interpolate as i

def f(x):
    return (1+x**3)*cos(pi*x)

x = array([0,0.7,1.1,1.8,2.7,3.1,4])

print (x)

y=f(x)
print (y)

pl=i.lagrange(x,y)
print (pl)

print ('apsolutna greska aproksimacije u x=1: ', abs(f(1)-pl(1)))