q<-0.1 u<-seq(0,1,0.05) #a #1 mi1<-q/(1-u*q) plot(u,mi1,type = 'l') #2 mi2<-rep(-log(1-q),length(u)) lines(u,mi2,col='red') #3 mi3<-q/(1-(1-u)*q) lines(u,mi3,col='blue') #b #1 qx1<-u*q plot(u,qx1,type = 'l') #2 qx2<-1-(1-q)^u lines(u,qx2,col='blue') #3 qx3<-1-(1-q)/(1-(1-u)*q) lines(u,qx3,col='red')