Archivo: Cr orbit 3.png


Este es un archivo de Wikimedia Commons
De Wikipedia, la enciclopedia libre
Saltar a navegación Saltar a búsqueda

Cr_orbit_3.png (640 × 480 píxeles, tamaño de archivo: 4 KB, tipo MIME: image / png )

Resumen

Ver también

  • Ó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:

Subtítulos

Add a one-line explanation of what this file represents

Items portrayed in this file

depicts

Historial del archivo

Haga clic en una fecha / hora para ver el archivo tal como apareció en ese momento.

Las siguientes páginas de la Wikipedia en inglés utilizan este archivo (no se enumeran las páginas de otros proyectos):

Uso de archivos global