https://github.com/cran/pracma
Raw File
Tip revision: 71455748623ef69836470c75c5f9384f6e872d45 authored by HwB on 28 June 2011, 00:00:00 UTC
version 0.6-3
Tip revision: 7145574
hilb.R
###
### HILB.R  Hilbert matrix
###


hilb <- function(n) {
    if (!is.numeric(n))
        stop("Input argument 'n' must be a numeric scalar.")
    if (n < 0) return(matrix(NA, nrow=0, ncol=0))
    if (length(n) > 1 || ceiling(n) != floor(n)) {
        n <- floor(n[1])
        warning("Size 'n' should be a single integer number.")
    }

    J <- matrix(rep(1:n, each=n), n, n)
    I <- t(J)
    E <- matrix(1, n, n)
    H <- E / (I + J - 1)
    return(H)
}
back to top