my_rdirichlet <- function (n, alpha)
{
l <- length(alpha)
x <- matrix(rgamma(l * n, alpha), ncol = l, byrow = TRUE)
sm <- x %*% rep(1, l)
x/as.vector(sm)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.