% Esempio di minimi quadrati
a = -1; 
b = 1; 
n = 80;
m = 5;
xp = linspace(a,b,n);
yp = cauchy(xp);
p = polyfit(xp,yp,m-1);
X = linspace(a,b);
Y = cauchy(X);
Yp = polyval(p,X);
plot(xp,yp,'o',X,Y,'--',X,Yp,'-')
legend(['Punti di interp. n = ', num2str(n)],'funz. di Cauchy',...
    ['Polin.min.quad. m = ', num2str(m)])

