scipy.linalg.

hilbert#

scipy.linalg.hilbert(n)[source]#

Create a Hilbert matrix of order n.

Returns the n by n array with entries h[i,j] = 1 / (i + j + 1).

Parameters:

nint: The size of the array to create.

Returns:

h(n, n) ndarray: The Hilbert matrix.

See also

invhilbert: Compute the inverse of a Hilbert matrix.

Notes

Added in version 0.10.0.

Examples

>>> from scipy.linalg import hilbert
>>> hilbert(3)
array([[ 1.        ,  0.5       ,  0.33333333],
       [ 0.5       ,  0.33333333,  0.25      ],
       [ 0.33333333,  0.25      ,  0.2       ]])