## Plot the normal density
x <- seq(from = -5, to = 5, by=0.01);
y <- (2*pi)^(-1/2)*exp(-x^2/2);
plot(x,y, type = "n") +
abline(h = 0) +
lines(x,y) +
title(main = "The Standard Normal Curve")