
function demoI

% funzione 
f=@(x) log(40+0.1*x);

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

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

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

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

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

% vettore di errori relativi
relerr=abs(v-u)./abs(v);

% massimo errore relativi
maxerr=max(relerr);

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

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

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