x40 <- seq(0, 40, by = 1)
y_dpois_4 <- dpois(x40, 4)
plot(x40,y_dpois_4)