function z=newton_eval(p,x,t)

% Auswertung eines Polynoms in Newton-Form
% p: Koeffizientenvektor p=[c0,c1,c2,...,cn] mit den
%    Koeffizienten der Newton-Polynome  1, x-x(1), ...
% x: Knotenvektor
% t: Vektor der Abszissen zur Auswertung
%
% Output: 
% z: Funktionswerte p(t)

% Hornerartige Auswertung
n=length(p)-1; % Polynomgrad
z=p(n+1)*ones(size(t));

for k=n:-1:1
    z=z.*(t-x(k))+p(k);
end