rMVN <- function(n, mean, sigma){
L <- mroot(sigma)
lL <- ncol(L)
t(mean + L%*%matrix(rnorm(lL*n), lL, n))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.