https://github.com/cran/pracma
Tip revision: e970b8e92676e0d1f2c84b116f0bad93b03b89c6 authored by Hans W. Borchers on 15 December 2019, 21:30:02 UTC
version 2.2.9
version 2.2.9
Tip revision: e970b8e
wilkinson.R
##
## w i l k i n s o n . R Wilkinson matrix
##
wilkinson <- function(n){
if (length(n) != 1 || n != round(n))
stop("Argument 'n' must be an integer.")
if (n <= 0) return(c())
m <- (n-1)/2.0
r <- rep(1, n-1)
Diag(abs(-m:m)) + Diag(r, 1) + Diag(r, -1)
}