
function t=L_k(k,x,s)

% Oggetto: 
% Funzione che valuta nel punto "s" il "k"-simo polinomio di Lagrange 
% relativamente al set di punti x(i) con i=1,...,m.
%
% Input:
% k: indice del polinomio di Lagrange
% x: ascisse di interpolazione;
% s: punto di valutazione;
%
% Output:
% t: valore assunto dal k-simo polinomio di Lagrange, nel punto "s".

m=length(x);
inum=setdiff(1:m,k);
xnum=x(inum);

argprod=s-xnum;
num=prod(argprod);

xk=x(k);
argprod=xk-xnum;
den=prod(argprod);

t=num/den;