R/triDiag.R

In SoDA: Functions and Examples for "Software for Data Analysis"

Defines functions triDiag

Documented in triDiag

triDiag <- function(diagonal, upper, lower,
                    nrow = length(diagonal), ncol = nrow) {
    value <- diag(diagonal, nrow, ncol)
    R <- row(value)
    C <- col(value)
    value[C == R + 1] <- upper
    value[C == R - 1] <- lower
    value
}

Try the SoDA package in your browser