function [mat] = hilb(dim)
for i = 1:dim
  for j=1:dim
   mat(i,j)=1/(i+j+1); 
end
end
//




