clear all;
close all;
%
% dati del problema:
%
% grado del polinomio
%
n=10;
%
x=linspace(-1,1,n+1);
%
% 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])
hold on;
%
p=polyfit(x,eval(nome),n);
fplot(['polyval(',mat2str(p),',x)'],[-1,1],'r')
plot(x,eval(nome),'*k')
