clear all;
close all;
i=1;
%
% dati del problema:
%
% intervallo e numero di punti in cui valutare l'errore
%
xx=(-1:2/10000:1);
%
% numero di valutazioni dell'errore da compiere
%
nn=(2:2:200);
%
% nome della funzione da considerare
%
%nome='sin(x)';
%nome='runge(x)';
nome='abs(x).*sin(1./max(abs(x),1.e-15))';
%
figure(1)
fplot(nome,[-1 1])
%
x=xx;
yy=eval(nome);
for n=nn
    x=(-1:2/n:1);
    y=eval(nome);
    zz=interp1(x,y,xx);
    errore(i)=max(abs(yy-zz));
    i=i+1;
end
figure(2)
loglog(nn,errore)
