% Input: a, vettore dei coefficienti a1,a2,a3,...,aN
%        x, vettore dei valori di x 
% Output, p, polinomio valutato nei punti x
%         p(x) = a1 + a2x + a3x^2 + ... + aNx^(N-1)
function p = calcola_polinomio(a, x)
  p = 0*x ; % inizializzazione polinomio a 0
  N = length(a) ;
  for i = N-1 : -1 : 0 
    p = a(i+1) + x .* p ; % gli indici dei vettori partono sempre da 1 
  end
end