# Distribution of the standard normal distribution
x <- (-4000:4000)/1000
plot(x, pnorm(x), type="l", xlab="x", ylab="P(X<=x)",
	main="Standard normal distribution")