function esempio3

fprintf('\n \t COSTANTI DI LEBESGUE, NODI CHEBYSHEV');
max_deg=50;
em=0.57721566490153286060651209008240243104215933593992; %EULERO-MASCHERONI
for deg=1:max_deg
    x=chebpts(deg+1); % PUNTI DI CHEB. IN CHEBFUN. PS: length(x)=deg+1.
    [L,Lconst]=lebesgue(x); % FUNZIONE E COSTANTE DI LEBESGUE.
    lebesgue_constants(deg)=Lconst;
    est(deg)=(2/pi)*(log(deg+1)+em+log(8/pi)); % STIMA COSTANTE LEBESGUE.
    fprintf('\n \t DEG: %3.0f LEBESGUE CONST.: %1.4e EST.:% 1.4e',...
        deg,Lconst,est(deg));
end

plot(1:max_deg,lebesgue_constants,'ro-');