R/Hilbert.R

Hilbert <- function(n)
{   ## generate the Hilbert matrix of dimension n
    n <- as.integer(n)
    i <- seq_len(n)
    new("dpoMatrix", x = c(1/outer(i - 1L, i, "+")), Dim = c(n,n))
}
bedatadriven/renjin-matrix documentation built on May 12, 2019, 10:05 a.m.