Revision 392ae21a013fb3f518e8f9eb8efb458a55a2eca2 authored by HwB on 09 April 2011, 00:00:00 UTC, committed by Gabor Csardi on 09 April 2011, 00:00:00 UTC
1 parent 162b332
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)
mdiag(abs(-m:m)) + mdiag(r, 1) + mdiag(r, -1)
}
Computing file changes ...