% esperimenti sull'interpolazione

% provate a cambiare il valore di n
n = 4;

% e a cambiare funzione, eventualmente aggiungendone altre
fun = @(x) sin(pi*x);
%fun = @(x) 1./(1+25*x.^2);

t = linspace(-1,1,n+1)';	% ascisse di interpolazione
y = fun(t);			% ordinate di interpolazione
x = linspace(-1,1,201)';	% ascisse per i grafici
p = lagint(x,t,y);		% calcolo del polinomio interpolante
f = fun(x);			% calcolo della funzione da interpolare

plot(x,f,x,p,'--',t,y,'o')
legend('f(x)','p_n(x)')

% per ampliare questo argomento si potrebbe implementare il calcolo del
% polinomio interpolante in base canonica e/o verificare che risultati si
% ottengono con i nodi di Chebychev

