swh:1:snp:81eadaa089e8253d8469bcef66aa332632c6c669
Tip revision: d5bfb0e3b7a3211ff2bbe2926e5ba31aa24e4df6 authored by HwB on 09 March 2012, 00:00:00 UTC
version 1.0.1
version 1.0.1
Tip revision: d5bfb0e
polyval.R
###
### POLYVAL.R Polynom
###
polyval <- function(p, x) {
if (length(x) == 0) return(c())
if (length(p) == 0) return(0 * x)
if (!is.vector(p, mode="numeric") && !is.vector(p, mode="complex"))
stop("Argument 'p' must be a real or complex vector.")
if (!is.vector(x) && !is.matrix(x))
stop("Argument 'x' must be a real or complex matrix.")
n <- length(p)
y <- outer(x[1:length(x)], (n-1):0, "^") %*% p
dim(y) <- dim(x)
return(y)
}