clear all;
close all;


n=5; a=-1; b=1;
x=linspace(a,b,n); %punti interpolazione equispaziati
%x=(b+a)/2-(b-a)/2*cos((2*[1:n]-1)*pi/(2*n)); %punti di chebyshev ordinati

%x=(b+a)/2-(b-a)/2*(cos((2*[1:n]-1)*pi/(2*n)))/cos(pi/(2*n)); %punti di cheb estesi
%x=-cos([0:n-1]*pi/(n-1));   %chebyshev-lobatto points


%delta=1.5/(n+1);
%delta=1.6/(n+1);
%delta=1.9/(n+1);
delta=2/(n);
xx=linspace(-1,1,5000); %punti di valutazione ATTENZIONE SIAMO fuori
                             %dall'intervallo [x(min), x(max)]
% costante di Lebesgue
[leb, L]=CostLebesgue(x,xx');

figure(1)
% plot dei polinomi di Lagrange
for k=1:size(L,1)
  plot(xx,L(k,:));
  hold on
end
plot(x,zeros(length(x),1),'ob');
title('Polinomi di Lagrange');
grid
hold off

figure(2)
% funzione di Lebesgue
for k=1:size(L,2)
  f_leb(k)=sum(abs(L(:,k)));
end
plot(xx,f_leb)
title(' funzione di Lebesgue');
grid

 figure(3)
% plot costante di Lebesgue
semilogy(1:length(leb), leb,'r-.o');
title('costante di Lebesgue')