function [y] = horner(a, z)
% Horner algoritmus
% a függvény kiszámolja a Horner-elrendezéssel az
%
% y = a(1)*z^n + a(2)*z^(n-1) + ... + a(n)*z + a(n+1)
%
% helyettesítési értéket és ezt adja vissza.
% 
n  = length(a) -1;
b  = zeros(n+1,1); % ebben számítjuk a Ruffini-sorozat tagjait
b(1) = a(1);
for j = 2 : n +1
    b(j) = a(j) + b(j-1)*z ;
end;
y = b(n +1);
end