n = 500;
j = 1:n;
A = zeros(n);
for i = 1:n
	A(i,j) = 1./(i + j - 1);
end
max(eig(A))