Pr[A<3] = 0.1247.

ppois(q=2, lambda=5)

mpois <- function(r) { floor(r + 1/3 - 0.02/r) }