% Esempio di interpolazione con spline
a = -1; 
b = 1; 
n=10;
for k=1:4
    xp = linspace(a,b,n);  % Nodi equidistanti
    yp = cauchy(xp);
    X = linspace(a,b,4*n);
    Y = cauchy(X);
    Ys = spline(xp,yp,X);
    N(k) = n;
    err(k) = max(abs(Y-Ys));
    plot(xp,yp,'o',X,Y,X,Ys)
    title(['spline con n = ',num2str(n), ' errore = ', num2str(err(k))])
    n = n*2;
    pause
end
loglog(N,err,N,err(1)*N(1)^4./N.^4,'--r');
