# przykład 3.1


###  dla p=0,6

div_cond<-function(l,p)
{
ss<-0
for (m in 0:floor(l/2)) ss<-ss+choose(l,m)*p^(l-m)*(1-p)^m
ss
}

ltab<-seq(3,300,by=2) # liczba modeli skladowych
etab<-c(0.4)  # blad pojedynczego modelu
tv<-matrix(0,nrow=length(etab),ncol=length(ltab))
colnames(tv)<-ltab
rownames(tv)<-etab
for (i in etab) for (j in ltab) tv[which(etab==i),which(ltab==j)]<-div_cond(j,1-i)

plot(ltab,tv,type="l",ylim=c(0.6,1),xlab="N",ylab="P(D*)",main="p=0,6")