Desde $\cos (2x) = 2 \cos^2 (x) - 1$, me pregunto acerca de la iteración $(x,y) \mapsto (2x\mod 2\pi, 2y^2-1)$. Se que convergen para el coseno del gráfico? Lo he probado en mathematica y tengo unos cuantos puntos, si me estafan y darle un valor inicial.
x = {1, 0.540302306};
step[{a_, b_}] := {Mod[2 a, 2 Pi], 2 b^2 - 1}
Table[(x = step[x]; x), {k, 1, 30}];
ListPlot[%, AspectRatio -> Automatic]
Aquí una variante basada en la suma de los ángulos de la fórmula da buena convergenc a la gráfica.
x = {1, 0.540302306};
step[{a_, b_}] := (c = Mod[a + 1, 2 Pi]; {c ,
b Cos[1] - If[c < Pi, 1, -1] Sin[1] Sqrt[1 - b^2]})
Table[(x = step[x]; x), {k, 1, 500}];
ListPlot[%, AspectRatio -> Automatic]
Estos dos algoritmos tienen diferencia de las tolerancias de error. En el primer algoritmo, si empiezo con un mal valor inicial como {1, 0.4} puedo obtener una caótica "en la nube", pero en el segundo caso, todavía converge al coseno del gráfico.
Supongo que me gustaría saber, ¿mi original iteración tiene cualquier otro puntos fijos además el coseno del gráfico?