prm <- gamma.params(mu=15.3, sigma=9.3)
prm[1]/prm[2] # mean
sqrt(prm[1]/(prm[2]^2)) # sd
plot(density(rgamma(n=1000, shape=prm[1], rate=prm[2])), main='')
abline(v=15.3, lty=2, col='red')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.