n=5;               % GRADO.
x=-5:10/n:5;        % NODI INTERP. (EQUISPAZ.).
y=cos(x);         % FUNZIONE NEI NODI EQUISP.

s=-5:10/(10*n):5;    % NODI TEST.     
t=interpol(x,y,s);  % INTERPOLANTE NEI NODI TEST.

plot(s,cos(s),'k-',s,t,'r-.',x,y,'ko'); hold on;  % PLOT INTP. VS RUNGE.
legend('cos(x)','interpolante(x)','Coppie da interpolare');
err_eqs=norm(runge(s)-t,inf);     % ABS. ERR. EQUISPAZ.
fprintf('\n \t [ABS.ERR.][EQS]: %2.2e [CHEB]: %2.2e',err_eqs);