# Density of the standard normal distribution x <- (-4000:4000)/1000 plot(x, dnorm(x), type="l", xlab="x", ylab="dens(X=x)", main="Standard normal density function")