function y = horner(a,x)
% Input:  Vektor a mit Polynomkoeffizienten.
%         Vektor x der Laenge k.
% Output: Vektor y mit Werten des Polynoms
%         an den Stellen x(1), ..., x(k).
n = length(a)-1; k = length(x);
y(1:k,1) = a(n+1);
for j = n:-1:1,
   y = a(j) + y.*x;
end