% matlab.pdf, p. 69
% vektorizace po radcich

function a = hilb3(n)
a=zeros(n);
for j=1:n;
a(j,1:n) = 1./((0:n-1)+j);
end;
