https://github.com/cran/pracma
Raw File
Tip revision: c79a04b5074656b36e591191eb8137b70a349932 authored by Hans W. Borchers on 30 June 2014, 00:00:00 UTC
version 1.7.0
Tip revision: c79a04b
moler.R
##
##  m o l e r . R
##


moler <- function(n) {
	if (length(n) != 1 || n != round(n))
		stop("Argument 'n' must be an integer.")
	if (n <= 0) return(c())

	A <- matrix(0, nrow = n, ncol = n)
	for (i in 1:n) {
		A[i, 1:i] <- (1:i) - 2
	}
	A <- A + t(A)
	diag(A) <- 1:n
	A
}
back to top