Bifurcation diagram for the logistic map

bifurcation_diagram_logistic_map.r

iters <- 1000
ns <- iters - 500
nl <- iters - ns + 1

rmin <- 2.4
rmax <- 4
rn   <- 5000

x <- rep(0,iters)

rs <- seq(rmin, rmax, length=rn)

plot(NA, NA, xlim=c(rmin, rmax), ylim=c(0,1), xlab="r", ylab="")

for (i in 1:rn) {

   r <- rs[i]

   x[1] <- 0.1

   for (j in 2:iters) {
      x[j] <- r*x[j-1]*(1-x[j-1])
   }

   points(rep(r,nl), x[ns:iters], pch=".", col=rgb(0,0,0,0.05))

}