
function A=hilbert(n)

if n <= 0
    A=[];
    return;
end

for i=1:n
    for j=1:n
        A(i,j)=1/(1+i+j);
    end
end