
function demoI

% funzione 
f=@(x) log(20+x.*sin(x));

% punti di Chebyshev estesi
x=chebyshev_estesi(100,-10,10);

% punti test
t=linspace(-10,10,1000);

% valutazione della funzione nei punti "x"
y=feval(f,x);

% valutazione interpolante spline nei punti "t"
u=interp1(x,y,t);

% valutazione della funzione nei punti "t"
v=feval(f,t);

% vettore di errori assoluti
abserr=abs(v-u);

% massimo errore assoluto
maxerr=max(abserr);

% stampa massimo errore
fprintf('\n \t maxerr: %1.4e \n \n',maxerr);

% salvataggio su file
fid=fopen('errori_assoluti.txt','w');
for k=1:1000
    fprintf(fid,'\n \t %1.15e | %1.15e',t(k),abserr(k));
end
fclose(fid);

% grafico
semilogy(t,abserr); 
hold on;
title('Errori assoluti');
hold off
print('grafico','-dpdf')