
function esempio5

clear all; clf; 

n=2; % GRADO INTP. TRIG. 
f=@(x) 1./(1+(x-pi).^2); % FUNZ. DI RUNGE SHIFTATA IN [0,2*pi] (PERIODICA).

fa=chebfun(f,[0,2*pi]); % CHEBFUN ALGEBRICA.
ff=chebfun(f,[0 2*pi],2*n+1, 'trig');  % SERIE DI FOUR. (INTP).

% ATTENZIONE.
% IN ATKINSON SI DICE CHE INTERPOLA NEI PUNTI j*2*pi/(2*n+1), j=0,..,2*n.
t=linspace(0,2*pi,2*n+2); % NODI EQUISPAZIATI 
AE=max(abs(f(t)-ff(t)));
fprintf('\n \t MAX. ERR. INTP. PTS.:%2.2e \n',AE);
plot(fa,'Linewidth',2); hold on; plot(ff,'Linewidth',2); 
plot(t,ff(t),'o','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',8); 
hold off;
legend('Funzione','Intp.Trig.');
