# For fractional powers...
mat.frac.pow <- function(x, p){
d <- eigen(x)
return(d$vectors %*% diag(d$values^p) %*% solve(d$vectors))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.