N = 194
p = 0.05
print(c(N*p, N*p^2))

Rs = seq(185, 190)
lhs = ppois(Rs-185, lambda=N*p)
plot(Rs, lhs, type='l')
abline(h=0.01, col='green')
grid()