function [H]=Hilb(n)
H=[];
//Hilb   Hilbert matrix.
//   Hilb(n) is the n by n matrix with elements 1/(i+j-1),
//   which is a famous example of a badly conditioned matrix.
//   See InvHilb for the exact inverse.
// 
 
J = 1:n;
J = J(ones(n,1),:);
I = J';
E = ones(n,n);
H = E ./ (I+J-1);
