https://github.com/cran/pracma
Tip revision: 3fdb68cc842f2ab3b59608f8e04d495763584f78 authored by Hans W. Borchers on 27 November 2015, 12:07:10 UTC
version 1.8.8
version 1.8.8
Tip revision: 3fdb68c
hankel.R
##
## h a n k e l . R
##
hankel <- function(a, b) {
if (!is.vector(a))
stop("Argument 'a' must be a numeric or complex vector.")
n <- length(a)
if (missing(b)) b <- c(a[n], rep(0, n-1))
if (!is.vector(b))
stop("Argument 'b' must be a numeric or complex vector.")
m <- length(b)
if (a[n] != b[1])
warning("a[n] not equal to b[1], b[1] set to a[n].")
H <- matrix(0, n, m)
for (i in 2:(n+m))
H[row(H)+col(H) == i] <- if (i <= n+1) a[i-1] else b[i-n]
return(H)
}