# consideriamo un numero di uova deposte da 3 fino a 100
uova <- c(3:100)

# probabilità di deporre 3,4,...,100 uova
fuova <- dpois(uova,5)

# probabilità che nascano 3 pulcini da 3,4,...,100 uova
fpulc <- dbinom(3,size=uova,prob=0.4)

# probabilità di avere esattamente 3 pulcini
sum(fuova*fpulc)