The square Hilbert matrix with elements H(i; j)=(i+j−1)−1 is widely used for testing inversion
algorithms and linear equation solvers. Not only is it symmetric positive definite, but it is also totally
positive. It is very ill-conditioned for even moderate values of order n (see Table 1).