function koef=Njutn1()
%Funkcija koef=Njutn1() vraca koeficijente koef I Njutnovog 
%interpolacionog polinoma po promenljivoj q konstruisanog koriscenjem
%svih vrednosti iz tablica.m

tablica;
n=length(X);

KR=zeros(n,n-1);%tablica konacnih razlika
for i=1:n-1
    KR(i,1)=F(i+1)-F(i);
end


for j=2:n-1
    for i=1:n-j
        KR(i,j)=KR(i+1,j-1)-KR(i,j-1);
    end
    
end
%disp([X' Y' KR]);

% I Njutnov interpolacioni polinom po promenljivoj q:
yk=F(1);
qk=[1 0];
for j=1:n-1
    yk=[0 yk]+qk*KR(1,j)/factorial(j);
    qk=conv(qk,[1 -j]);
end
koef=yk;