
function errore=controlla_Lk(k,x)

% Oggetto: 
% Funzione che valuta "quanto" L_k(x_j)=delta(i,j) 
%
% Input:
% k: indice del polinomio di Lagrange
% x: ascisse di interpolazione;
%
% Output:
% errore: massimo errore max(abs(L_k(x_j)-delta(i,j))).

m=length(x);

% vettore v
for i=1:m
    v(i)=L_k(k,x,x(i));
end

% vettore ek
ek=zeros(size(v));  ek(k)=1;

% errore in norma infinito tra "v" ed "ek"
errore=norm(v-ek);