mean <- 85
sd <- 15
cutoff <- 120
pop <- 1e+6

p <- 1-pnorm(cutoff, mean, sd)
s <- sprintf("%.2f%% of the population are above %d, making for a total of %d.", 
	         round(p*100,2), cutoff, as.integer(pop*p))
print(s)