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))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.