function H=matrice_hilbert(n)
% es. per n=3
%  1   1/2   1/3
% 1/2 1/3  1/4
% 1/3 1/4  1/5
%
% H(i,j)=1/(i+j-1)

H=zeros(n);
for i=1:n
    for j=1:n
        H(i,j)=1/(i+j-1);
    end
end



% secondo esercizio: grafico di cond(H_n) per n=1,...,20
% usate semilogy(x,y) invece di plot(x,y) 
