% esempio di interpolazione dei dati
% fenomeno di RUNGE
close all
clear all
clc


% dati lungo x
x = linspace(-5,5,9);

% dati lungo y
y = 1./(1+x.^2);

% costruzione della interpolante polinomiale
p = polyfit(x,y,length(x)-1);


% visualizzazione dei dati
xpoly = linspace(-5,5,1000); % punti su cui valutare il polinomio
ypoly = polyval(p,xpoly);    % valutazione del polinomio su un insieme di punti

% valutazione della funzione che vado ad interpolare nei nodi del polinomio
yval = 1./(1+xpoly.^2);

% disegno polinomio e i sui punti dove ho interpolato
plot(xpoly,ypoly,'b-');
hold on;
plot(xpoly,yval,'r-');
plot(x,y,'o');
hold off;
