# NOT RUN {
hilbert <- function(n) { i <- seq_len(n); 1/outer(i - 1L, i, "+") }
H12 <- hilbert(12)
rankMM(H12) # -> 11 - numerically more realistic
rankMM(H12, tol=0) # -> 12
B <- orth(H12)
t(B) %*% B
## pracma::subspace(H12, B)
# }
Run the code above in your browser using DataLab