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, 107.1), labels=c("60", "107,1"), tick=FALSE)
par(new=TRUE)
xx <-c(60.1:107.1, 107.1:60.1)
yy <-c(c(dnorm(c(60.1:107.1), mean=100, sd=10)),c(rep(0,48)))
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=107.1, col="black", lwd=2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.