drootinvgamma <-
function(x,shape,rate,log=F) {
output = log(2)+shape*log(rate)-lgamma(shape)-(2*shape+1)*log(x)-rate/x^2
if (log==F) output = exp(output)
return(output)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.