function izvod1(x)
% Poziv funkcije: izvod1(2.1)

tablica;

n=length(X);
% TABLICA KONACNIH RAZLIKA

kRazlike = zeros(n,4); 
for i = 1:n-1
 kRazlike(i,1) = Y(i+1)-Y(i);
end;
for j = 2:4
    for i = 1:n-j
        kRazlike(i,j) = kRazlike(i+1,j-1) - kRazlike(i,j-1);
    end;
end;

disp([X' Y' kRazlike]);

% Izvod I Njutnovog interpolacionog polinoma
% y' = (kRazlike(1,1) + (q-0.5)*kRazlike(1,2) + 
% (3*q^2-6*q+2)*kRazlike(1,3)/6+(4*q^3-18*q^2+22*q-6)/24*krazlike(1,4))/h

h=X(2)-X(1);
Q = (x-X(1))/h;
yi = (kRazlike(1,1) + (Q-0.5)*kRazlike(1,2) + (3*Q^2-6*Q+2)*kRazlike(1,3)/6+(4*Q^3-18*Q^2+22*Q-6)/24*kRazlike(1,4))/h;

disp('Vrednost prvog izvoda funkcije f u tacki ');
x
disp('je: ');
yi
