mu <- 50; sd <- sqrt(40)
pl <- pnorm(60, mu, sd)
pl
p2 <- 1 - pnorm(100, mu, sd)
p2
p3 <- pnorm(20, mu, sd)
p3
p4 <- qnorm(0.293, mu, sd)
p4