Find an orbit of period 5 for the chaotic logistic function $g(x)=4x(1-x)$.
Note: There's a couple of ways to go on this. One way would be to use a computer to solve $g^5(x)=x$. Another way must be doable without a calculator??
Find an orbit of period 5 for the chaotic logistic function $g(x)=4x(1-x)$.
Note: There's a couple of ways to go on this. One way would be to use a computer to solve $g^5(x)=x$. Another way must be doable without a calculator??
I utilized Mathematica in order to numerically solve for orbits of period 5 for the chaotic logistic function $g(x)=4x(1-x)$.
I first defined the function using the following code:
g[c_][x_] = 4 x - 4 x^2 + c;
g[n_, c_][x_] := Nest[g[c], x, n];
Having defined the function, I next solved for values of x where after 5 iterations, we have returned to the initial value.
per5Points = NSolve[g[5, 0][x] == x, x]
Finally, I extracted and tested one of the nonzero output solutions, taking it through two cycles to be quite certain of the periodic behavior.
p5 = x /. per5Points[[2]]
NestList[g[0], p5, 10]
Which produced the cycle $0.00903565, 0.035816, 0.138133, 0.476209, 0.997736, 0.00903565,
0.035816, 0.138133, 0.476209, 0.997736, 0.00903565$ , which does return to the initial value every 5th iteration.
The full list of solutions generated is:
{x -> 0.}, {x -> 0.00903565}, {x -> 0.010235}, {x -> 0.035816}, {x ->
0.0405211}, {x -> 0.0793732}, {x -> 0.0896183}, {x ->
0.138133}, {x -> 0.155517}, {x -> 0.209972}, {x -> 0.235524}, {x ->
0.292363}, {x -> 0.329169}, {x -> 0.377134}, {x -> 0.411438}, {x ->
0.445821}, {x -> 0.493}, {x -> 0.511038}, {x -> 0.578354}, {x ->
0.615607}, {x -> 0.662237}, {x -> 0.715041}, {x -> 0.716322}, {x ->
0.772761}, {x -> 0.826939}, {x -> 0.928658}, {x -> 0.941514}, {x ->
1.04274}, {x -> 1.04331}, {x -> 1.09213}, {x -> 1.10036}, {x ->
1.12518}
I think this solution by @SomeCallMeTim is awesome! I wonder if anyone else has a different approach? One that could be done by hand??
Using @SomeCallMeTim's values, we could draw a graph that looks like this
and imagine the cobweb plot superimposed onto it, but it's not obvious that a starting value of $0.00903565$ produces a period $5$ orbit.
To find a period 5 orbit of g, we can use the binary representation and choose values for c such that
$0_2.\overline{c_1c_2c_3c_4c_5}=0_2\overline{.00001}=\sum_{k=1}^{\infty}\frac{1}{32^k}=\sum_{k=0}^{\infty}\frac{1}{32^{k+1}}=\frac{1}{32}*\sum_{k=0}^{\infty}\frac{1}{32^{k}}=\frac{1}{32}
\frac{1}{1-\frac{1}{32}}=\frac{31}{32}\frac{32}{31}=\frac{1}{31}$.
So, the orbit containing $\frac{1}{31}$ is a period 5 orbit of g.
@Rick
Be sure to check your work. For example, if you think that $1/31$ is an orbit of period 5 for $g$, then compute 5 iterates and see where you are.
In reality, I think you've missed a crucial step.
I've added Example 7.9 to our notes on real iteration that might help.