hilb.RdGenerate Hilbert matrix of dimension n
hilb(n)Generate the Hilbert matrix H of dimension n with elements
H[i, j] = 1/(i+j-1).
(Note: This matrix is ill-conditioned, see e.g. det(hilb(6)).)
matrix of dimension n
hilb(5)
#> [,1] [,2] [,3] [,4] [,5]
#> [1,] 1.0000000 0.5000000 0.3333333 0.2500000 0.2000000
#> [2,] 0.5000000 0.3333333 0.2500000 0.2000000 0.1666667
#> [3,] 0.3333333 0.2500000 0.2000000 0.1666667 0.1428571
#> [4,] 0.2500000 0.2000000 0.1666667 0.1428571 0.1250000
#> [5,] 0.2000000 0.1666667 0.1428571 0.1250000 0.1111111