## This is the R implementation of the Riemann zeta function.
## NOTE: Note so well construct.
##' @export
zeta <- function(s, k = 10)
{
sVec <- as.vector(s)
nObs <- length(s)
kSeries <- 1:k
kMat <- matrix(1/kSeries, nObs, k)
out <- rowSums(kMat^sVec)
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.