% Musterlösung
%
% Numerische Mathematik, SS 2013
%
% Interpolationspolynom auswertn mit Neville-Aitken
%
% Ingolf Busch, Juli 2013

% Input:  x Vektor der Stützstellen [x_0 x_1 ... x_n]
%         y Vektor der Funktionswerten [y_0 y_1 ... y_n]
%         t Vektor der Auswertungspunkte


% Output: ft  Vektor der Funtionswerte p_n(t) 

function yt = my_neville_aitken(x,y,t)

n = length(x);
p=cell(size(y));
for i=1:length(y)
    p{i} = y(i)*ones(size(t));
end

for k = (n-1):-1:1
    for j=1:k
        p{j} = (( t - x(j) ).*p{j+1} - ( t - x(n-k+j) ).*p{j}) ./ ( x(n-k+j) - x(j) );
    end
end

yt = p{1};

