Distance between points along orbit for various z-points for fc(z)=z*z -0.390540870218399 -0.586787907346969*i
Long description
Maxima code
Here is batch file for Maxima:
c:-0.71+0.1*%i;
f(z,c):=z*z+c;
iMax:300;
z0:0+0*%i;
dist(z_prev,z_next):=
block(
dr:realpart(z_next)-realpart(z_prev),
di:imagpart(z_next)-imagpart(z_prev),
d:sqrt(dr*dr+di*di),
return(d)
);
/* */
z1:f(z0,c);
d1:dist(z0,z1);
/* save the z values to 2 lists */
xx:makelist (1, i, 1, 1); /* list of numbers of iterations */
yy:makelist (d1, i, 1, 1); /* list of distance(z(n),z(n+1)) */
/* */
z_prev:z1;
for i:2 thru iMax step 1 do
block
(
z_next:f(z_prev,c),
d:dist(z_prev,z_next),
xx:cons(i,xx),
yy:cons(d,yy),
z_prev:z_next
);
load(draw);
draw2d(
file_name = "distance",
terminal = 'png,
yrange = [0,0.8],
xrange = [0,iMax],
pic_width = 1000,
pic_height = 1000,
title= "distance between points zn and z(n+1) where z(n+1)=f(zn)= zn*zn+c and z0=0 is critical point",
key = concat("distance for c = ",string(c)),
xlabel = "number of iteration",
ylabel = "distance",
point_type = filled_circle,
point_size = 0.5,
color = red,
points(xx,yy)
);
רישיון
אני, בעל זכויות היוצרים על היצירה הזאת, מפרסם אותה בזאת תחת הרישיונות הבאים:
מוענקת בכך הרשות להעתיק, להפיץ או לשנות את המסמך הזה, לפי תנאי הרישיון לשימוש חופשי במסמכים של גנו, גרסה 1.2 או כל גרסה מאוחרת יותר שתפורסם על־ידי המוסד לתוכנה חופשית; ללא פרקים קבועים, ללא טקסט עטיפה קדמית וללא טקסט עטיפה אחורית. עותק של הרישיון כלול בפרק שכותרתו הרישיון לשימוש חופשי במסמכים של גנו.http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue
ייחוס – יש לתת ייחוס הולם, לתת קישור לרישיון, ולציין אם נעשו שינויים. אפשר לעשות את זה בכל צורה סבירה, אבל לא בשום צורה שמשתמע ממנה שמעניק הרישיון תומך בך או בשימוש שלך.
שיתוף זהה – יצירת רמיקס, שינוי או בנייה על סמך החומר הזה, תטיל עליך חובה להפיץ את התרומות שלך לפי תנאי רישיון זהה או תואם למקור.
{{Information |Description=Distance between point of orbit of discrete dynamical system |Source=self-made with Maxima |Date= |Author= Adam majewski |Permission= |other_versions= }}