gnuplot / fractal / selfsq (E)

Julia Set (Self-Squared Fractal)

Same as the Mandelbrot set. We calculate the recurrence relation
with a complex variable A, and look for a set of numbers
with which the value of |z(n)|^2 does not diverge for
a large number of n. The differences from the Mandelbrot set
are that the initial value of the recurrence relation is the (x,y)
coordinate on the complex plane, and the number of A is
arbitrary.

z(0)   = X + iY
z(n+1) = z(n)*z(n) + A

The function calculate the function above is just the same as our
mandel function defined in the
previous section
. The complex value A is given outside
the function, and the (X,Y) coordinate is used as the initial
value. Here we employ A=-0.37-0.612 i.

gnuplot> set xrange [-0.5:0.5]
gnuplot> set yrange [-0.5:0.5]
gnuplot> set logscale z
gnuplot> set isosample 50
gnuplot> set hidden3d
gnuplot> set contour
gnuplot> a= -0.37
gnuplot> b= -0.612
gnuplot> splot mandel(a,b,complex(x,y),0) notitle
selfsq1

You can make a great variety of images by changing the complex
constant A. In fact, a very small change in the constant
result in a completely different picture. For example, the value of
A was -0.37 -0.612 i in the figure above, but if you
change the imaginary part into 0.6, you get:

selfsq2


up

The following CGs of the self-squared fractal are made with the
parameters above.

selfsquared1
selfsquared2
up