R/hilb.R

Defines functions hilb

Documented in hilb

###
### $Id: hilb.R 29 2022-05-30 23:02:22Z proebuck $
###
### Hilbert matrix.
###


##-----------------------------------------------------------------------------
hilb <- function(n) {
    if (!is.numeric(n)) {
        stop(sprintf("argument %s must be numeric", sQuote("n")))
    } else if (!(length(n) == 1)) {
        stop(sprintf("argument %s must be of length 1", sQuote("n")))
    }

    n <- floor(n)
    if (n < 0) {
        return(matrix(NA, nrow=0, ncol=0))
    }

    i <- seq_len(n)
    H <- 1 / outer(i-1, i, "+")

    H
}

Try the matlab package in your browser

Any scripts or data that you put into this service are public.

matlab documentation built on July 1, 2024, 5:07 p.m.