% Esempio di minimi quadrati
a = -1; 
b = 1; 
n1 = 400;
m = 24;
xp1 = linspace(a,b,n1);
yp1 = cauchy(xp1);
p1 = polyfit(xp1,yp1,m-1);

%n2 = 16;
%xp2 = linspace(a,b,n2);
%yp2 = cauchy(xp2);
%p2 = polyfit(xp2,yp2,m-1);

X = linspace(a,b);
Y = cauchy(X);
Yp1 = polyval(p1,X);
%Yp2 = polyval(p2,X);
plot(xp1,yp1,'o',X,Y,'--',X,Yp1,'-r')
legend(['Punti di interp. n1 = ', num2str(n1)],'funz. di Cauchy',...
    ['Polin.min.quad. m = ', num2str(m)])

