Özet

AçıklamaCr orbit 3.png	critical orbit. Here one can see that point z=0 tends to 3 point attracting cycle.
Kaynak	self-made with Maxima program ( code below)   Bu PNG grafik Gnuplot ile oluşturuldu.
Yazar	Adam majewski

See also

• 
  Critical orbit in Siegel disc case : Golden mean
• 
  3D view of critical orbit tending to parabolic fixed point
• 
  distance between points of critical orbit in case of attracting fixed point

Maxima source code

/* this is batch file for Maxima. Save it with mac extension.
   It computes orbit of Z0 under f(z,c), finds period and draws orbit */
c:-0.2-0.7*%i;
f(z,c):=z*z+c;
z0:0;
iMax:100;
eps:0.01;
/* first point */
z:z0;
orbit:[z];
/* ------------------- compute forward orbit -----------------------------*/
for i:1 thru iMax step 1 do
block
(
 z:f(z,c),
 orbit:endcons(z,orbit)
 /* disp(rectform(orbit[i])) */
 );
/* -----------  find period of orbit --------------------------------------*/
IsEqual(c1,c2,eps):=
   if  abs(realpart(c1)-realpart(c2))<=eps and abs(imagpart(c1)-imagpart(c2))<=eps
    then true
    else false;
GivePeriod(orbit,eps):=
block
 (
   period:0,
   iLast:length(orbit),
   i:iLast,
   block
    (
     loop,
     i:i-1,
     if not IsEqual(orbit[iLast],orbit[i],eps) then go(loop),
     period:iLast-i	
    )
  );
period:GivePeriod(orbit,eps);
print(period);
/*-------------- draw orbit on the screen ----------------------------- */
/* save the z values to 2 lists */
xx:makelist (realpart(z0), i, 1, 1); /* list of re(z) */
yy:makelist (imagpart(z0), i, 1, 1); /* list of im(z) */
for i:2 thru length(orbit) step 1 do
 block
 (
 xx:cons(realpart(orbit[i]),xx), 
 yy:cons(imagpart(orbit[i]),yy)
 );
load(draw);
draw2d(
   file_name = "orbit",
   terminal  = 'png,
   yrange = [-1,1],
   xrange = [-1,1],
   title= concat("Period ",string(period)," orbit of z0=",string(z0)," for c=",string(c)," f(z,c):=z*z+c "),
   key = "point of orbit",
   xlabel     = "Z.re ",
   ylabel     = "Z.im",
   point_type    = filled_circle,
   point_size    = 0.5,
   color         = red,
   points(xx,yy)
 );

Lisanslama

Ben, bu işin telif sahibi, burada işi aşağıdaki lisans altında yayımlıyorum:

Bu dosya, Creative Commons Atıf-Benzer Paylaşım 3.0 Taşınmamış lisansı ile lisanslanmıştır

Şu seçeneklerde özgürsünüz:

• paylaşım – eser paylaşımı, dağıtımı ve iletimi
• içeriği değiştirip uyarlama – eser adaptasyonu

Aşağıdaki koşullar geçerli olacaktır:

• atıf – Esere yazar veya lisans sahibi tarafından belirtilen (ancak sizi ya da eseri kullanımınızı desteklediklerini ileri sürmeyecek bir) şekilde atıfta bulunmalısınız.
• benzer paylaşım – Maddeyi yeniden karıştırır, dönüştürür veya inşa ederseniz, katkılarınızı orijinal olarak aynı veya uyumlu lisans altında dağıtmanız gerekir.

https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue