
function esercizio1_correzione

d = [-1, 1]; % INTERVALLO.
ff = @(x) 1./(1+25*x.^2); % FZ. RUNGE
f=chebfun(ff,d);
echeb=[];
nn=1:1:100;
for n=nn
    fc = chebfun(ff, n+1); % INTERPOLAZIONE IN (n+1) NODI DI CHEBYSHEV.
    echeb=[echeb norm(f-fc,inf)];
end

mygamma=sqrt(echeb(100)/echeb(98));
% mygamma=sqrt(echeb(end)/echeb(end-2));

C=echeb(end)/(mygamma^100);

fprintf('\n \t gamma: %1.5e',mygamma);
fprintf('\n \t C: %1.5e \n \n',C);

