dnormR <- function(x, mu=0, s=1, log. = FALSE) {
out <- -(log(2*pi*s*s) + ((x - mu)/s)^2)/2
ifelse(log., out, exp(out))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.