n=10;
x=cos(linspace(0,pi,n+1))';  % n+1 punti di Chebyshev estesi in [-1,1]
V=ones(n+1,n+1);
V(:,2)=x;
for j=3:(n+1), V(:,j)=V(:,j-1).*x; end  % colonna j = x * colonna precedente
%
% V e' la matrice di Vandermonde
%
t=linspace(-1,1,1000)'; % punti di valutazione
A=ones(1000,n+1);
A(:,2)=t;
for j=3:(n+1), A(:,j)=t.*A(:,j-1); end
%
% A e' una matrice 1000 X (n+1)
% La riga i di A consiste da [1, t(i) t(i)^2 t(i)^3 .... t(i)^n]
%
B=A*inv(V);
%
% la riga i di B consiste da [L0(t(i)) L1(t(i)) L2(t(i)) ... Ln(t(i))]
% dove Lj = il j-esimo polinomio di Lagrange
%
y=sum(abs(B')); 
%
% y e' un vettore.
% il componente i di y e' la funzione di Lebesgue per x=t(i)
%
plot(t,y,x,1,'o')