#Plots fŸr die Gleichung g(z)=z
z<-seq(0,1,by=0.01)
par(mfrow=c(2,1))

#mu ˛ 1
g<-exp(-0.8*(1-z))
plot(z,g,type="l",ylim=c(0,1),xlim=c(0,1),ylab="g(z)",main="(a) g'(1) < 1 und q=1")
lines(z,z)
abline(h=1,lty=2)
abline(v=1,lty=2)

#mu > 1
g<-exp(-2*(1-z))
plot(z,g,type="l",ylim=c(0,1),xlim=c(0,1),ylab="g(z)",main="(b) g'(1) > 1 und 0 < q < 1")
lines(z,z)
abline(h=1,lty=2)
abline(v=1,lty=2)

#Plots fŸr Aussterbewahrscheinlichkeiten
#Geometrische
xg<-seq(1,5,0.01)
pg<-1/xg

#Poisson
pp<-seq(0,1,0.01)
xp<--log(pp)/(1-pp)

#Plot
plot(xg,pg,xlim=c(0,5),ylim=c(0,1),main="Aussterbewahrscheinlichkeiten",sub="__ Poisson, -- Geometrische",type="l",lty=2,xlab="E[Kinderzahl]",ylab="")
lines(xp,pp)

#Horizontale Linie hinzufŸgen
A<-c(0,1)
B<-c(1,1)
lines(A,B)