# Inverse distribution of the standard normal distribution p <- (1:9999)/10000 plot(p, qnorm(p), type="l", xlab="probability p", ylab="x such that p = P(X<=x)", main="Inverse of the standard normal distribution")