%Spline 1 - esempio di interpolazione con le spline
% Altra funzione
figure(1)
clf
n = 8;
x = linspace(-1,1,n+1);
y = 1./(1+25*x.^2);
X = linspace(-1,1,1000);
Y = 1./(1+25*X.^2);
plot(x,y,'o',X,Y)
hold on
pp = spline(x,y);
Ys = ppval(pp,X);
plot(X,Ys,'--r')
legend('Nodi','Cauchy','Spline')
axis([-1 1 -0.5 1])
disp(['massimo errore  con ',num2str(n),' intervalli = ',num2str(max(abs(Y-Ys)))])