% Hilbert matrices (double precision)
% Aij=1/(i+j-1)

function H=HilbertMatrix(N)

H=ones(N,N);
for k=1:N
    H(k,:)=1./(k-1+(1:N));
end
