support <- 114:125 b <- function(k) exp(lchoose(10000,k)+k*log(.01)+ (10000-k)*log(.99)) plot(support, sapply(support, function(d) 1-sum(sapply(0:d, b))), pch=1) p <- function(k) exp(k*log(100)-lfactorial(k)-100) points(support, sapply(support, function(d) 1-sum(sapply(0:d, p))), pch=2) n <- function(d) 1-pnorm(d, mean=100, sd=sqrt(99)) points(support, sapply(support, n), pch=3)