Mandelbrot set

Initial image of a Mandelbrot set zoom sequence with continuously coloured environment

In mathematics, the Mandelbrot set, named after Benoît Mandelbrot, is a set of points in the complex plane, the boundary of which forms a fractal. Mathematically, the Mandelbrot set can be defined as the set of complex c-values for which the orbit of 0 under iteration of the complex quadratic polynomial x_n+1 = x_n² + c remains bounded.

Eg. c = 1 gives the sequence 0, 1, 2, 5, 26… which leads to infinity. As this sequence is unbounded, 1 is not an element of the Mandelbrot set.

On the other hand, c = i gives the sequence 0, i, (−1 + i), −i, (−1 + i), −i…, which is bounded, and so it belongs to the Mandelbrot set.

When computed and graphed on the complex plane, the Mandelbrot Set is seen to have an elaborate boundary, which does not simplify at any given magnification. This qualifies the boundary as a fractal.

The Mandelbrot set has become popular outside mathematics both for its aesthetic appeal and for being a complicated structure arising from a simple definition. Benoît Mandelbrot and others worked hard to communicate this area of mathematics to the public.

History

The Mandelbrot set has its place in complex dynamics, a field first investigated by the French mathematicians Pierre Fatou and Gaston Julia at the beginning of the 20th century. The first pictures of it were drawn in 1978 by Robert Brooks and Peter Matelski as part of a study of Kleinian Groups.

Mandelbrot studied the parameter space of quadratic polynomials in an article that appeared in 1980. The mathematical study of the Mandelbrot set really began with work by the mathematicians Adrien Douady and John H. Hubbard, who established many of its fundamental properties and named the set in honour of Mandelbrot.

The mathematicians Heinz-Otto Peitgen and Peter Richter became well-known for promoting the set with glossy photographs, books, and a touring gallery.

The cover article of the August 1985 Scientific American featured an image created by Mandelbrot, Peitgen, and Hubbard.

The work of Douady and Hubbard coincided with a huge increase in interest in complex dynamics and abstract mathematics, and the study of the Mandelbrot set has been a centerpiece of this field ever since. An exhaustive list of all the mathematicians who have contributed to the understanding of this set since then is beyond the scope of this article, but such a list would notably include Mikhail Lyubich, , Curt McMullen, John Milnor, Mitsuhiro Shishikura, and Jean-Christophe Yoccoz.

Formal definition

The Mandelbrot set $M$ is defined by a family of complex quadratic polynomials

$P_{c}:\mathbb {C} \to \mathbb {C}$

given by

$P_{c}:z\mapsto z^{2}+c$,

where $c$ is a complex parameter. For each $c$, one considers the behavior of the sequence $(0,P_{c}(0),P_{c}(P_{c}(0)),P_{c}(P_{c}(P_{c}(0))),\ldots )$ obtained by iterating $P_{c}(z)$ starting at critical point $z=0\,$, which either escapes to infinity or stays within a disk of some finite radius. The Mandelbrot set is defined as the set of all points $c$ such that the above sequence does not escape to infinity.

A mathematician's depiction of the Mandelbrot set M, a point c is coloured black if it belongs to the set, and white if not. Re[c] and Im[c] denote the real and imaginary parts of c.

More formally, if $P_{c}^{\circ n}(z)$ denotes the nth iterate of $P_{c}(z)$ (i.e. $P_{c}(z)$ composed with itself n times), the Mandelbrot set is the subset of the complex plane given by

$M=\left\{c\in \mathbb {C} :\sup _{n\in \mathbb {N} }|P_{c}^{\circ n}(0)|<\infty \right\}.$

Mathematically, the Mandelbrot set is just a set of complex numbers. A given complex number c either belongs to M or it does not. A picture of the Mandelbrot set can be made by colouring all the points $c$ which belong to M black, and all other points white. The more colourful pictures usually seen are generated by colouring points not in the set according to how quickly or slowly the sequence $|P_{c}^{\circ n}(0)|$ diverges to infinity. See the section on computer drawings below for more details.

The Mandelbrot set can also be defined as the connectedness locus of the family of polynomials $P_{c}(z)$. That is, it is the subset of the complex plane consisting of those parameters $c$ for which the Julia set of $P_{c}$ is connected.

Basic properties

The Mandelbrot set is a compact set, contained in the closed disk of radius 2 around the origin. In fact, a point $c$ belongs to the Mandelbrot set if and only if $|P_{c}^{\circ n}(0)|\leq 2$ for all $n\geq 0$. In other words, if the absolute value of $P_{c}^{\circ n}(0)$ ever becomes larger than 2, the sequence will escape to infinity.

Correspondence between the mandelbrot set and the logistic map

The intersection of $M$ with the real axis is precisely the interval $[-2,0.25]\,$. The parameters along this interval can be put in one-to-one correspondence with those of the real logistic family,

$z\mapsto \lambda z(1-z),\quad \lambda \in [1,4].\,$

The correspondence is given by

$c={\frac {1-(\lambda -1)^{2}}{4}}.$

In fact, this gives a correspondence between the entire parameter space of the logistic family and that of the Mandelbrot set.

The area of the Mandelbrot set is estimated to be 1.506 591 77 ± 0.000 000 08.

Douady and Hubbard have shown that the Mandelbrot set is connected. In fact, they constructed an explicit conformal isomorphism between the complement of the Mandelbrot set and the complement of the closed unit disk. Mandelbrot had originally conjectured that the Mandelbrot set is disconnected. This conjecture was based on computer pictures generated by programs which are unable to detect the thin filaments connecting different parts of $M$. Upon further experiments, he revised his conjecture, deciding that $M$ should be connected.

The dynamical formula for the uniformisation of the complement of the Mandelbrot set, arising from Douady and Hubbard's proof of the connectedness of $M$, gives rise to external rays of the Mandelbrot set. These rays can be used to study the Mandelbrot set in combinatorial terms and form the backbone of the Yoccoz parapuzzle.

The boundary of the Mandelbrot set is exactly the bifurcation locus of the quadratic family; that is, the set of parameters $c$ for which the dynamics changes abruptly under small changes of $c.$ It can be constructed as the limit set of a sequence of plane algebraic curves, the Mandelbrot curves, of the general type known as polynomial lemniscates. The Mandelbrot curves are defined by setting p₀=z, p_n=p_n-1²+z, and then interpreting the set of points |p_n(z)|=1 in the complex plane as a curve in the real Cartesian plane of degree 2ⁿ⁺¹ in x and y.

Other properties

The main cardioid and period bulbs

Periods of hyperbolic components

Upon looking at a picture of the Mandelbrot set, one immediately notices the large cardioid-shaped region in the center. This main cardioid is the region of parameters $c$ for which $P_{c}$ has an attracting fixed point. It consists of all parameters of the form

$c={\frac {1-(\mu -1)^{2}}{4}}$

for some $\mu \,$ in the open unit disk.

To the left of the main cardioid, attached to it at the point $c=-3/4$, a circular-shaped bulb is visible. This bulb consists of those parameters $c\,$ for which $P_{c}$ has an attracting cycle of period 2. This set of parameters is an actual circle, namely that of radius 1/4 around -1.

There are infinitely many other bulbs tangent to the main cardioid: for every rational number ${\frac {p}{q}}$, with p and q coprime, there is such a bulb that is tangent at the parameter

$c_{\frac {p}{q}}={\frac {1-\left(e^{2\pi i{\frac {p}{q}}}-1\right)^{2}}{4}}.$

Attracting cycle in 2/5-bulb plotted over Julia set (animation)

This bulb is called the ${\frac {p}{q}}$-bulb of the Mandelbrot set. It consists of parameters which have an attracting cycle of period $q$ and combinatorial rotation number ${\frac {p}{q}}$. More precisely, the $q$ periodic Fatou components containing the attracting cycle all touch at a common point (commonly called the $\alpha \,$-fixed point). If we label these components $U_{0},\dots ,U_{q-1}$ in counterclockwise orientation, then $P_{c}$ maps the component $U_{j}$ to the component $U_{j+p\,(\operatorname {mod} q)}$.

Attracting cycles and Julia sets for parameters in the 1/2, 3/7, 2/5, 1/3, 1/4, and 1/5 bulbs

The change of behavior occurring at $c_{\frac {p}{q}}$ is known as a bifurcation: the attracting fixed point "collides" with a repelling period q-cycle. As we pass through the bifurcation parameter into the ${\frac {p}{q}}$-bulb, the attracting fixed point turns into a repelling fixed point (the $\alpha$-fixed point), and the period q-cycle becomes attracting.

Hyperbolic components

All the bulbs we encountered in the previous section were interior components of the Mandelbrot set in which the maps $P_{c}\,$ have an attracting periodic cycle. Such components are called hyperbolic components.

It is conjectured that these are the only interior regions of $M$. This problem, known as density of hyperbolicity, may be the most important open problem in the field of complex dynamics. Hypothetical non-hyperbolic components of the Mandelbrot set are often referred to as "queer" components.

For real quadratic polynomials, this question was answered positively in the 1990s independently by Lyubich and by Graczyk and Świątek. (Note that hyperbolic components intersecting the real axis correspond exactly to periodic windows in the Feigenbaum diagram. So this result states that such windows exist near every parameter in the diagram.)

Not every hyperbolic component can be reached by a sequence of direct bifurcations from the main cardioid of the Mandelbrot set. However, such a component can be reached by a sequence of direct bifurcations from the main cardioid of a little Mandelbrot copy (see below).

Local connectivity

It is conjectured that the Mandelbrot set is locally connected. This famous conjecture is known as MLC (for Mandelbrot Locally Connected). By the work of Adrien Douady and John H. Hubbard, this conjecture would result in a simple abstract "pinched disk" model of the Mandelbrot set. In particular, it would imply the important hyperbolicity conjecture mentioned above.

The celebrated work of Jean-Christophe Yoccoz established local connectivity of the Mandelbrot set at all finitely-renormalizable parameters; that is, roughly speaking those which are contained only in finitely many small Mandelbrot copies. Since then, local connectivity has been proved at many other points of $M$, but the full conjecture is still open.

Self-similarity

Self similarity in the Mandelbrot set shown by zooming in on a round feature while panning in the negative-x direction. The display center pans from (−1, 0) to (−1.31, 0) while the view magnifies from 0.5 × 0.5 to 0.12 × 0.12 to approximate the Feigenbaum ratio $\delta$.

Self-similarity around Misiurewicz point −0.1011 + 0.9563i.

The Mandelbrot set is self-similar under magnification in the neighborhoods of the Misiurewicz points. It is also conjectured to be self-similar around generalized Feigenbaum points (e.g. −1.401155 or −0.1528 + 1.0397i), in the sense of converging to a limit set.

Quasi-self-similarity in the Mandelbrot set

The Mandelbrot set in general is not strictly self-similar but it is quasi-self-similar, as small slightly different versions of itself can be found at arbitrarily small scales.

The little copies of the Mandelbrot set are all slightly different, mostly because of the thin threads connecting them to the main body of the set.

Further results

The Hausdorff dimension of the boundary of the Mandelbrot set equals 2 as determined by a result of Mitsuhiro Shishikura. It is not known whether the boundary of the Mandelbrot set has positive planar Lebesgue measure.

In the Blum-Shub-Smale model of real computation, the Mandelbrot set is not computable, but its complement is computably enumerable. However, many simple objects (e.g., the graph of exponentiation) are also not computable in the BSS model. At present it is unknown whether the Mandelbrot set is computable in models of real computation based on computable analysis, which correspond more closely to the intuitive notion of "plotting the set by a computer." Hertling has shown that the Mandelbrot set is computable in this model if the hyperbolicity conjecture is true.

Relationship with Julia sets

An "embedded Julia set"

As a consequence of the definition of the Mandelbrot set, there is a close correspondence between the geometry of the Mandelbrot set at a given point and the structure of the corresponding Julia set.

This principle is exploited in virtually all deep results on the Mandelbrot set. For example, Shishikura proves that, for a dense set of parameters in the boundary of the Mandelbrot set, the Julia set has Hausdorff dimension two, and then transfers this information to the parameter plane. Similarly, Yoccoz first proves the local connectivity of Julia sets, before establishing it for the Mandelbrot set at the corresponding parameters. Adrien Douady phrases this principle as

Plough in the dynamical plane, and harvest in parameter space.

Geometry

cycle periods and antennae

Recall that, for every rational number ${\frac {p}{q}}$, where p and q are relatively prime, there is a hyperbolic component of period q bifurcating from the main cardioid. The part of the Mandelbrot set connected to the main cardioid at this bifurcation point is called the p/q-limb. Computer experiments suggest that the diameter of the limb tends to zero like ${\frac {1}{q^{2}}}$. The best current estimate known is the famous Yoccoz-inequality, which states that the size tends to zero like 1 / q.

A period-q limb will have q − 1 "antennae" at the top of its limb. We can thus determine the period of a given bulb by counting these antennas.

Image gallery of a zoom sequence

The Mandelbrot set shows more intricate detail the closer one looks or magnifies the image, usually called "zooming in". The following example of an image sequence zooming to a selected c value gives an impression of the infinite richness of different geometrical structures, and explains some of their typical rules.

The magnification of the last image relative to the first one is about 60,000,000,000 to 1. Relating to an ordinary monitor, it represents a section of a Mandelbrot set with a diameter of 20 million kilometres. Its border would show an inconceivable number of different fractal structures.

Start	Step 1	Step 2	Step 3	Step 4
Step 5	Step 6	Step 7	Step 8	Step 9
Step 10	Step 11	Step 12	Step 13	Step 14

Start: Mandelbrot set with continuously coloured environment.

1. Gap between the "head" and the "body" also called the "seahorse valley".
2. On the left double-spirals, on the right "seahorses".
3. "Seahorse" upside down, its "body" is composed by 25 "spokes" consisting of 2 groups of 12 "spokes" each and one "spoke" connecting to the main cardioid; these 2 groups can be attributed by some kind of metamorphosis to the 2 "fingers" of the "upper hand" of the Mandelbrot set, therefore, the number of "spokes" increases from one "seahorse" to the next by 2; the "hub" is a so-called Misiurewicz point; between the "upper part of the body" and the "tail" a distorted small copy of the Mandelbrot set called satellite may be recognized.
4. The central endpoint of the "seahorse tail" is also a Misiurewicz point.
5. Part of the "tail" — there is only one path consisting of the thin structures that leads through the whole "tail"; this zigzag path passes the "hubs" of the large objects with 25 "spokes" at the inner and outer border of the "tail"; it makes sure that the Mandelbrot set is a so-called simply connected set, which means there are no islands and no loop roads around a hole.
6. Satellite. The two "seahorse tails" are the beginning of a series of concentric crowns with the satellite in the center.
7. Each of these crowns consists of similar "seahorse tails"; their number increases with powers of 2, a typical phenomenon in the environment of satellites, the unique path to the spiral center mentioned in zoom step 5 passes the satellite from the groove of the cardioid to the top of the "antenna" on the "head".
8. "Antenna" of the satellite. Several satellites of second order may be recognized.
9. The "seahorse valley" of the satellite. All the structures from the image of zoom step 1 reappear.
10. Double-spirals and "seahorses" - unlike the image of zoom step 2 they have appendices consisting of structures like "seahorse tails"; this demonstrates the typical linking of n+1 different structures in the environment of satellites of the order n, here for the simplest case n=1.
11. Double-spirals with satellites of second order - analog to the "seahorses" the double-spirals may be interpreted as a metamorphosis of the "antenna".
12. In the outer part of the appendices islands of structures may be recognized; they have a shape like Julia sets J_c; the largest of them may be found in the center of the "double-hook" on the right side.
13. Part of the "double-hook".
14. At first sight, these islands seem to consist of infinitely many parts like Cantor sets, as is actually the case for the corresponding Julia set J_c. Here they are connected by tiny structures so that the whole represents a simply connected set. These tiny structures meet each other at a satellite in the center that is too small to be recognized at this magnification. The value of c for the corresponding J_c is not that of the image center but, relative to the main body of the Mandelbrot set, has the same position as the center of this image relative to the satellite shown in zoom step 7.

Generalizations

Multibrot set with d changing from 0 to 8

Sometimes the connectedness loci of families other than the quadratic family are also referred to as the Mandelbrot sets of these families.

The connectedness loci of the unicritical polynomial families ƒ_c = z^d + c for d > 2 are often called Multibrot sets.

Multibrot sets of degrees 3 and 4

For general families of holomorphic functions, the boundary of the Mandelbrot set generalizes to the bifurcation locus, which is a natural object to study even when the connectedness locus is not useful.

It is also possible to consider similar constructions in the study of non-analytic mappings. Of particular interest is the tricorn, the connectedness locus of the anti-holomorphic family

$z\mapsto {\bar {z}}^{2}+c.\,$

The tricorn (also sometimes called the Mandelbar set) was encountered by Milnor in his study of parameter slices of real cubic polynomials. It is not locally connected. This property is inherited by the connectedness locus of real cubic polynomials.

Computer drawings

Buddhabrot method

Still image of a movie of increasing magnification on 0.001643721971153 + 0.822467633298876i

Algorithms :

• Escape time algorithm
  • boolean version ( draws M-set and its exterior using 2 colours ) = Mandelbrot algorithm
  • discrete (integer) version = level set method ( LSM/M ); draws Mandelbrot set and colour bands in its exterior
  • continuous version
  • level curves version = draws lemniscates of Mandelbrot set = boundaries of Level Sets
  • decomposition of exterior of Mandelbrot set
• complex potential
  • Hubbard-Douady (real) potential of Mandelbrot set (CPM/M) - radial part of complex potential
  • external angle of Mandelbrot set - angular part of complex potential
  • abstract M-set
• Distance estimation method for Mandelbrot set
  • exterior distance estimation= Milnor algorithm (DEM/M)
  • interior distance estimation
• algorithm used to explore interior of Mandelbrot set
  • period of hyperbolic components
  • multiplier of periodic orbit ( internal rays(angle) and intenal radius )
  • bof61 and bof60

Every algorithm can be implemented in sequential or parallel version. Mirror symmetry can be used to speed-up calculations.

Escape time algorithm

The simplest algorithm for generating a representation of the Mandelbrot set is known as the "escape time" algorithm. A repeating calculation is performed for each x, y point in the plot area and based on the behaviour of that calculation, a colour is chosen for that pixel.

The x and y locations of each point are used as starting values in a repeating, or iterating calculation (described in detail below). The result of each iteration is used as the starting values for the next. The values are checked during each iteration to see if they have reached a critical 'escape' condition. If that condition is reached, the calculation is stopped, the pixel is drawn, and the next x, y point is examined. For some starting values, escape occurs quickly, after only a small number of iterations. For other starting values, it may take hundreds or thousands of iterations to escape. For values within the Mandelbrot set, escape will never occur. The programmer or user must choose how much iteration, or 'depth,' they wish to examine. The higher the maximum number of iterations, the more detail and subtlety emerge in the final image, but the longer time it will take to calculate the picture.

The colour of each point represents how quickly the values reached the escape point. Often black is used to show values that fail to escape before the iteration limit, and gradually brighter colours are used for points that escape. This gives a visual representation of how many cycles were required before reaching the escape condition.

For programmers

The definition of the Mandelbrot set, together with its basic properties, suggests a simple algorithm for drawing a picture of the Mandelbrot set. The region of the complex plane we are considering is subdivided into a certain number of pixels. To colour any such pixel, let $c\,$ be the midpoint of that pixel. We now iterate the critical value $c\,$ under $P_{c}\,$, checking at each step whether the orbit point has modulus larger than 2.

If this is the case, we know that the midpoint does not belong to the Mandelbrot set, and we colour our pixel. (Either we colour it white to get the simple mathematical image or colour it according to the number of iterations used to get the well-known colourful images). Otherwise, we keep iterating for a certain (large, but fixed) number of steps, after which we decide that our parameter is "probably" in the Mandelbrot set, or at least very close to it, and colour the pixel black.

In pseudocode, this algorithm would look as follows.

For each pixel on the screen do:
{
  x = x0 = x co-ordinate of pixel
  y = y0 = y co-ordinate of pixel

  iteration = 0
  max_iteration = 1000
 
  while ( x*x + y*y <= (2*2)  AND  iteration < max_iteration ) 
  {

    xtemp = x*x - y*y + x0
    y = 2*x*y + y0

    x = xtemp

    iteration = iteration + 1
  }
 
  if ( iteration == max_iteration ) 
  then 
    colour = black
  else 
    colour = iteration

  plot(x0,y0,colour)
}

where, relating the pseudocode to $c\,,z\,$ and $P_{c}\,$:

• $z=x+iy$
• $z^{2}=x^{2}+i2xy-y^{2}$
• $c=x0+iy0$

and so, as can be seen in the pseudocode in the computation of x and y:

• $x=Re(z^{2}+c)=x^{2}-y^{2}+x0$ and y = $Im(z^{2}+c)=2xy+y0$

To get colourful images of the set, the assignment of a colour to each value of the number of executed iterations can be made using one of a variety of functions (linear, exponential, etc). One practical way to do it, without slowing down the calculations, is to use the number of executed iterations as an entry to a look-up colour palette table initialized at startup. If the colour table has, for instance, 500 entries, then you can use n mod 500, where n is the number of iterations, to select the colour to use. You can initialize the colour palette matrix in various different ways, depending on what special feature of the escape behavior you want to emphasize graphically.

Continuous (smooth) coloring

File:Escape Time Algorithm.png

This image was rendered with the Escape Time Algorithm. Notice the very obvious "bands" of color.

This image was rendered with the Normalized Iteration Count Algorithm. Notice the bands of color have been replaced by a smooth gradient.

Another example of the Normalized Iteration Count Algorithm. Notice that there is no banding effect; all of the colors flow into each other. Also, the colors take on the same pattern that would be observed if the Escape Time Algorithm was used.

The Escape Time Algorithm is popular for its simplicity. However, it creates bands of color, which, as a type of aliasing can detract from an image's aesthetic value. This can be improved using the Normalized Iteration Count Algorithm, which provides a smooth transition of colors between iterations. The equation is

$n+{\frac {\ln(2\ln(B))-\ln(\ln(|z|))}{\ln(P)}}$

where n is the number of iterations for z, B is the bailout radius (it is normally 2 for a Mandelbrot set, but it can be changed), and P is the power for which z is raised to in the Mandelbrot set equation (z_n+1 = zn^P + c, P is generally 2). Another equation for this is

$n+1-{\frac {\ln(\ln(|z|))}{\ln(2)}}.$

Note that this new equation is simpler than the first, but it only works for Mandelbrot sets with a bailout radius of 2 and a power of 2.
While this algorithm is relatively simple to implement (using either equation), there are a few things that need to be taken into consideration. First, the two equations return a continuous stream of numbers. However, it is up to you to decide on how the return values will be converted into a color. Some type of method for casting these numbers onto a gradient should be developed. Second, it is recommended that a few extra iterations are done so that z can grow. If you stop iterating as soon as z escapes, there is the possibility that the smoothing algorithm will not work.

Distance estimates

One can compute distance from point c ( in exterior or interior ) to nearest point on the boundary of Mandelbrot set.

Exterior distance estimation

The proof of the connectedness of the Mandelbrot set in fact gives a formula for the uniformizing map of the complement of $M$ (and the derivative of this map). By the Koebe 1/4 theorem, one can then estimate the distance between the mid-point of our pixel and the Mandelbrot set up to a factor of 4.

In other words, provided that the maximal number of iterations is sufficiently high, one obtains a picture of the Mandelbrot set with the following properties:

1. Every pixel which contains a point of the Mandelbrot set is coloured black.
2. Every pixel which is colored black is close to the Mandelbrot set.

Exterior distance estimate may be used to color whole complement of Mandelbrot set

The distance estimate of a pixel c (a complex number) from the Mandelbrot set is given by

$b=\lim _{n\to \infty }2\cdot \ln(\mid {P_{c}^{\circ n}(c)}\mid )\cdot {\frac {\mid {P_{c}^{\circ n}(c)}\mid }{\mid {\frac {\partial }{\partial {c}}}P_{c}^{\circ n}(c)\mid }}$

where

• $P_{c}(z)\,$ stands for complex quadratic polynomial
• $P_{c}^{\circ n}(c)$ stands for n iterations of $P_{c}(z)\to z$ or $z^{2}+c\to z$, starting with z=c: $P_{c}^{\circ 0}(c)=c$, $P_{c}^{\circ n+1}(c)=P_{c}^{\circ n}(c)^{2}+c$;
• ${\frac {\partial }{\partial {c}}}P_{c}^{\circ n}(c)$ is the derivative of $P_{c}^{\circ n}(c)$ with respect to c. This derivative can be found by starting with ${\frac {\partial }{\partial {c}}}P_{c}^{\circ 0}(c)=1$ and then ${\frac {\partial }{\partial {c}}}P_{c}^{\circ n+1}(c)=2\cdot {}P_{c}^{\circ n}(c)\cdot {\frac {\partial }{\partial {c}}}P_{c}^{\circ n}(c)+1$. This can easily be verified by using the chain rule for the derivative.

From a mathematician's point of view, this formula only works in limit where n goes to infinity, but very reasonable estimates can be found with just a few additional iterations after the main loop exits.

Once b is found, by the Koebe 1/4-theorem, we know there's no point of the Mandelbrot set with distance from c smaller than b/4.

Interior distance estimation

Pixels colored according to the estimated interior distance

It is also possible to estimate the distance of a limitly periodic (i.e., inner) point to the boundary of the Mandelbrot set. The estimate is given by

$b={\frac {1-\mid {{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})}\mid ^{2}}{\mid {{\frac {\partial }{\partial {c}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})+{\frac {\partial }{\partial {z}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0}){\frac {{\frac {\partial }{\partial {c}}}P_{c}^{\circ p}(z_{0})}{1-{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})}}}\mid }}$

where

• $p$ is the period,
• $c$ is the point to be estimated,
• $P_{c}(z)$ is the complex quadratic polynomial $P_{c}(z)=z^{2}+c$
• $P_{c}^{\circ p}(z_{0})$ is $p$ compositions of $P_{c}(z)\to z$, starting with $P_{c}^{\circ 0}(z)=z_{0}$
• $z_{0}$ is any of the $p$ points that make the attractor of the iterations of $P_{c}(z)\to z$ starting with $P_{c}^{\circ 0}(z)=c$; $z_{0}$ satisfies $z_{0}=P_{c}^{\circ p}(z_{0})$,
• ${\frac {\partial }{\partial {c}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})$, ${\frac {\partial }{\partial {z}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})$, ${\frac {\partial }{\partial {c}}}P_{c}^{\circ p}(z_{0})$ and ${\frac {\partial }{\partial {z}}}P_{c}^{\circ p}(z_{0})$ are various derivatives of $P_{c}^{\circ p}(z)$, evaluated at $z_{0}$.

Given $p$ and $z_{0}$, $P_{c}^{\circ p}(z_{0})$ and its derivatives can be evaluated by:

• ${\begin{aligned}{\frac {\partial }{\partial {c}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ 0}(z_{0})=&0\\{\frac {\partial }{\partial {z}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ 0}(z_{0})=&0\\{\frac {\partial }{\partial {c}}}P_{c}^{\circ 0}(z_{0})=&0\\{\frac {\partial }{\partial {z}}}P_{c}^{\circ 0}(z_{0})=&1\\P_{c}^{\circ 0}(z_{0})=&z_{0}\end{aligned}}$
• ${\begin{aligned}{\frac {\partial }{\partial {c}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ n+1}(z_{0})&=2\cdot ({\frac {\partial }{\partial {c}}}P_{c}^{\circ n}(z_{0})\cdot {\frac {\partial }{\partial {z}}}P_{c}^{\circ n}(z_{0})+P_{c}^{\circ n}(z_{0})\cdot {\frac {\partial }{\partial {c}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ n}(z_{0}))\\{\frac {\partial }{\partial {z}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ n+1}(z_{0})&=2\cdot ({\frac {\partial }{\partial {z}}}P_{c}^{\circ n}(z_{0})^{2}+P_{c}^{\circ n}(z_{0})\cdot {\frac {\partial }{\partial {z}}}{\frac {\partial }{\partial {z}}}P_{c}^{\circ n}(z_{0}))\\{\frac {\partial }{\partial {c}}}P_{c}^{\circ n+1}(z_{0})&=2\cdot P_{c}^{\circ n}(z_{0})\cdot {\frac {\partial }{\partial {c}}}P_{c}^{\circ n}(z_{0})+1\\{\frac {\partial }{\partial {z}}}P_{c}^{\circ n+1}(z_{0})&=2\cdot P_{c}^{\circ n}(z_{0})\cdot {\frac {\partial }{\partial {z}}}P_{c}^{\circ n}(z_{0})\\P_{c}^{\circ n+1}(z_{0})&=P_{c}^{\circ n}(z_{0})^{2}+c\end{aligned}}$.

Analogous to the exterior case, once b is found, we know that all points within the distance of b/4 from c are inside the Mandelbrot set.

There are two practical problems with the interior distance estimate: first, we need to find $z_{0}$ precisely, and second, we need to find $p$ precisely. The problem with $z_{0}$ is that the convergence to $z_{0}$ by iterating $P_{c}(z)$ requires, theoretically, an infinite number of operations. The problem with period is that, sometimes, due to rounding errors, a period is falsely identified to be an integer multiple of the real period (e.g., a period of 86 is detected, while the real period is only 43=86/2). In such case, the distance is overestimated, i.e., the reported radius could contain points outside the Mandelbrot set.

3D view : smallest absolute value of the orbit of the interior points of the Mandelbrot set

Optimizations

One way to improve calculations is to find out beforehand whether the given point lies within the cardioid or in the period-2 bulb.

To prevent having to do huge numbers of iterations for other points in the set, one can do "periodicity checking"—which means check if a point reached in iterating a pixel has been reached before. If so, the pixel cannot diverge, and must be in the set. This is most relevant for fixed-point calculations, where there is a relatively high chance of such periodicity—a full floating-point (or higher-accuracy) implementation would rarely go into such a period.

Periodicity checking is, of course, a trade-off: The need to remember points costs memory and data management instructions, whereas it saves computational instructions.

In popular culture

The folk-rock musician Jonathan Coulton wrote a song about the Mandelbrot Set, with the same title. Pisut Wisessing, a Cornell University student, created a music video for the song with chalkboard animation.

References

Further reading

• Lennart Carleson and Theodore W. Gamelin, Complex Dynamics, Springer 1993, ISBN 0-387-97942-5
• John W. Milnor, Dynamics in One Complex Variable (Third Edition), Annals of Mathematics Studies 160, Princeton University Press 2006, ISBN 0-691-12488-4 (First appeared in 1990 as a Stony Brook IMS Preprint, available as arXiV:math.DS/9201272.)
• Nigel Lesmoir-Gordon. "The Colours of Infinity: The Beauty, The Power and the Sense of Fractals." ISBN 1-904555-05-5 (The book comes with a related DVD of the Arthur C. Clarke documentary introduction to the fractal concept and the Mandelbrot set.

See also

• Burning Ship fractal
• Newton fractal
• Pickover stalk
• Orbit portrait
• XaoS

External links

• Template:Dmoz
• The Mandelbrot Set and Julia Sets by Michael Frame, Benoit Mandelbrot, and Nial Neger