curve(from=60, to=140, dnorm(x, mean=100, sd=10), xaxt="n", ylab="f(x)", xlab="x", col="blue", ylim=c(0.0,0.04), lty=1, lwd=4, font.lab=2, "xaxs"="i" ,"yaxs"="i", bty="l")
axis( side=1, at=c(60, 125), tick=FALSE)
par(new=TRUE)
xx <-c(60:125, 125:60)
yy <-c(c(dnorm(c(60:125), mean=100, sd=10)),c(rep(0,66)))
polygon(xx, yy, col="grey", border=NA)
par(new=TRUE)
curve(from=60, to=140, dnorm(x, mean=100, sd=10), xaxt="n", ylab="f(x)", xlab="x", col="blue", ylim=c(0.0,0.04), lty=1, lwd=4, font.lab=2, "xaxs"="i" ,"yaxs"="i", bty="l")
abline(v=125, col="black", lwd=2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.