Órbita crítica en caso de disco Siegel: media dorada
Vista 3D de la órbita crítica que tiende a un punto fijo parabólico
distancia entre puntos de órbita crítica en caso de atraer un punto fijo
Código fuente de Maxima
/ * este es un archivo por lotes para Maxima. Guárdelo con la extensión mac. Calcula la órbita de Z0 bajo f (z, c), encuentra el período y dibuja la órbita * /
c: -0,2-0,7 *% i;
f (z, c): = z * z + c;
z0: 0;
iMax: 100;
eps: 0,01;
/* primer punto */
z: z0;
órbita: [z];
/ * ------------------- calcular la órbita de avance -------------------------- --- * /
para i: 1 a iMax paso 1 hacer
cuadra
( z: f (z, c), órbita: endcons (z, órbita) / * disp (rectform (órbita [i])) * / );
/ * ----------- encontrar período de órbita --------------------------------- ----- * /
IsEqual (c1, c2, eps): = if abs (realpart (c1) -realpart (c2)) <= eps y abs (imagpart (c1) -imagpart (c2)) <= eps entonces cierto más falso;
GivePeriod (órbita, eps): =
cuadra ( período: 0, iLast: longitud (órbita), yo: iLast, cuadra ( lazo, yo: i-1, si no es IsEqual (órbita [iLast], órbita [i], eps) entonces vaya (bucle), período: iLast-i ) );
período: GivePeriod (órbita, eps);
imprimir (punto);
/ * -------------- dibujar órbita en la pantalla ----------------------------- * /
/ * guarda los valores z en 2 listas * /
xx: makelist (parte real (z0), i, 1, 1); / * lista de re (z) * /
yy: makelist (imagpart (z0), i, 1, 1); / * lista de im (z) * /
para i: 2 a través de la longitud (órbita) paso 1 hacer cuadra ( xx: contras (parte real (órbita [i]), xx), yy: contras (imagpart (órbita [i]), yy) );
cargar (dibujar);
dibujar2d file_name = "órbita", terminal = 'png, rango y = [-1,1], rango x = [-1,1], title = concat ("Período", cadena (período), "órbita de z0 =", cadena (z0), "para c =", cadena (c), "f (z, c): = z * z + c "), key = "punto de órbita", xlabel = "Z.re", ylabel = "Z.im", tipo_punto = círculo_lleno, tamaño_punto = 0.5, color = rojo, puntos (xx, yy) );
Licencia
Yo, el titular de los derechos de autor de este trabajo, lo publico bajo la siguiente licencia:
compartir - copiar, distribuir y transmitir el trabajo
remezclar - adaptar el trabajo
Bajo las siguientes condiciones:
atribución : debe otorgar el crédito correspondiente, proporcionar un enlace a la licencia e indicar si se realizaron cambios. Puede hacerlo de cualquier manera razonable, pero no de ninguna manera que sugiera que el licenciante lo respalda a usted o su uso.
compartir por igual : si remezcla, transforma o construye sobre el material, debe distribuir sus contribuciones bajo la misma licencia o una licencia compatible que la original.