
function fxstar = hornerN (x,c,xstar)

% HORNERN Calcola il valore del polinomio interpolatore in x^*
% utilizzando la forma di Newton e l'algoritmo di Horner
%
% Uso:
% fxstar = hornerN (x,c,xstar)
%
% Dati di ingresso:
% x vettore dei nodi
% c vettore dei coefficienti della forma di Newton
% ordinati per indici crescenti (c_0, c_1, ... )
% xstar valore in cui si vuole valutare il polinomio
%
% Dati di uscita:
% fxstar valore di P(x^*)

% Modificato: 26 giugno 2018.

n1=length(x);
n2=length(c);
if abs(n1-n2) > 0
    beep;
    error('Inconsistente numero di componenti. Stop.');
end
u=c(n1);
for j=n1-1:-1:1
    u=u.*(xstar-x(j))+c(j);
end
fxstar=u;
