
clear all;

a=-1; b=1;
grado_vett=5:5:50;

targets=(linspace(a,b,1000))';

for index=1:length(grado_vett)
    grado=grado_vett(index);
    
    % Equispaziati: grado+1 punti.
    xeq=(linspace(a,b,grado+2))'; 
    lebeq(index,1)=lebesgue_constant(xeq,targets,grado,a,b);   
    
    % Chebyshev: grado+1 punti.
    xch=cheb(a,b,grado+1);
    lebch(index,1)=lebesgue_constant(xch,targets,grado,a,b);

    % Gauss-Legendre-Lobatto: grado+1 punti.
    xw=lobatto_jacobi(grado,a,b); 
    xll=xw(:,1);  
    lebgll(index,1)=lebesgue_constant(xll,targets,grado,a,b); 
    
end

format short e
[lebeq lebch lebgll]