logistic function double period bifurcation points (r's) on figures
change iteration range
make r range same as in lyapunov expo script