HilbertMat(N: INT): MAT over QQ
/**/ HilbertMat(3); matrix(QQ, [[1, 1/2, 1/3], [1/2, 1/3, 1/4], [1/3, 1/4, 1/5]])