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