# W&J Example page 27.

binom<-function(t,s,prob) {return (prob^s * (1-prob)^(t-s)) } #Distribuição Binomial

dat<-seq(0,10)
y<-binom(5,3,dat/dat[length(dat)])
y<-y/sum(y)

plot(dat,y,ylim=c(0,0.6),ylab="Probabilidade",xlab="Número verdadeiro de bolas vermelhas") 
lines(dat,y)

y2<-binom(50,30,dat/dat[length(dat)])
y2<-y2/sum(y2)

points(dat,y2,col="red")
lines(dat,y2,lty=("dashed"),col="red")