% Musterlösung
%
% Numerische Mathematik, SS 2013
%
% Interpolationspolynom mit dividierten Differenzen auswerten
%
% Ingolf Busch, Juli 2013

% Input:  f     Vektor der dividierten Differenzen [f_0, f_01, f_01..n]
%         x     Vektor der Stützstellen[x_0 x_1 ... x_n]
%         t     Vektor der Auswertungspunkte

% Output: yt  Vektor der Funktionswerte p_n(t) des Interpolationspolynoms

function yt = my_horner_newton(f,x,t)

yt = f(end);

for i = (length(f)-1):-1:1
    
    yt = f(i) + (t-x(i)).*yt;

end