# Mandelbrot-Menge

Die Mandelbrot-Menge, benannt nach Benoît Mandelbrot, ist die Menge der komplexen Zahlen 
  
    
      
        c
      
    
    {\displaystyle c}
  
, für welche die durch die iterative Vorschrift 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  
 mit dem Anfangswert 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
 definierte Folge 
  
    
      
        
          z
          
            0
          
        
        ,
        
          z
          
            1
          
        
        ,
        
          z
          
            2
          
        
        ,
        
          z
          
            3
          
        
        ,
        …
      
    
    {\displaystyle z_{0},z_{1},z_{2},z_{3},\ldots }
  
 endlich bleibt, d. h. beschränkt ist.
Interpretiert man die Mandelbrot-Menge (eine Teilmenge der Gaußschen Zahlenebenen) als geometrische Figur, so ergibt sie ein Fraktal, das im allgemeinen Sprachgebrauch oft Apfelmännchen genannt wird.
Bilder berechnet man, indem man jedem Pixel 
  
    
      
        (
        x
        ,
        y
        )
      
    
    {\displaystyle (x,y)}
  
 eines Bildes eine komplexe Zahl zuordnet 
  
    
      
        (
        c
        =
        
          c
          
            0
          
        
        +
        a
        
        ⋅
        
        x
        +
        b
        i
        
        ⋅
        
        y
        )
      
    
    {\displaystyle (c=c_{0}+a\!\cdot \!x+bi\!\cdot \!y)}
  
 und beginnend mit 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
 untersucht, ob und wann die Iterationen anfangen, zu „explodieren“. Bleiben die Werte klein, wird das Pixel häufig schwarz gefärbt, kommt es zu einer „Explosion“ der Zahlenwerte, wird die Anzahl der dafür notwendigen Iterationen als Farbe kodiert.
Die ersten mit einem Computer generierten Darstellungen wurden 1978 von Robert W. Brooks und Peter Matelski vorgestellt. 1980 veröffentlichte Benoît Mandelbrot eine Arbeit über das Thema. Später wurde sie von Adrien Douady und John Hamal Hubbard in einer Reihe grundlegender mathematischer Arbeiten systematisch untersucht. Die mathematischen Grundlagen dafür wurden bereits 1905 von dem französischen Mathematiker Pierre Fatou erarbeitet.


## Definition
### Definition über Rekursion
Die Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 ist die Menge aller komplexen Zahlen 
  
    
      
        c
      
    
    {\displaystyle c}
  
, für welche die rekursiv definierte Folge komplexer Zahlen 
  
    
      
        
          z
          
            0
          
        
        ,
        
          z
          
            1
          
        
        ,
        
          z
          
            2
          
        
        ,
        …
      
    
    {\displaystyle z_{0},z_{1},z_{2},\dotsc }
  
 mit dem Anfangsglied

  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  

und dem Bildungsgesetz

  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  

beschränkt bleibt. Das heißt, eine komplexe Zahl 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ist Element der Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
, wenn die Beträge der mit diesem 
  
    
      
        c
      
    
    {\displaystyle c}
  
 berechneten 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 nicht über jede Grenze wachsen, unabhängig davon, wie groß 
  
    
      
        n
      
    
    {\displaystyle n}
  
 wird. Dies lässt sich wie folgt schreiben:

  
    
      
        
          z
          
            0
          
        
        =
        0
        ,
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
        ;
        
        
        c
        ∈
        
          M
        
        
        ⟺
        
        ∃
        g
        ∈
        
          R
        
        :
        
          lim sup
          
            n
            →
            ∞
          
        
        
          |
        
        
          z
          
            n
          
        
        
          |
        
        ≤
        g
        .
      
    
    {\displaystyle z_{0}=0,z_{n+1}=z_{n}^{2}+c;\;\;c\in \mathbb {M} \iff \exists g\in \mathbb {R} :\limsup _{n\to \infty }|z_{n}|\leq g.}
  

Es lässt sich zeigen, dass der Betrag der 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 über jede Grenze wächst, sobald ein 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 mit 
  
    
      
        
          |
        
        
          z
          
            n
          
        
        
          |
        
        >
        2
      
    
    {\displaystyle |z_{n}|>2}
  
 auftritt, somit ist diese Definition gleichbedeutend mit:

  
    
      
        
          z
          
            0
          
        
        =
        0
        ,
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
        ;
        
        
        c
        ∈
        
          M
        
        
        ⟺
        
        
          lim sup
          
            n
            →
            ∞
          
        
        
          |
        
        
          z
          
            n
          
        
        
          |
        
        ≤
        2.
      
    
    {\displaystyle z_{0}=0,z_{n+1}=z_{n}^{2}+c;\;\;c\in \mathbb {M} \iff \limsup _{n\to \infty }|z_{n}|\leq 2.}
  


### Definition über komplexe quadratische Polynome
Die Mandelbrot-Menge lässt sich auch über komplexe quadratische Polynome beschreiben:

  
    
      
        
          P
          
            c
          
        
        :
        
          C
        
        →
        
          C
        
        ,
        
        z
        ↦
        
          z
          
            2
          
        
        +
        c
      
    
    {\displaystyle P_{c}:\mathbb {C} \to \mathbb {C} ,\,z\mapsto z^{2}+c}
  

mit einem komplexen Parameter 
  
    
      
        c
      
    
    {\displaystyle c}
  
. Für jedes 
  
    
      
        c
      
    
    {\displaystyle c}
  
 wird die Folge

  
    
      
        (
        
          P
          
            c
          
          
            0
          
        
        (
        0
        )
        ,
        
        
          P
          
            c
          
          
            1
          
        
        (
        0
        )
        ,
        
        
          P
          
            c
          
          
            2
          
        
        (
        0
        )
        ,
        
        ⋯
        )
        =
        (
        
        
          P
          
            c
          
          
            n
          
        
        (
        0
        )
        
        
          )
          
            n
            ∈
            
              N
            
          
        
      
    
    {\displaystyle (P_{c}^{0}(0),\;P_{c}^{1}(0),\,P_{c}^{2}(0),\,\dotsb )=(\,P_{c}^{n}(0)\,)_{n\in \mathbb {N} }}
  

iterativ berechnet, wobei 
  
    
      
        
          P
          
            c
          
          
            n
          
        
      
    
    {\displaystyle P_{c}^{n}}
  
 die 
  
    
      
        n
      
    
    {\displaystyle n}
  
-fache Hintereinanderausführung der Iteration bedeutet, also

  
    
      
        
          P
          
            c
          
          
            0
          
        
        (
        z
        )
        =
        z
      
    
    {\displaystyle P_{c}^{0}(z)=z}
  

  
    
      
        
          P
          
            c
          
          
            n
            +
            1
          
        
        (
        z
        )
        =
        
          P
          
            c
          
        
        (
        
          P
          
            c
          
          
            n
          
        
        (
        z
        )
        )
        ,
        
        n
        ∈
        
          N
        
        .
      
    
    {\displaystyle P_{c}^{n+1}(z)=P_{c}(P_{c}^{n}(z)),\;n\in \mathbb {N} .}
  

In Abhängigkeit vom Wert des Parameters 
  
    
      
        c
      
    
    {\displaystyle c}
  
 wächst diese Folge dann entweder unbeschränkt, sodass also 
  
    
      
        c
      
    
    {\displaystyle c}
  
 kein Element der Mandelbrot-Menge ist, oder sie verbleibt innerhalb eines Bereichs um den Ursprung der Zahlenebene, und 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ist Element der Mandelbrot-Menge.
Die Mandelbrot-Menge ist eine Untermenge der komplexen Zahlen mit der Definition:

  
    
      
        
          M
        
        =
        {
        c
        ∈
        
          C
        
        
        
          |
        
        
        ∃
        s
        ∈
        
          R
        
        :
        ∀
        n
        ∈
        
          N
        
        :
        
          |
        
        
          P
          
            c
          
          
            n
          
        
        (
        0
        )
        
          |
        
        ≤
        s
        }
      
    
    {\displaystyle \mathbb {M} =\lbrace c\in \mathbb {C} \;|\;\exists s\in \mathbb {R} :\forall n\in \mathbb {N} :|P_{c}^{n}(0)|\leq s\rbrace }
  

oder gleichbedeutend:

  
    
      
        
          M
        
        =
        {
        c
        ∈
        
          C
        
        
        
          |
        
        
        ∀
        n
        ∈
        
          N
        
        :
        
          |
        
        
          P
          
            c
          
          
            n
          
        
        (
        0
        )
        
          |
        
        ≤
        2
        }
        .
      
    
    {\displaystyle \mathbb {M} =\lbrace c\in \mathbb {C} \;|\;\forall n\in \mathbb {N} :|P_{c}^{n}(0)|\leq 2\rbrace .}
  

Zur Erläuterung werden einige Eigenschaften und Beispiele angeführt:

Aufgrund der zuvor beschriebenen Feststellung kann 
  
    
      
        s
        =
        2
      
    
    {\displaystyle s=2}
  
 gesetzt werden. Dabei gibt der Wert 
  
    
      
        s
        =
        2
      
    
    {\displaystyle s=2}
  
 den Radius um den Ursprung an, innerhalb dessen ein Element von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 liegen kann. Außerhalb dieses Kreises sind keine Elemente von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 zu finden.
Wegen der Betragsfunktion ist 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 symmetrisch zur reellen Achse.
Um die Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 grafisch darzustellen, müssen die Werte des Parameters 
  
    
      
        c
      
    
    {\displaystyle c}
  
 alle einzeln bis zu einer selbstbestimmten Anzahl von Iterationen berechnet werden.
Für 
  
    
      
        c
        =
        −
        2
         
      
    
    {\displaystyle c=-2\ }
  
 lautet die Folge 
  
    
      
        0
        ,
        
        −
        2
        ,
         
        2
        ,
         
        2
        ,
         
        2
        ,
         
        2
        ,
        …
      
    
    {\displaystyle 0,\,-2,\ 2,\ 2,\ 2,\ 2,\dotsc }
  
 und ist beschränkt. Daher ist 
  
    
      
        c
        =
        −
        2
         
      
    
    {\displaystyle c=-2\ }
  
 Element von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
.
Für 
  
    
      
        c
        =
        +
        2
         
      
    
    {\displaystyle c=+2\ }
  
 lautet die Folge 
  
    
      
        0
        ,
         
        2
        ,
         
        6
        ,
         
        38
        ,
         
        1446
        ,
        …
      
    
    {\displaystyle 0,\ 2,\ 6,\ 38,\ 1446,\dotsc }
  
 und ist divergent. Daher ist 
  
    
      
        c
        =
        +
        2
         
      
    
    {\displaystyle c=+2\ }
  
 kein Element von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
.


### Definition über Julia-Mengen
Die Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 wurde von Benoît Mandelbrot ursprünglich zur Klassifizierung von Julia-Mengen eingeführt, die bereits Anfang des 20. Jahrhunderts von den französischen Mathematikern Gaston Maurice Julia und Pierre Fatou untersucht wurden. Die Julia-Menge 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 zu einer bestimmten komplexen Zahl 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ist definiert als der Rand der Menge aller Anfangswerte 
  
    
      
        
          z
          
            0
          
        
      
    
    {\displaystyle z_{0}}
  
, für die die obige Zahlenfolge beschränkt bleibt. Es kann bewiesen werden, dass die Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 genau die Menge der Werte 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ist, für die die zugehörige Julia-Menge 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 zusammenhängend ist.
Dieses Prinzip wird in vielen Resultaten über das Verhalten der Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 vertieft. So zeigt Shishikura, dass der Rand der Mandelbrot-Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 ebenso wie die zugehörige Julia-Menge 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 die Hausdorff-Dimension 2 hat. Ein unveröffentlichtes Manuskript von Jean-Christophe Yoccoz diente John Hamal Hubbard als Grundlage für seine Ergebnisse über lokal zusammenhängende Julia-Mengen 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 und lokal zusammenhängende Mandelbrot-Mengen 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
.


### Bezug zur Chaostheorie
Das Bildungsgesetz, das der Folge zugrunde liegt, ist die einfachste nichtlineare Gleichung, anhand deren sich der Übergang von Ordnung zu Chaos durch Variation eines Parameters provozieren lässt. Dazu genügt es, reelle Zahlenfolgen zu betrachten.
Für Werte 
  
    
      
        −
        0
        
          ,
        
        75
        ≤
        c
        ≤
        0
        
          ,
        
        25
      
    
    {\displaystyle -0{,}75\leq c\leq 0{,}25}
  
, das heißt innerhalb der Kardioide, konvergiert die Folge. Auf der „Antenne“, die bis 
  
    
      
        c
        =
        −
        2
      
    
    {\displaystyle c=-2}
  
 reicht, verhält sich die Folge chaotisch. Der Übergang zu chaotischem Verhalten erfolgt nun über ein Zwischenstadium mit periodischen Grenzzyklen. Dabei nimmt die Periode zum chaotischen Bereich hin stufenweise um den Faktor zwei zu, ein Phänomen, das als Periodenverdopplung und Bifurkation bezeichnet wird. Jeder 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Bereich zu einer bestimmten Periode entspricht dabei einer der kreisförmigen „Knospen“ auf der 
  
    
      
        x
      
    
    {\displaystyle x}
  
-Achse.
Die Periodenverdopplung beginnt mit dem „Kopf“ und setzt sich in der Folge der „Knospen“ zur „Antenne“ hin fort. Das Verhältnis der Längen aufeinander folgender Parameterintervalle und damit das der Knospendurchmesser zu unterschiedlicher Periode strebt dabei gegen die Feigenbaum-Konstante 
  
    
      
        δ
        ≈
        4,669
      
    
    {\displaystyle \delta \approx 4{,}669}
  
, eine fundamentale Konstante der Chaostheorie. Dieses Verhalten ist typisch für den Übergang realer Systeme zu chaotischer Dynamik. Die auffälligen Lücken im chaotischen Bereich entsprechen Inseln mit periodischem Verhalten, denen in der komplexen Ebene die Satelliten auf der „Antenne“ zugeordnet sind.
Für gewisse komplexe 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte stellen sich Grenzzyklen ein, die auf einer geschlossenen Kurve liegen, deren Punkte jedoch nicht periodisch, sondern chaotisch abgedeckt werden. Eine solche Kurve ist in der Chaostheorie als sogenannter seltsamer Attraktor bekannt.
Die Mandelbrot-Menge ist daher ein elementares Objekt für die Chaostheorie, an dem sich fundamentale Phänomene studieren lassen. Sie wird aus diesem Grund hinsichtlich ihrer Bedeutung für die Chaostheorie gelegentlich mit der von Geraden für die euklidische Geometrie verglichen.


### Geometrische und mathematische Eigenschaften
Die Mandelbrot-Menge ist abgeschlossen (da ihr Komplement offen ist) und in der abgeschlossenen Scheibe mit dem Radius 2 um den Ursprung enthalten und somit kompakt.
Es sei 
  
    
      
        
          P
          
            c
          
        
        (
        z
        )
        =
        
          z
          
            2
          
        
        +
        c
      
    
    {\displaystyle P_{c}(z)=z^{2}+c}
  
 und 
  
    
      
        
          P
          
            c
          
          
            n
          
        
        (
        z
        )
      
    
    {\displaystyle P_{c}^{n}(z)}
  
 bezeichne die 
  
    
      
        n
      
    
    {\displaystyle n}
  
-te Iteration. Ein Punkt 
  
    
      
        c
      
    
    {\displaystyle c}
  
 gehört genau dann zur Mandelbrot-Menge, falls 
  
    
      
        
          |
        
        
          P
          
            c
          
          
            n
          
        
        (
        0
        )
        
          |
        
        ≤
        2
      
    
    {\displaystyle |P_{c}^{n}(0)|\leq 2}
  
 für alle 
  
    
      
        n
        ≥
        0
      
    
    {\displaystyle n\geq 0}
  
.
Wird der Betrag von 
  
    
      
        
          P
          
            c
          
          
            n
          
        
        (
        0
        )
      
    
    {\displaystyle P_{c}^{n}(0)}
  
 größer als 2, dann entkommt die Iteration ins Unendliche, der Betrag wächst über jede Grenze. 
  
    
      
        c
      
    
    {\displaystyle c}
  
 gehört dann nicht zur Mandelbrot-Menge.
Der ungeheure Formenreichtum der Mandelbrot-Menge erschließt sich aus ihrem Bezug zu Julia-Mengen. Julia-Mengen zur Iteration 
  
    
      
        z
        →
        
          z
          
            2
          
        
        +
        c
      
    
    {\displaystyle z\to z^{2}+c}
  
 sind Fraktale, außer für einige 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte wie 
  
    
      
        c
        =
        −
        2
      
    
    {\displaystyle c=-2}
  
 (Strecke) oder 
  
    
      
        c
        =
        0
      
    
    {\displaystyle c=0}
  
 (Kreis). Die Formen dieser fraktalen Strukturen sind innerhalb einer Julia-Menge stets die gleichen, umspannen aber für Julia-Mengen zu verschiedenem Parameter 
  
    
      
        c
      
    
    {\displaystyle c}
  
 einen enormen Formenreichtum. Es zeigt sich, dass die Strukturen der Mandelbrot-Menge in der Umgebung eines bestimmten Wertes 
  
    
      
        c
      
    
    {\displaystyle c}
  
 genau die Strukturen der zugehörigen Julia-Menge 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 wiedergeben. Damit enthält die Mandelbrot-Menge den kompletten Formenreichtum der unendlich vielen Julia-Mengen (s. u.).
In den fraktalen Strukturen am Rand finden sich verkleinerte ungefähre Kopien der gesamten Mandelbrot-Menge, die Satelliten. Jeder Bildausschnitt der Mandelbrot-Menge, der sowohl Punkte aus 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 als auch solche außerhalb von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 umfasst, enthält unendlich viele dieser Satelliten. Unmittelbar am Rand eines Satelliten treten fast die gleichen Strukturen auf wie an den entsprechenden Stellen des Originals. Diese Strukturen sind jedoch nach weiter außen hin mit den Strukturen kombiniert, die für die größere Umgebung des Satelliten typisch sind.
Da jeder Satellit wiederum mit Satelliten höherer Ordnung bestückt ist, lässt sich immer eine Stelle finden, an der eine beliebige Anzahl beliebiger verschiedener Strukturen in beliebiger Reihenfolge kombiniert auftritt. Diese Strukturen sind allerdings nur bei extremer Vergrößerung erkennbar.
Die Mandelbrot-Menge ist spiegelsymmetrisch zur reellen Achse. Sie ist zusammenhängend (das heißt, sie bildet keine Inseln), wie Adrien Douady und John Hamal Hubbard 1984 bewiesen, und es wird vermutet (Douady/Hubbard), dass sie lokal zusammenhängend ist (MLC-Vermutung). Dies ist eine der großen offenen Fragen in der komplexen Dynamik und bisher unbewiesen (obwohl es Teilresultate zum Beispiel von Jean-Christophe Yoccoz gibt, der lokalen Zusammenhang für bestimmte Werte von 
  
    
      
        c
      
    
    {\displaystyle c}
  
 bewies, für die endlich-renormalisierbaren Punkte). Die MLC erlaubt weitreichende Folgerungen über die Topologie der Mandelbrot-Menge. Beispielsweise würde daraus die Hyperbolizitätsvermutung folgen, dass jede offene Menge in der Mandelbrot-Menge (also das Innere der Mandelbrot-Menge) aus Punkten mit attraktiven Zyklen besteht. Die Mandelbrot-Menge ist zwar selbstähnlich, aber nicht exakt, denn keine zwei Teilstrukturen ihres Randes sind exakt gleich; aber in der Nähe vieler Randpunkte bilden sich bei fortgesetzter Ausschnittvergrößerung im Grenzfall periodische Strukturen. An speziellen Punkten hat die Mandelbrot-Menge Selbstähnlichkeit (vermutet von John Milnor und bewiesen von Mikhail Lyubich 1999).
Da die Mandelbrot-Menge Kardioid- und Kreisflächen enthält, hat sie die fraktale Dimension 2. Der Rand der Mandelbrot-Menge hat eine unendliche Länge, und seine Hausdorff-Dimension beträgt nach Arbeiten von Mitsuhiro Shishikura ebenfalls 2; das impliziert, dass die Box-Dimension den Wert 2 hat. Es ist denkbar, dass der Rand der Mandelbrot-Menge einen positiven (notwendig endlichen) Flächeninhalt hat; andernfalls wäre dieser Flächeninhalt null. Der Flächeninhalt der Mandelbrot-Menge ist nicht bekannt und beträgt nach numerischen Schätzungen etwa 1,5065918849.
Die Mandelbrot-Menge enthält deformierte Kopien aller Julia-Mengen, wie Tan Lei 1990 für die Misiurewicz-Punkte der Mandelbrot-Menge bewiesen hat, die dicht im Rand der Mandelbrot-Menge liegen. Das ist ein weiterer Beleg für die enge Verwandtschaft der Struktur von Julia- und Mandelbrot-Mengen. So wurden in den Beweisen von Yoccoz für lokalen Zusammenhang der Mandelbrot-Menge bei endlich renormalisierbaren Punkten und von Shishikura über die fraktale Dimension des Randes der Mandelbrot-Menge zuerst die entsprechenden Eigenschaften bei den zum Parameterwert gehörigen Julia-Mengen untersucht und dann auf die Mandelbrot-Menge übertragen.
Die Frage, ob die Mandelbrot-Menge entscheidbar ist, ergibt zunächst keinen Sinn, da 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 überabzählbar ist. Einen Ansatz, den Begriff der Entscheidbarkeit auf überabzählbare Mengen zu verallgemeinern, stellt das Blum-Shub-Smale-Modell dar. Innerhalb dessen ist die Mandelbrot-Menge nicht entscheidbar.


## Bildergalerien
### Galerie der Iterationen
Die folgende Galerie gibt einen Überblick über die Nullstellen einiger Polynome 
  
    
      
        
          z
          
            n
          
        
        (
        c
        )
      
    
    {\displaystyle z_{n}(c)}
  
 der Iteration 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  
 mit 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
. Der Realteil 
  
    
      
        x
      
    
    {\displaystyle x}
  
 von 
  
    
      
        c
        =
        x
        +
        i
        y
      
    
    {\displaystyle c=x+iy}
  
 erstreckt sich in den Bildern von −2,2 bis 1, der Imaginärteil 
  
    
      
        y
      
    
    {\displaystyle y}
  
 von −1,4 bis 1,4.


### Galerie einer Zoomfahrt
Die folgende Galerie zeigt Bilder einer Zoomfahrt in die Mandelbrot-Menge. Jedes Bild ist dabei die Vergrößerung eines kleinen Ausschnitts des vorangegangenen Bildes. Dazu wird in den 15 Bildern jeweils um einen Faktor 4 bis 8 immer tiefer in die Mandelbrot-Menge hineingezoomt. Das letzte Bild erreicht auf diese Weise einen Vergrößerungsfaktor von etwa 60 Milliarden. Die Sequenz gibt einen Eindruck vom geometrischen Formenreichtum und erläutert gewisse typische Strukturelemente.


### Galerie ausgewählter Kunstwerke
## Verhalten der Zahlenfolge
Die verschiedenen Strukturelemente von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 stehen in engem Zusammenhang mit bestimmten Verhaltensweisen der Zahlenfolge, die 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 zugrunde liegt. Je nach Wert von 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ergibt sich eine der folgenden vier Möglichkeiten:

_    Sie konvergiert gegen einen Fixpunkt.
____ Sie konvergiert gegen einen periodischen Grenzzyklus, der aus zwei oder mehr Werten besteht. Dazu zählen auch die Fälle, in denen sich die Folge von Anfang periodisch verhält.
_    Sie wiederholt sich nie, bleibt aber beschränkt. Manche Werte zeigen chaotisches Verhalten mit Wechsel zwischen fast periodischen Grenzyklen und scheinbar zufälligem Verhalten.
_    Sie divergiert gegen Unendlich (bestimmte Divergenz).
Alle 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte, die nicht bestimmt divergieren, gehören zu 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
.
Die folgende Tabelle zeigt Beispiele für diese vier Grenzverhalten der Iteration 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  
 für 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
:


### Geometrische Zuordnung
Konvergenz liegt genau für die Werte von 
  
    
      
        c
      
    
    {\displaystyle c}
  
 vor, die das Innere der Kardioide bilden, den „Körper“ von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
, sowie für abzählbar viele ihrer Randpunkte. Periodische Grenzzyklen finden sich in den (angenähert) kreisförmigen „Knospen“ wie im „Kopf“, in den Kardioiden der Satelliten sowie ebenfalls auf abzählbar vielen Randpunkten dieser Komponenten. Eine fundamentale Vermutung besagt, dass es für alle inneren Punkte der Mandelbrot-Menge einen Grenzzyklus gibt. Die Folge ist echt vorperiodisch für abzählbar viele Parameter, die oft Misiurewicz-Thurston-Punkte genannt werden (nach Michał Misiurewicz und William Thurston). Dazu gehören die „Antennenspitzen“ wie der Punkt 
  
    
      
        z
        =
        −
        2
      
    
    {\displaystyle z=-2}
  
 ganz links und Verzweigungspunkte der Mandelbrot-Menge.
In den überabzählbar vielen übrigen Punkten der Mandelbrot-Menge kann sich die Folge auf viele verschiedene Weisen verhalten, die jeweils sehr unterschiedliche dynamische Systeme erzeugen und die teilweise Gegenstand intensiver Forschung sind. Je nach Definition des Wortes lässt sich „chaotisches“ Verhalten finden.


### Periodisches Verhalten
#### Die kreisförmigen Strukturen
Jede kreisförmige „Knospe“ und jede Satelliten-Kardioide zeichnet sich durch eine bestimmte Periodizität des Grenzzyklus aus, gegen den die Folge für die zugehörigen 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte strebt. Die Anordnung der „Knospen“ an der zugehörigen Kardioide folgt dabei den folgenden Regeln, aus denen sich unmittelbar die Periodizitäten ablesen lassen. Jede „Knospe“ berührt genau einen Basiskörper, nämlich eine größere „Knospe“ oder eine Kardioide.
Die Periodizität einer „Knospe“ ist die Summe der Periodizitäten der beiden nächsten größeren „Nachbarknospen“ in beide Richtungen am selben Basiskörper, sofern es solche gibt. Gibt es am Rand des Basiskörpers bis zur Kontaktstelle mit dessen Basiskörper oder bis zur Kerbe der Kardioide nur kleinere „Knospen“, so trägt anstelle der Periodizität einer „Nachbarknospe“ die des Basiskörpers selbst zur Summe bei. Daraus leiten sich unmittelbar die folgenden Eigenschaften ab:

Tendenziell sind die „Knospen“ oder Kardioiden umso kleiner, je größer ihre Periodizität ist.
Die Periodizität der größten „Knospe“ an einem Basiskörper beträgt stets das Doppelte, wie der „Dutt“ mit der Periode 4 am „Kopf“.
Die Periodizität einer „Knospe“ eines Satelliten ist das Produkt der Periodizität der Satelliten-Kardioide und der der korrespondierenden „Knospe“ der Hauptkardioide.
Ferner erklärt diese Regel das Auftreten bestimmter Folgen von „Knospen“ wie vom „Kopf“ zur Kardioidkerbe hin mit einer Periodizitätszunahme zur nächsten „Knospe“ hin um den Wert 1 oder vom „Arm“ zum „Kopf“ hin um den Wert 2.


#### Attraktive Zyklen
Gibt es für ein 
  
    
      
        c
      
    
    {\displaystyle c}
  
 ein Folgenglied mit der Eigenschaft 
  
    
      
        
          z
          
            n
          
        
        =
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{n}=z_{0}=0}
  
, so wiederholt sich die Folge von Anfang an streng periodisch und zwar mit der Periode 
  
    
      
        n
      
    
    {\displaystyle n}
  
. Da sich 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 durch 
  
    
      
        n
      
    
    {\displaystyle n}
  
-malige Anwendung der Iterationsvorschrift ergibt, wobei bei jedem Schritt quadriert wird, lässt es sich als Polynom von 
  
    
      
        c
      
    
    {\displaystyle c}
  
 vom Grad 
  
    
      
        
          2
          
            n
            −
            1
          
        
      
    
    {\displaystyle 2^{n-1}}
  
 formulieren. Die 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte für periodische Folgen der Periode 
  
    
      
        n
      
    
    {\displaystyle n}
  
 werden daher über die 
  
    
      
        
          2
          
            n
            −
            1
          
        
      
    
    {\displaystyle 2^{n-1}}
  
 Nullstellen dieses Polynoms erhalten. Es zeigt sich, dass jede Zahlenfolge gegen diesen Zahlenzyklus konvergiert, sofern eins ihrer Folgenglieder hinreichend nahe an diesen Zyklus gerät, die werden Attraktoren genannt. Das führt dazu, dass alle Zahlenfolgen zu einer gewissen Umgebung des 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Wertes, der den Attraktor repräsentiert, gegen einen stabilen Zyklus der Periode 
  
    
      
        n
      
    
    {\displaystyle n}
  
 konvergieren. Jede kreisförmige „Knospe“ und jede Kardioide eines Satelliten repräsentiert genau eine solche Umgebung. Exemplarisch seien die Gebiete mit den Perioden 1 bis 3 aufgeführt:

Periode 1: Die Kardioide des Hauptapfelmännchens. Der Rand dieser Kardioide ist gegeben durch Punkte der Form 
  
    
      
        c
        =
        
          
            
              1
              2
            
          
        
        z
        −
        
          
            
              1
              4
            
          
        
        
          z
          
            2
          
        
      
    
    {\displaystyle c={\tfrac {1}{2}}z-{\tfrac {1}{4}}z^{2}}
  
 mit 
  
    
      
        
          |
        
        z
        
          |
        
        =
        1
      
    
    {\displaystyle |z|=1}
  
.
Periode 2: Der „Kopf“. Die 2. Nullstelle 
  
    
      
        c
        =
        0
      
    
    {\displaystyle c=0}
  
 entspricht der Hauptkardioide, die wegen der Periode 1 natürlich bei der Ermittlung aller höherer Perioden als Nullstelle auftritt. Diese Überlegung zeigt, dass die Zahl der Attraktoren mit der Periode 
  
    
      
        n
        >
        1
      
    
    {\displaystyle n>1}
  
 maximal 
  
    
      
        
          2
          
            n
            −
            1
          
        
        −
        1
      
    
    {\displaystyle 2^{n-1}-1}
  
 betragen kann, und das nur dann, wenn 
  
    
      
        n
      
    
    {\displaystyle n}
  
 eine Primzahl ist. Der Kopf selbst ist eine Kreisscheibe mit Mittelpunkt 
  
    
      
        −
        1
      
    
    {\displaystyle -1}
  
 und Radius 
  
    
      
        
          
            
              1
              4
            
          
        
      
    
    {\displaystyle {\tfrac {1}{4}}}
  
, d. h., der Rand dieser Kreisscheibe ist gegeben durch Punkte der Form 
  
    
      
        c
        =
        −
        1
        +
        
          
            
              1
              4
            
          
        
        z
      
    
    {\displaystyle c=-1+{\tfrac {1}{4}}z}
  
 mit 
  
    
      
        
          |
        
        z
        
          |
        
        =
        1
      
    
    {\displaystyle |z|=1}
  
.
Periode 3: Die „Knospen“, die den „Armen“ entsprechen und die Kardioide des größten Satelliten auf der „Kopfantenne“. Die vierte Nullstelle 
  
    
      
        c
        =
        0
      
    
    {\displaystyle c=0}
  
 entfällt wieder.
Die Anzahl der anziehenden Zyklen mit der genauen Periode 
  
    
      
        n
      
    
    {\displaystyle n}
  
, d. h. 
  
    
      
        
          z
          
            n
          
        
        =
        
          z
          
            0
          
        
      
    
    {\displaystyle z_{n}=z_{0}}
  
 und 
  
    
      
        n
      
    
    {\displaystyle n}
  
 ist minimal mit dieser Eigenschaft, ist die Folge A000740 in OEIS.


#### Repulsive Zyklen
Neben attraktiven Zyklen gibt es repulsive, die sich dadurch auszeichnen, dass Zahlenfolgen in ihrer Umgebung sich zunehmend von ihnen entfernen. Sie lassen sich jedoch erreichen, da jedes 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 abgesehen von der Situation 
  
    
      
        
          z
          
            n
            −
            1
          
        
        =
        0
      
    
    {\displaystyle z_{n-1}=0}
  
 wegen des Quadrats in der Iterationsvorschrift zwei potenzielle Vorgänger in der Folge hat, die sich nur durch ihr Vorzeichen unterscheiden. 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte, für die die zugehörige Folge irgendwann über einen solchen zweiten Vorläufer eines Periodenmitgliedes in einen derartigen instabilen Zyklus mündet, sind beispielsweise die „Naben“ der rad- oder spiralförmigen Strukturen sowie die Endpunkte der weitverbreiteten antennenartigen Strukturen, die sich formal als „Naben“ von „Rädern“ oder Spiralen mit einer einzigen Speiche interpretieren lassen. Derartige 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte werden als Misiurewicz-Punkte bezeichnet.
Ein Misiurewicz-Punkt 
  
    
      
        c
      
    
    {\displaystyle c}
  
 hat ferner die Eigenschaft, dass 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 in seiner näheren Umgebung nahezu deckungsgleich mit demselben Ausschnitt der zugehörige Julia-Menge 
  
    
      
        
          J
          
            c
          
        
      
    
    {\displaystyle J_{c}}
  
 ist. Je weiter sich dem Misiurewicz-Punkt genähert wird, umso besser wird die Übereinstimmung. Da Julia-Mengen für 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte innerhalb von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 zusammenhängend sind und außerhalb von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 Cantor-Mengen aus unendlich vielen Inseln mit der Gesamtfläche null, sind sie in der Übergangszone am Rand von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 besonders filigran. Jeder Misiurewicz-Punkt ist aber gerade ein Randpunkt von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
, und jeder Ausschnitt der Randzone von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
, der sowohl Punkte in 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 als auch außerhalb davon enthält, enthält unendlich viele davon. Damit ist der gesamte Formenreichtum sämtlicher Julia-Mengen dieses filigranen Typs in der Umgebung der Misiurewicz-Punkte in 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 repräsentiert.


#### Satelliten
Ein weiteres Strukturelement, das den Formenreichtum der Mandelbrot-Menge begründet, sind die verkleinerten Kopien ihrer selbst, die sich in den filigranen Strukturen ihres Randes befinden. Dabei korrespondiert das Verhalten der Zahlenfolgen innerhalb eines Satelliten in folgender Weise mit dem der Folgen im Hauptkörper. Innerhalb eines Satelliten konvergieren alle Zahlenfolgen gegen Grenzzyklen, deren Perioden sich von denen an den entsprechenden Stellen im Hauptkörper von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 um einen Faktor 
  
    
      
        p
      
    
    {\displaystyle p}
  
 unterscheiden. Wird für einen bestimmten 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Wert aus dem Satelliten nur jedes 
  
    
      
        p
      
    
    {\displaystyle p}
  
-te Folgenglied betrachtet, so ergibt sich eine Folge, die bis auf einen räumlichen Maßstabsfaktor nahezu identisch ist mit derjenigen, die sich für den entsprechenden 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Wert im Hauptkörper von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 ergibt. Die mathematische Begründung hierfür ist tiefliegend; sie entstammt den Arbeiten von Douady und Hubbard über „polynomartige Abbildungen“.
Die zusätzlichen Strukturelemente in der unmittelbaren Umgebung eines Satelliten sind eine Folge davon, dass zwischen zweien der betrachteten Folgenglieder mit dem Indexabstand 
  
    
      
        p
      
    
    {\displaystyle p}
  
 sich eines mit dem Wert 
  
    
      
        
          z
          
            n
          
        
        =
        0
      
    
    {\displaystyle z_{n}=0}
  
 befinden kann, das damit einen periodischen Verlauf mit der Periode 
  
    
      
        n
      
    
    {\displaystyle n}
  
 begründet. Die entsprechende Folge außerhalb des Hauptkörpers divergiert jedoch, da sie keine solchen Zwischenglieder besitzt.
Es handelt sich bei der Mandelbrot-Menge selbst um eine universelle Struktur, die bei völlig anderen nichtlinearen Systemen und Klassifizierungsregeln in Erscheinung treten kann. Grundvoraussetzung ist jedoch, dass die beteiligten Funktionen winkeltreu sind. Werden solche Systeme betrachtet, die von einem komplexen Parameter 
  
    
      
        c
      
    
    {\displaystyle c}
  
 abhängen, und klassifiziert man ihr Verhalten bezüglich einer bestimmten Eigenschaft der Dynamik in Abhängigkeit von 
  
    
      
        c
      
    
    {\displaystyle c}
  
, dann werden unter bestimmten Umständen in der Parameter-Ebene kleine Kopien der Mandelbrot-Menge gefunden. Ein Beispiel ist die Frage, für welche Polynome dritten Grades das iterative Newton-Verfahren zur Bestimmung von Nullstellen mit einem bestimmten Startwert versagt und für welche nicht.
Wie im Bild kann die Mandelbrot-Menge dabei verzerrt auftreten, zum Beispiel sitzen dort die Armknospen an etwas anderer Stelle. Ansonsten ist die Mandelbrot-Menge jedoch vollkommen intakt, inklusive aller Knospen, Satelliten, Filamente und Antennen. Der Grund für das Auftauchen der Mandelbrot-Menge ist, dass die betrachtete Funktionenfamilien in bestimmten Gebieten – abgesehen von Drehungen und Verschiebungen – recht gut mit der Funktionenfamilie 
  
    
      
        {
        
          f
          
            c
          
        
        :
        z
        ↦
        
          z
          
            2
          
        
        +
        c
        ∣
        c
        ∈
        
          C
        
        }
      
    
    {\displaystyle \{f_{c}\colon z\mapsto z^{2}+c\mid c\in \mathbb {C} \}}
  
, welche die Mandelbrot-Menge definiert, übereinstimmen. Dabei sind in einem gewissen Rahmen Abweichungen zulässig, und trotzdem kristallisiert sich die Mandelbrot-Menge heraus. Dieses Phänomen wird als strukturelle Stabilität bezeichnet und ist im Endeffekt verantwortlich für das Auftreten der Satelliten in der Umgebung von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
, weil Teilfolgen der iterierten Funktionen lokal das gleiche Verhalten zeigen wie die Gesamtfamilie.


#### Intermediär wechselhaftes Verhalten
Durch die Möglichkeit der Zahlenfolge, wiederholt in die unmittelbare Umgebung eines repulsiven Zyklus zu geraten, und bei dem anschließend tendenziell divergenten oder chaotischen Verhalten wiederum beinahe in einen anderen Zyklus zu geraten, können sich intermediär sehr komplizierte Verhaltensweisen der Folge ausbilden, bis sich der endgültige Charakter der Folge zeigt, wie die beiden Abbildungen demonstrieren. Die Umgebung der zugehörigen 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Werte in 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 ist entsprechend strukturreich.
Die Darstellung der Folgepunkte selbst in der komplexen Ebene zeigt in diesen Fällen eine größere Komplexität. Das quasiperiodische Verhalten in der Nachbarschaft eines repulsiven Zyklus führt in diesen Fällen oft zu spiralförmigen Strukturen mit mehreren Armen, wobei die Folgepunkte das Zentrum umkreisen, während der Abstand zu ihm zunimmt. Die Anzahl der Arme entspricht daher der Periode. Die Punktanhäufungen an den Enden der Spiralarme in der zweiten Abbildung sind die Folge der beiden zugehörigen Beinahe-Einfänge durch repulsive (instabile) Zyklen.


### Dichteverteilung der Folgenglieder
Das Bild zeigt in der komplexen Ebene die Dichteverteilung der Folgenglieder, die sich durch Auswertung von 60 Millionen Folgen ergibt, wobei die Helligkeit ein Maß dafür ist, wie viele Orbitale durch den Punkt verlaufen. Blaue Bereiche kennzeichnen Folgenglieder mit kleinem Index, während eine gelbliche Färbung Folgenglieder mit hohen Indizes anzeigt. Folgen aus der großen Kardioide der Mandelbrot-Menge tendieren zu einer Konvergenz zu einem 
  
    
      
        c
      
    
    {\displaystyle c}
  
-Wert auf einem Kreis um den Ursprung, der als runder Bereich mit sehr hoher Dichte zu erkennen ist. Die kleineren Gebilde nahe der imaginären Achse markieren die konjugierten Bereiche, zwischen denen Folgenglieder vieler Folgen hin und her springen (Annäherung an einen Grenzzyklus mit der Periode 2 für große Folgenindizes).


## Praktische Berechnung
Die grafische Darstellung der Mandelbrot-Menge und ihrer Strukturen im Randbereich ist nur mittels Computer durch sogenannte Fraktalgeneratoren möglich. Dabei entspricht jeder Bildpunkt einer Zahl 
  
    
      
        c
      
    
    {\displaystyle c}
  
 der komplexen Ebene. Der Computer ermittelt für alle Bildpunkte, ob die zugehörigen Folgen divergieren oder nicht. Sobald der Betrag 
  
    
      
        
          |
        
        
          z
          
            n
          
        
        
          |
        
      
    
    {\displaystyle |z_{n}|}
  
 eines Folgengliedes die Fluchtgrenze 
  
    
      
        s
        =
        2
      
    
    {\displaystyle s=2}
  
 überschreitet, ist klar, dass die Folge divergiert; die Iterationszahl 
  
    
      
        n
      
    
    {\displaystyle n}
  
 kann als Maß für den Divergenzgrad dienen. Die Bildpunkte werden nach einer Farbtabelle eingefärbt, die jeder Iterationszahl eine Farbe zuordnet.
Für harmonische Farbübergänge wird in der Praxis ein Fluchtradius 
  
    
      
        r
        >
        2
      
    
    {\displaystyle r>2}
  
 gewählt, da andernfalls die Breite der Farbstreifen zu sehr oszilliert. Je größer der Fluchtradius gewählt wird, desto besser entsprechen die Farbgrenzen den Äquipotentiallinien, die sich ergeben, wenn die Mandelbrot-Menge als elektrisch geladener Leiter interpretiert wird. Für kontinuierliche Farbverläufe, wie in der obigen Galerie einer Zoomfahrt, ist es erforderlich, entweder die genaue Fluchtzeit 
  
    
      
        t
      
    
    {\displaystyle t}
  
 zu ermitteln, bei der der Fluchtradius 
  
    
      
        r
      
    
    {\displaystyle r}
  
 gerade erreicht wird, oder den Abstand 
  
    
      
        d
      
    
    {\displaystyle d}
  
 zum Rand der Menge abzuschätzen.
Da die Anzahl der Iterationsschritte 
  
    
      
        n
      
    
    {\displaystyle n}
  
, nach der die Fluchtgrenze 
  
    
      
        s
        =
        2
      
    
    {\displaystyle s=2}
  
 überschritten wird, beliebig groß sein kann, muss als Abbruchkriterium eine maximale Iterationszahl 
  
    
      
        
          n
          
            max
          
        
      
    
    {\displaystyle n_{\text{max}}}
  
 festgelegt werden. Alle Zahlen 
  
    
      
        c
      
    
    {\displaystyle c}
  
, deren Folgen danach einen Fluchtradius 
  
    
      
        r
        >
        2
      
    
    {\displaystyle r>2}
  
 noch nicht erreicht haben, werden zur Menge 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 gezählt. Je geringer der Abstand von 
  
    
      
        c
      
    
    {\displaystyle c}
  
 zum Rand der Menge ist, desto größer ist in der Regel die Iterationszahl 
  
    
      
        n
      
    
    {\displaystyle n}
  
, bei der 
  
    
      
        r
      
    
    {\displaystyle r}
  
 erreicht oder überschritten wird. Je stärker die Vergrößerung ist, mit der der Rand dargestellt wird, desto größer muss die maximale Iterationszahl 
  
    
      
        
          n
          
            max
          
        
      
    
    {\displaystyle n_{\text{max}}}
  
 gewählt werden, und umso mehr Rechenzeit ist nötig.
Ein erheblicher Teil der Rechenzeit wird dort benötigt, wo die Zahlenfolge nicht divergiert. Moderne Programme versuchen daher, die Rechenzeit an diesen Stellen durch verschiedene Verfahren zu verkürzen. Eine Möglichkeit besteht darin, die Berechnung abzubrechen, sobald die Zahlenfolge konvergiert oder in einem periodischen Zyklus gefangen ist. Andere Programme nutzen aus, dass das Innere einer geschlossenen Kurve automatisch in 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 liegt, wenn der Rand bereits zu 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 gehört.
Grafisch besonders reizvoll ist die Darstellung des Randes von 
  
    
      
        
          M
        
      
    
    {\displaystyle \mathbb {M} }
  
 mit seinem Formenreichtum. Je stärker die gewählte Vergrößerung ist, umso komplexere Strukturen lassen sich dort finden. Mit entsprechenden Computerprogrammen lässt sich dieser Rand wie mit einem Mikroskop mit beliebiger Vergrößerung darstellen. Zur Untersuchung interessanter Strukturen sind häufig Vergrößerungen erforderlich, die mit hardwareunterstützten Datentypen aufgrund deren limitierter Genauigkeit nicht berechnet werden können. Manche Programme enthalten daher Langzahlarithmetik-Datentypen mit beliebig wählbarer Genauigkeit. Damit sind (fast) beliebige Vergrößerungsfaktoren möglich.


### Interpolation und Distanzschätzung
Wenn der Fluchtradius 
  
    
      
        r
      
    
    {\displaystyle r}
  
 groß genug ist, kann zur Ermittlung der Fluchtzeit 
  
    
      
        t
      
    
    {\displaystyle t}
  
 anstelle der eigentlichen Folge 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  
 mit 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
 die vereinfachte Folge 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        0
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+0}
  
 mit 
  
    
      
        
          z
          
            0
          
        
        =
        c
      
    
    {\displaystyle z_{0}=c}
  
 betrachtet werden. Die ersten Folgenglieder sind dann 
  
    
      
        
          z
          
            0
          
        
        =
        c
      
    
    {\displaystyle z_{0}=c}
  
, 
  
    
      
        
          z
          
            1
          
        
        =
        
          c
          
            2
          
        
      
    
    {\displaystyle z_{1}=c^{2}}
  
, 
  
    
      
        
          z
          
            2
          
        
        =
        
          c
          
            4
          
        
      
    
    {\displaystyle z_{2}=c^{4}}
  
 und 
  
    
      
        
          z
          
            3
          
        
        =
        
          c
          
            8
          
        
      
    
    {\displaystyle z_{3}=c^{8}}
  
, ganz allgemein ergibt sich 
  
    
      
        
          z
          
            n
          
        
        =
        
          c
          
            
              2
              
                n
              
            
          
        
      
    
    {\displaystyle z_{n}=c^{2^{n}}}
  
. Betrachtet man eine Iterationszahl 
  
    
      
        n
      
    
    {\displaystyle n}
  
, bei der die Folge einen Fluchtradius 
  
    
      
        r
        ≫
        2
      
    
    {\displaystyle r\gg 2}
  
 erreicht oder überschreitet, so gilt 
  
    
      
        
          |
        
        
          z
          
            n
          
        
        
          |
        
        =
        
          |
        
        c
        
          
            |
          
          
            
              2
              
                n
              
            
          
        
        ≥
        r
      
    
    {\displaystyle |z_{n}|=|c|^{2^{n}}\geq r}
  
 und 
  
    
      
        n
        =
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                
                  |
                
              
              
                ln
                ⁡
                
                  |
                
                c
                
                  |
                
              
            
          
        
        )
      
    
    {\displaystyle n=\log _{2}({\tfrac {\ln |z_{n}|}{\ln |c|}})}
  
. Bestimmt man nun die Fluchtzeit 
  
    
      
        t
        ≤
        n
      
    
    {\displaystyle t\leq n}
  
 aus der Interpolation 
  
    
      
        
          |
        
        c
        
          
            |
          
          
            
              2
              
                t
              
            
          
        
        =
        r
      
    
    {\displaystyle |c|^{2^{t}}=r}
  
, so ergibt sich mit 
  
    
      
        t
        =
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                r
              
              
                ln
                ⁡
                
                  |
                
                c
                
                  |
                
              
            
          
        
        )
      
    
    {\displaystyle t=\log _{2}({\tfrac {\ln r}{\ln |c|}})}
  
:

  
    
      
        n
        −
        t
        =
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                
                  |
                
              
              
                ln
                ⁡
                
                  |
                
                c
                
                  |
                
              
            
          
        
        )
        −
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                r
              
              
                ln
                ⁡
                
                  |
                
                c
                
                  |
                
              
            
          
        
        )
        =
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                
                  |
                
              
              
                ln
                ⁡
                r
              
            
          
        
        )
        .
      
    
    {\displaystyle n-t=\log _{2}({\tfrac {\ln |z_{n}|}{\ln |c|}})-\log _{2}({\tfrac {\ln r}{\ln |c|}})=\log _{2}({\tfrac {\ln |z_{n}|}{\ln r}}).}
  

Für 
  
    
      
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        ≥
        r
        ≫
        2
      
    
    {\displaystyle |z_{n}(c)|\geq r\gg 2}
  
 erhält man die Fluchtzeit 
  
    
      
        t
        (
        c
        )
        =
        n
        −
        
          log
          
            2
          
        
        ⁡
        (
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
              
              
                ln
                ⁡
                r
              
            
          
        
        )
      
    
    {\displaystyle t(c)=n-\log _{2}({\tfrac {\ln |z_{n}(c)|}{\ln r}})}
  
 aus einer Korrektur (Normalisierung) der Iterationszahl 
  
    
      
        n
      
    
    {\displaystyle n}
  
. Die äquivalente Gleichung 
  
    
      
        
          
            
              1
              
                2
                
                  n
                
              
            
          
        
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        =
        
          2
          
            −
            t
            (
            c
            )
          
        
        ln
        ⁡
        r
      
    
    {\displaystyle {\tfrac {1}{2^{n}}}\ln |z_{n}(c)|=2^{-t(c)}\ln r}
  
 weist auf die enge Beziehung der Potentialfunktionen 
  
    
      
        
          Φ
          
            n
          
        
        (
        c
        )
        =
        
          
            
              1
              
                2
                
                  n
                
              
            
          
        
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
      
    
    {\displaystyle \Phi _{n}(c)={\tfrac {1}{2^{n}}}\ln |z_{n}(c)|}
  
 und 
  
    
      
        
          Ψ
          
            r
          
        
        (
        c
        )
        =
        
          2
          
            −
            t
            (
            c
            )
          
        
        ln
        ⁡
        r
      
    
    {\displaystyle \Psi _{r}(c)=2^{-t(c)}\ln r}
  
 hin. Häufig wird das Potential als Grenzwert 
  
    
      
        V
        (
        c
        )
        =
        
          lim
          
            n
            →
            ∞
          
        
        
          Φ
          
            n
          
        
        (
        c
        )
      
    
    {\displaystyle V(c)=\lim _{n\to \infty }\Phi _{n}(c)}
  
 definiert; die Potentialfunktion 
  
    
      
        
          Ψ
          
            e
          
        
        (
        c
        )
        =
        
          2
          
            −
            t
            (
            c
            )
          
        
        =
        
          e
          
            −
            t
            (
            c
            )
            ln
            ⁡
            2
          
        
      
    
    {\displaystyle \Psi _{e}(c)=2^{-t(c)}=e^{-t(c)\ln 2}}
  
 zur Eulerschen Zahl mit 
  
    
      
        ln
        ⁡
        e
        =
        1
      
    
    {\displaystyle \ln e=1}
  
 wird nach George Green mit 
  
    
      
        G
        (
        c
        )
      
    
    {\displaystyle G(c)}
  
 bezeichnet.
Nähert sich 
  
    
      
        c
      
    
    {\displaystyle c}
  
 dem Rand der Menge, so wird das Potential immer kleiner und erreicht am Rand den Grenzwert null. Daher kann der kürzeste Vektor 
  
    
      
        
          
            
              v
              →
            
          
        
      
    
    {\displaystyle {\vec {v}}}
  
 vom Punkt 
  
    
      
        c
      
    
    {\displaystyle c}
  
 zum Rand näherungsweise aus der Gleichung 
  
    
      
        
          Φ
          
            n
          
        
        (
        c
        )
        +
        
          
            
              v
              →
            
          
        
        ⋅
        ∇
        
          Φ
          
            n
          
        
        (
        c
        )
        =
        0
      
    
    {\displaystyle \Phi _{n}(c)+{\vec {v}}\cdot \nabla \Phi _{n}(c)=0}
  
 durch ein Skalarprodukt mit dem Gradientenvektor 
  
    
      
        
          
            
              u
              →
            
          
        
        =
        ∇
        
          Φ
          
            n
          
        
        (
        c
        )
      
    
    {\displaystyle {\vec {u}}=\nabla \Phi _{n}(c)}
  
 ermittelt werden. Man erhält so eine Lösung 
  
    
      
        
          
            
              v
              →
            
          
        
      
    
    {\displaystyle {\vec {v}}}
  
 mit der Länge 
  
    
      
        
          |
        
        
          
            
              v
              →
            
          
        
        
          |
        
        =
        
          
            
              
                
                  Φ
                  
                    n
                  
                
                (
                c
                )
              
              
                
                  |
                
                ∇
                
                  Φ
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
              
            
          
        
        =
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
              
              
                
                  |
                
                ∇
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
                
                  |
                
              
            
          
        
      
    
    {\displaystyle |{\vec {v}}|={\tfrac {\Phi _{n}(c)}{|\nabla \Phi _{n}(c)|}}={\tfrac {\ln |z_{n}(c)|}{|\nabla \ln |z_{n}(c)||}}}
  
. Aus 
  
    
      
        ∇
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        =
        
          
            
              1
              
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
              
            
          
        
        ∇
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
      
    
    {\displaystyle \nabla \ln |z_{n}(c)|={\tfrac {1}{|z_{n}(c)|}}\nabla |z_{n}(c)|}
  
 und 
  
    
      
        
          |
        
        ∇
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        
          |
        
        =
        
          |
        
        
          z
          
            n
          
        
        
          
          ′
        
        (
        c
        )
        
          |
        
      
    
    {\displaystyle |\nabla |z_{n}(c)||=|z_{n}\!'(c)|}
  
 folgt die Distanzformel:

  
    
      
        
          d
          
            n
          
        
        (
        c
        )
        =
        
          
            
              
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
              
              
                
                  |
                
                ∇
                ln
                ⁡
                
                  |
                
                
                  z
                  
                    n
                  
                
                (
                c
                )
                
                  |
                
                
                  |
                
              
            
          
        
        =
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        :
        
          |
        
        
          
            
              
                
                  z
                  
                    n
                  
                
                
                  
                  ′
                
                (
                c
                )
              
              
                
                  z
                  
                    n
                  
                
                (
                c
                )
              
            
          
        
        
          |
        
        =
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        (
        c
        )
        
          |
        
        ⋅
        
          |
        
        
          
            
              
                
                  z
                  
                    n
                  
                
                (
                c
                )
              
              
                
                  z
                  
                    n
                  
                
                
                  
                  ′
                
                (
                c
                )
              
            
          
        
        
          |
        
        .
      
    
    {\displaystyle d_{n}(c)={\tfrac {\ln |z_{n}(c)|}{|\nabla \ln |z_{n}(c)||}}=\ln |z_{n}(c)|:|{\tfrac {z_{n}\!'(c)}{z_{n}(c)}}|=\ln |z_{n}(c)|\cdot |{\tfrac {z_{n}(c)}{z_{n}\!'(c)}}|.}
  

Die ersten Polynome sind 
  
    
      
        
          z
          
            0
          
        
        (
        c
        )
        =
        0
      
    
    {\displaystyle z_{0}(c)=0}
  
, 
  
    
      
        
          z
          
            1
          
        
        (
        c
        )
        =
        c
      
    
    {\displaystyle z_{1}(c)=c}
  
, 
  
    
      
        
          z
          
            2
          
        
        (
        c
        )
        =
        
          c
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{2}(c)=c^{2}+c}
  
 und 
  
    
      
        
          z
          
            3
          
        
        (
        c
        )
        =
        (
        
          c
          
            2
          
        
        +
        c
        
          )
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{3}(c)=(c^{2}+c)^{2}+c}
  
, die zugehörigen Ableitungen sind 
  
    
      
        
          z
          
            0
          
        
        
          
          ′
        
        (
        c
        )
        =
        0
      
    
    {\displaystyle z_{0}\!'(c)=0}
  
, 
  
    
      
        
          z
          
            1
          
        
        
          
          ′
        
        (
        c
        )
        =
        1
      
    
    {\displaystyle z_{1}\!'(c)=1}
  
, 
  
    
      
        
          z
          
            2
          
        
        
          
          ′
        
        (
        c
        )
        =
        2
        c
        +
        1
      
    
    {\displaystyle z_{2}\!'(c)=2c+1}
  
 und 
  
    
      
        
          z
          
            3
          
        
        
          
          ′
        
        (
        c
        )
        =
        2
        (
        
          c
          
            2
          
        
        +
        c
        )
        (
        2
        c
        +
        1
        )
        +
        1
      
    
    {\displaystyle z_{3}\!'(c)=2(c^{2}+c)(2c+1)+1}
  
. Alle weiteren Folgenglieder ergeben sich aus der Iterationsvorschrift 
  
    
      
        
          z
          
            n
            +
            1
          
        
        (
        c
        )
        =
        
          z
          
            n
          
          
            2
          
        
        (
        c
        )
        +
        c
      
    
    {\displaystyle z_{n+1}(c)=z_{n}^{2}(c)+c}
  
 und der Ableitungsregel 
  
    
      
        
          z
          
            n
            +
            1
          
        
        
          

          
          ′
        
        (
        c
        )
        =
        2
        
          z
          
            n
          
        
        (
        c
        )
        
          z
          
            n
          
        
        
          
          ′
        
        (
        c
        )
        +
        1
      
    
    {\displaystyle z_{n+1}{}'(c)=2z_{n}(c)z_{n}\!'(c)+1}
  
.
Um die Distanzformel einzusehen, kann wieder die vereinfachte Folge 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        0
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+0}
  
 mit 
  
    
      
        
          z
          
            0
          
        
        =
        c
      
    
    {\displaystyle z_{0}=c}
  
 betrachtet werden. Die zugehörige Julia-Menge ist genau der Rand des Einheitskreises, und mit 
  
    
      
        
          z
          
            n
          
        
        =
        
          c
          
            
              2
              
                n
              
            
          
        
      
    
    {\displaystyle z_{n}=c^{2^{n}}}
  
 und 
  
    
      
        
          z
          
            n
          
        
        
          
          ′
        
        =
        
          2
          
            n
          
        
        
          c
          
            
              2
              
                n
              
            
            −
            1
          
        
      
    
    {\displaystyle z_{n}\!'=2^{n}c^{2^{n}-1}}
  
 ergibt sich

  
    
      
        
          d
          
            n
          
        
        =
        ln
        ⁡
        
          |
        
        
          c
          
            
              2
              
                n
              
            
          
        
        
          |
        
        ⋅
        
          |
        
        
          
            
              
                c
                
                  
                    2
                    
                      n
                    
                  
                
              
              
                
                  2
                  
                    n
                  
                
                
                  c
                  
                    
                      2
                      
                        n
                      
                    
                    −
                    1
                  
                
              
            
          
        
        
          |
        
        =
        
          2
          
            n
          
        
        ln
        ⁡
        
          |
        
        c
        
          |
        
        ⋅
        
          |
        
        
          
            
              c
              
                2
                
                  n
                
              
            
          
        
        
          |
        
        =
        ln
        ⁡
        
          |
        
        c
        
          |
        
        ⋅
        
          |
        
        c
        
          |
        
        ,
      
    
    {\displaystyle d_{n}=\ln |c^{2^{n}}|\cdot |{\tfrac {c^{2^{n}}}{2^{n}c^{2^{n}-1}}}|=2^{n}\ln |c|\cdot |{\tfrac {c}{2^{n}}}|=\ln |c|\cdot |c|,}
  

was in der Nähe des Einheitskreises eine gute Näherung für den Abstand 
  
    
      
        d
        =
        
          |
        
        c
        
          |
        
        −
        1
      
    
    {\displaystyle d=|c|-1}
  
 zum Rand ist: 
  
    
      
        
          |
        
        c
        
          |
        
        =
        1
        
          ,
        
        1
      
    
    {\displaystyle |c|=1{,}1}
  
 ergibt 
  
    
      
        
          d
          
            n
          
        
        =
        ln
        ⁡
        (
        1
        
          ,
        
        1
        )
        ⋅
        1
        
          ,
        
        1
        ≈
        0
        
          ,
        
        1
      
    
    {\displaystyle d_{n}=\ln(1{,}1)\cdot 1{,}1\approx 0{,}1}
  
 und 
  
    
      
        
          |
        
        c
        
          |
        
        =
        1
        
          ,
        
        01
      
    
    {\displaystyle |c|=1{,}01}
  
 ergibt 
  
    
      
        
          d
          
            n
          
        
        =
        ln
        ⁡
        (
        1
        
          ,
        
        01
        )
        ⋅
        1
        
          ,
        
        01
        ≈
        0
        
          ,
        
        01
      
    
    {\displaystyle d_{n}=\ln(1{,}01)\cdot 1{,}01\approx 0{,}01}
  
. Wenn der wirkliche Abstand zum Rand gleich 
  
    
      
        d
      
    
    {\displaystyle d}
  
 ist, so ergibt die Distanzformel den Wert

  
    
      
        
          d
          
            n
          
        
        =
        ln
        ⁡
        
          |
        
        
          z
          
            n
          
        
        
          |
        
        ⋅
        
          |
        
        
          
            
              
                z
                
                  n
                
              
              
                
                  z
                  
                    n
                  
                
                
                  
                  ′
                
              
            
          
        
        
          |
        
        ≈
        ln
        ⁡
        (
        1
        +
        d
        )
        ⋅
        (
        1
        +
        d
        )
        ≈
        d
        .
      
    
    {\displaystyle d_{n}=\ln |z_{n}|\cdot |{\tfrac {z_{n}}{z_{n}\!'}}|\approx \ln(1+d)\cdot (1+d)\approx d.}
  

Es stimmt aber nicht, dass der Grenzwert 
  
    
      
        
          lim
          
            n
            →
            ∞
          
        
        
          d
          
            n
          
        
      
    
    {\displaystyle \lim _{n\to \infty }d_{n}}
  
 gegen den wirklichen Abstand 
  
    
      
        d
      
    
    {\displaystyle d}
  
 konvergiert, tatsächlich gelten nur schwächere Ungleichungen. Die Abschätzung kann mithilfe der lambertschen W-Funktion jedoch etwas verbessern werden, und mit 
  
    
      
        W
        (
        
          d
          
            n
          
        
        )
        ≈
        ln
        ⁡
        (
        1
        +
        d
        )
      
    
    {\displaystyle W(d_{n})\approx \ln(1+d)}
  
 erhält man:

  
    
      
        
          d
          
            n
          
        
        ≈
        ln
        ⁡
        (
        1
        +
        d
        )
        
          e
          
            ln
            ⁡
            (
            1
            +
            d
            )
          
        
        
        ⟺
        
        d
        ≈
        
          e
          
            W
            (
            
              d
              
                n
              
            
            )
          
        
        −
        1.
      
    
    {\displaystyle d_{n}\approx \ln(1+d)e^{\ln(1+d)}\iff d\approx e^{W(d_{n})}-1.}
  


### Deep Zoom und Störungsrechnung
Prozessoren unterstützen in Hardware die Berechnung mit typischerweise 32-bit- und 64-bit-Gleitkommazahlen. Dies ermöglicht es, in die Mandelbrot-Menge etwa um den Faktor 1012 bis 1013 hineinzuzoomen, bevor benachbarte Pixel keine unterschiedlichen Startwerte 
  
    
      
        c
      
    
    {\displaystyle c}
  
 mehr annehmen können (und damit Fehler durch Blockartefakte offensichtlich werden).
Ab diesem Punkt reicht die Genauigkeit der in Hardware implementierten Arithmetik nicht mehr aus. Man ist gezwungen, auf in Software implementierte Routinen für hochgenaue Arithmetik auszuweichen, was bei aktuellen Prozessoren zu einem immensen Geschwindigkeitseinbruch führt.
Man kann diesen deutlich durch Nutzen von Störungsrechnung reduzieren. Diese Methode wurde 2013 von Kevin Martin beschrieben und 2021 von Zhuoran Yu um eine automatische Fehlerkorrektur erweitert.
Dazu wird ein (möglichst nicht divergenter) Referenzpunkt 
  
    
      
        c
      
    
    {\displaystyle c}
  
 in hochgenauer Arithmetik durchgerechnet und die zugehörige Referenzserie 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 in einfacher oder doppelter Genauigkeit aufgezeichnet:

  
    
      
        
          z
          
            0
          
        
        ,
         
        
          z
          
            1
          
        
        ,
         
        
          z
          
            2
          
        
        ,
         
        
          z
          
            3
          
        
        ,
        …
        ,
        
          z
          
            n
          
        
        .
      
    
    {\displaystyle z_{0},\ z_{1},\ z_{2},\ z_{3},\ldots ,z_{n}.}
  

Für in der direkten Nähe liegende Punkte 
  
    
      
        
          c
          ′
        
        =
        c
        +
        δ
      
    
    {\displaystyle c'=c+\delta }
  
 würde sich die Serie 
  
    
      
        
          z
          
            n
          
          ′
        
      
    
    {\displaystyle z_{n}'}
  
 ergeben:

  
    
      
        
          z
          
            0
          
          ′
        
        ,
         
        
          z
          
            1
          
          ′
        
        ,
         
        
          z
          
            2
          
          ′
        
        ,
         
        
          z
          
            3
          
          ′
        
        ,
        …
        ,
        
          z
          
            n
          
          ′
        
        .
      
    
    {\displaystyle z_{0}',\ z_{1}',\ z_{2}',\ z_{3}',\ldots ,z_{n}'.}
  

Man will nun die Abweichung 
  
    
      
        
          ϵ
          
            n
          
        
        =
        
          z
          
            n
          
          ′
        
        −
        
          z
          
            n
          
        
      
    
    {\displaystyle \epsilon _{n}=z_{n}'-z_{n}}
  
 indirekt ermitteln, um die direkte Berechnung der Serie 
  
    
      
        
          z
          
            n
          
          ′
        
      
    
    {\displaystyle z_{n}'}
  
 zu vermeiden. Dabei gilt:

  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
         
        +
        c
         
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}\ +c\ }
  
 für die Berechnung des Referenzpunktes und

  
    
      
        
          z
          
            n
            +
            1
          
          ′
        
        =
        
          z
          
            n
          
          ′
        
        
          
          
            2
          
        
        +
        
          c
          ′
        
         
      
    
    {\displaystyle z_{n+1}'=z_{n}'\!^{2}+c'\ }
  
 mit 
  
    
      
        
          z
          
            n
            +
            1
          
          ′
        
        =
        
          z
          
            n
            +
            1
          
        
        +
        
          ϵ
          
            n
            +
            1
          
        
      
    
    {\displaystyle z_{n+1}'=z_{n+1}+\epsilon _{n+1}}
  
 und 
  
    
      
        
          z
          
            n
          
          ′
        
        =
        
          z
          
            n
          
        
        +
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle z_{n}'=z_{n}+\epsilon _{n}}
  
 für die Berechnung des Nachbarpunktes.
Setzt man ein, ergibt sich

  
    
      
        
          z
          
            n
            +
            1
          
        
        +
        
          ϵ
          
            n
            +
            1
          
        
        =
        
          z
          
            n
            +
            1
          
          ′
        
        =
        
          z
          
            n
          
          ′
        
        
          
          
            2
          
        
        +
        
          c
          ′
        
        =
        
          z
          
            n
          
          ′
        
        
          
          
            2
          
        
        +
        (
        c
        +
        δ
        )
        =
        (
        
          z
          
            n
          
        
        +
        
          ϵ
          
            n
          
        
        
          )
          
            2
          
        
        +
        (
        c
        +
        δ
        )
        .
      
    
    {\displaystyle z_{n+1}+\epsilon _{n+1}=z_{n+1}'=z_{n}'\!^{2}+c'=z_{n}'\!^{2}+(c+\delta )=(z_{n}+\epsilon _{n})^{2}+(c+\delta ).}
  

Man betrachtet die beiden äußeren Terme

  
    
      
        
          z
          
            n
            +
            1
          
        
        
        +
        
          ϵ
          
            n
            +
            1
          
        
        =
        (
        
          z
          
            n
          
        
        +
        
          ϵ
          
            n
          
        
        
          )
          
            2
          
        
        +
        (
        c
        +
        δ
        )
         
      
    
    {\displaystyle z_{n+1}\quad +\epsilon _{n+1}=(z_{n}+\epsilon _{n})^{2}+(c+\delta )\ }
  
 und multipliziert aus

  
    
      
        
          z
          
            n
            +
            1
          
        
        
        +
        
          ϵ
          
            n
            +
            1
          
        
        =
        (
        
          z
          
            n
          
          
            2
          
        
        +
        2
        
          z
          
            n
          
        
        
          ϵ
          
            n
          
        
        +
        
          ϵ
          
            n
          
          
            2
          
        
        )
        +
        (
        c
        +
        δ
        )
         
      
    
    {\displaystyle z_{n+1}\quad +\epsilon _{n+1}=(z_{n}^{2}+2z_{n}\epsilon _{n}+\epsilon _{n}^{2})+(c+\delta )\ }
  
 und setzt 
  
    
      
        
          z
          
            n
            +
            1
          
        
        =
        
          z
          
            n
          
          
            2
          
        
        +
        c
      
    
    {\displaystyle z_{n+1}=z_{n}^{2}+c}
  

  
    
      
        
          z
          
            n
          
          
            2
          
        
         
        +
        c
        +
        
          ϵ
          
            n
            +
            1
          
        
        =
        (
        
          z
          
            n
          
          
            2
          
        
        +
        2
        
          z
          
            n
          
        
        
          ϵ
          
            n
          
        
        +
        
          ϵ
          
            n
          
          
            2
          
        
        )
        +
        (
        c
        +
        δ
        )
         
      
    
    {\displaystyle z_{n}^{2}\ +c+\epsilon _{n+1}=(z_{n}^{2}+2z_{n}\epsilon _{n}+\epsilon _{n}^{2})+(c+\delta )\ }
  
 und löst die Klammern auf und kürzt

  
    
      
        
        
          
            
              z
              
                n
              
              
                2
              
            
          
        
        
        +
        
        
          
            c
          
        
        
        +
        
          ϵ
          
            n
            +
            1
          
        
        =
        
        
          
            
              z
              
                n
              
              
                2
              
            
          
        
        
        +
        2
        
          z
          
            n
          
        
        
          ϵ
          
            n
          
        
        +
        
          ϵ
          
            n
          
          
            2
          
        
        +
        
        
          
            c
          
        
        
        +
        δ
         
      
    
    {\displaystyle \!{\cancel {z_{n}^{2}}}\!+\!{\cancel {c}}\!+\epsilon _{n+1}=\!{\cancel {z_{n}^{2}}}\!+2z_{n}\epsilon _{n}+\epsilon _{n}^{2}+\!{\cancel {c}}\!+\delta \ }
  
 und erhält

  
    
      
        
          ϵ
          
            n
            +
            1
          
        
        
        =
        
        2
        
          z
          
            n
          
        
        
          ϵ
          
            n
          
        
        +
        
          ϵ
          
            n
          
          
            2
          
        
        +
        δ
        .
      
    
    {\displaystyle \epsilon _{n+1}\quad =\quad 2z_{n}\epsilon _{n}+\epsilon _{n}^{2}+\delta .}
  

Wegen 
  
    
      
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{0}=0}
  
 und 
  
    
      
        
          z
          
            0
          
          ′
        
        =
        0
      
    
    {\displaystyle z_{0}'=0}
  
 ist auch 
  
    
      
        
          ϵ
          
            0
          
        
        =
        
          z
          
            0
          
          ′
        
        −
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle \epsilon _{0}=z_{0}'-z_{0}=0}
  
. Mit der abgespeicherten Referenzserie 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 können die Abweichungen 
  
    
      
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle \epsilon _{n}}
  
 der Reihe nach berechnet werden:

  
    
      
        
          ϵ
          
            1
          
        
        
        =
        
        2
        
          z
          
            0
          
        
        
          ϵ
          
            0
          
        
        +
        
          ϵ
          
            0
          
          
            2
          
        
        +
        δ
        =
        δ
        ,
      
    
    {\displaystyle \epsilon _{1}\quad =\quad 2z_{0}\epsilon _{0}+\epsilon _{0}^{2}+\delta =\delta ,}
  

  
    
      
        
          ϵ
          
            2
          
        
        
        =
        
        2
        
          z
          
            1
          
        
        
          ϵ
          
            1
          
        
        +
        
          ϵ
          
            1
          
          
            2
          
        
        +
        δ
        ,
      
    
    {\displaystyle \epsilon _{2}\quad =\quad 2z_{1}\epsilon _{1}+\epsilon _{1}^{2}+\delta ,}
  

  
    
      
        
          ϵ
          
            3
          
        
        
        =
        
        2
        
          z
          
            2
          
        
        
          ϵ
          
            2
          
        
        +
        
          ϵ
          
            2
          
          
            2
          
        
        +
        δ
        ,
      
    
    {\displaystyle \epsilon _{3}\quad =\quad 2z_{2}\epsilon _{2}+\epsilon _{2}^{2}+\delta ,}
  

  
    
      
         
        ⋮
      
    
    {\displaystyle \ \vdots }
  

  
    
      
        
          ϵ
          
            n
          
        
        
        =
        
        2
        
          z
          
            n
            −
            1
          
        
        
          ϵ
          
            n
            −
            1
          
        
        +
        
          ϵ
          
            n
            −
            1
          
          
            2
          
        
        +
        δ
        .
      
    
    {\displaystyle \epsilon _{n}\quad =\quad 2z_{n-1}\epsilon _{n-1}+\epsilon _{n-1}^{2}+\delta .}
  

Solange diese Terme klein bleiben, verhält sich 
  
    
      
        
          z
          
            n
          
          ′
        
        =
        
          z
          
            n
          
        
        +
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle z_{n}'=z_{n}+\epsilon _{n}}
  
 ähnlich wie 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
. Insbesondere bleibt die Serie 
  
    
      
        
          z
          
            n
          
          ′
        
      
    
    {\displaystyle z_{n}'}
  
 beschränkt, wenn die Referenzserie 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
 und die Abweichungen 
  
    
      
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle \epsilon _{n}}
  
 beschränkt sind. Wenn die Abweichungen jedoch zu groß werden, kann das zu Fehlern in der Darstellung (Glitches) führen. Im Prinzip kann man dann auf einen günstigeren Referenzpunkt innerhalb der Referenzserie ausweichen: Liegt ein Punkt 
  
    
      
        z
        =
        
          z
          
            i
          
        
        +
        
          ϵ
          
            i
          
        
      
    
    {\displaystyle z=z_{i}+\epsilon _{i}}
  
 näher am Referenzpunkt 
  
    
      
        
          z
          
            j
          
        
      
    
    {\displaystyle z_{j}}
  
 als am aktuellen Referenzpunkt 
  
    
      
        
          z
          
            i
          
        
      
    
    {\displaystyle z_{i}}
  
, so kann er auch durch 
  
    
      
        z
        =
        
          z
          
            j
          
        
        +
        
          ϵ
          
            j
          
        
      
    
    {\displaystyle z=z_{j}+\epsilon _{j}}
  
 dargestellt und die neue Abweichung aus der Differenz 
  
    
      
        
          ϵ
          
            j
          
        
        =
        z
        −
        
          z
          
            j
          
        
      
    
    {\displaystyle \epsilon _{j}=z-z_{j}}
  
 berechnet werden. Allerdings verliert man bei einer solchen Subtraktion aufgrund der Auslöschung die notwendige Genauigkeit. Nur beim ersten Referenzpunkt 
  
    
      
        
          z
          
            j
          
        
        =
        
          z
          
            0
          
        
        =
        0
      
    
    {\displaystyle z_{j}=z_{0}=0}
  
 passiert das nicht, und die Abweichung ergibt sich zu 
  
    
      
        
          ϵ
          
            j
          
        
        =
        z
        −
        0
        =
        z
      
    
    {\displaystyle \epsilon _{j}=z-0=z}
  
. Man kann also jederzeit zum Anfang der Referenzserie zurückkehren, wenn der Abstand zu 
  
    
      
        
          z
          
            0
          
        
      
    
    {\displaystyle z_{0}}
  
 kleiner als der Abstand zu 
  
    
      
        
          z
          
            i
          
        
      
    
    {\displaystyle z_{i}}
  
 ist, d. h. wenn 
  
    
      
        
          |
        
        z
        
          |
        
        <
        
          |
        
        
          ϵ
          
            i
          
        
        
          |
        
      
    
    {\displaystyle |z|<|\epsilon _{i}|}
  
 ist. Durch dieses gelegentliche Zurücksetzen der Referenzserie (Rebasing) werden die Abweichungen 
  
    
      
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle \epsilon _{n}}
  
 möglichst klein gehalten.

Praktische Anwendung

Das Programm „Kalles Fraktaler“ basierte auf diesem Ansatz, und neuere Entwicklungen haben zu einer Generation noch leistungsfähigerer Fraktalgeneratoren geführt.

Probleme

Der verwendete Referenzpunkt sollte möglichst innerhalb der zu untersuchenden Iterationen nicht divergieren, da die Nachbarpunkte maximal bis zu dem Abbruch der Referenzserie untersucht werden können. Alle Punkte mit einer größeren Iterationszahl müssen erneut berechnet werden. Die Auswahl und das Verwerfen von sich als ungeeignet herausgestellten Referenzpunkten stellt eine Herausforderung dar.
Für die Berechnung der Abweichungen 
  
    
      
        
          ϵ
          
            n
          
        
      
    
    {\displaystyle \epsilon _{n}}
  
 wird keine hohe Genauigkeit benötigt, allerdings ist der Wertebereich von 64-bit-Gleitkommazahlen arg limitiert. Ab etwa 10300 muss man leider mit eigenen Mathematikroutinen arbeiten, die zwar nur eine 52/64 bit-Mantisse benötigen, aber deutlich langsamer als hartkodierte Gleitkomma-Operationen sind (die aber wesentlich schneller als Berechnungen mit langen Mantissen sind).
Man benötigt Speicherbandbreite für das Lesen der Referenzserie 
  
    
      
        
          z
          
            n
          
        
      
    
    {\displaystyle z_{n}}
  
. Die erste Verlangsamung gibt es, wenn diese nicht mehr in den L2-Cache passt, sehr viel langsamer wird es, wenn der L3-Cache nicht mehr ausreicht und auf den RAM zugegriffen werden muss.
Es gibt Abweichungen zwischen der nativen Iteration mit 
  
    
      
        
          c
          ′
        
      
    
    {\displaystyle c'}
  
 und der Störungsrechnung basierend auf 
  
    
      
        
          c
          ′
        
        =
        c
        +
        δ
      
    
    {\displaystyle c'=c+\delta }
  
. Sieht man sich das ganze genauer an, ist die Berechnung basierend auf 
  
    
      
        
          c
          ′
        
        =
        c
        +
        δ
      
    
    {\displaystyle c'=c+\delta }
  
 meistens genauer, wenn die Referenzpunktberechnung hinreichend genau war (was mit Langzahlarithmetik für diesen einen Wert aber problemlos machbar ist).


### Berechnungsgeschwindigkeit
Berechnungen nahmen beim Stand der Technik Ende der 1980er Jahre viel Zeit in Anspruch. Vergleichsweise zeigt dies eine Aufstellung durch die Anzahl der Iterationen, die verschiedene CPUs pro Sekunde durchführen konnten.

Eingesetzte Software

BASIC-Interpreter des KC 85/3
Selbst geschriebener Assemblercode
Borland C++ mit integriertem Assembler, die Gleitkomma-Emulation dieses Compilers war vergleichsweise langsam
Visual Studio C++ mit Intel Intrinsics, Nutzung von AVX2 und Fused multiply-add, multithreaded, 2 Threads pro Core
Eingesetzte Prozessoren

Z80: 8-Bit-Prozessor ohne Hardware-Multiplikation
80386: 32-Bit-Prozessor mit Hardware-Multiplikation und Gleitkomma-Koprozessor
Pentium: 32-Bit-Prozessor mit Hardware-Multiplikation und mit eingebautem Gleitkomma-Koprozessor
Xeon E5-2623: Vektorbefehle, 4 Kerne/8 Threads pro Sockel


## Rezeption in der Öffentlichkeit
Außerhalb der Fachwelt wurde die Mandelbrot-Menge vor allem durch den ästhetischen Wert der Computergrafiken bekannt, der durch künstlerische Farbgestaltung des Außenbereichs, der nicht zur Menge gehört, unterstützt wird. Sie erlangte durch Publikationen von Bildern in den Medien Ende der 1980er Jahre einen für ein mathematisches Thema dieser Art ungewöhnlich großen Bekanntheitsgrad und dürfte bis heute (Stand 2023) das populärste Fraktal sein, sowie auch eines der populärsten Objekte der zeitgenössischen Mathematik.
Die Mandelbrot-Menge wird als das formenreichste geometrische Gebilde bezeichnet. Sie hat Computerkünstler inspiriert und zu einem Aufschwung fraktaler Konzepte beigetragen. Dabei finden zahlreiche Modifikationen des Algorithmus Anwendung, welcher der Mandelbrot-Menge zugrunde liegt.
Ein weiterer Aspekt ist der extreme Kontrast zwischen diesem und der Einfachheit des zugrunde liegenden Algorithmus, der an biologische Systeme erinnert, bei denen nach naturwissenschaftlicher Sicht ebenfalls aus einer vergleichsweise geringen Zahl von Regeln äußerst komplexe Systeme entstehen können, sowie die Nähe zur Chaosforschung, die ebenfalls in der Öffentlichkeit großes Interesse geweckt hatte.
Die Bezeichnung ,Apfelmännchen‘ leitet sich von der geometrischen Grobform einer um 90 Grad im Uhrzeigersinn gedrehten Mandelbrot-Menge her.
Der US-amerikanische Musiker Jonathan Coulton hat ein Lied über die Mandelbrot-Menge veröffentlicht, in dem Benoît Mandelbrot dafür gedankt wird, dass er Ordnung in das Chaos gebracht habe.


## Verwandte Themen
Das Buddhabrot ist eine gefilterte Dichteverteilung der Folgenglieder, bei der alle beschränkten Folgen ignoriert werden.
Die Mandelknolle und die Mandelbox sind mögliche Verallgemeinerungen der Mandelbrot-Menge für drei- oder mehrdimensionale Räume.


## Anmerkungen
## Literatur
Benoît Mandelbrot: Die fraktale Geometrie der Natur. ISBN 3-7643-2646-8.
John Briggs, F. David Peat: Die Entdeckung des Chaos. ISBN 3-446-15966-5.
Heinz-Otto Peitgen, Peter H. Richter: The Beauty of Fractals. ISBN 0-387-15851-0.
Heinz-Otto Peitgen, Dietmar Saupe: The Science of Fractal Images. ISBN 0-387-96608-0.
Karl Günter Kröber: Das Märchen vom Apfelmännchen – 1. Wege in die Unendlichkeit. ISBN 3-499-60881-2.
Karl Günter Kröber: Das Märchen vom Apfelmännchen – 2. Reise durch das malumitische Universum. ISBN 3-499-60882-0.
Dierk Schleicher: On Fibers and Local Connectivity of Mandelbrot and Multibrot Sets. In: M.Lapidus, M. van Frankenhuysen (eds): Fractal Geometry and Applications: A Jubilee of Benoît Mandelbrot. Proceedings of Symposia in Pure Mathematics 72, American Mathematical Society (2004), S. 477–507, 1999, PDF; 0,4 MB.


## Weblinks
Eric W. Weisstein: Mandelbrot Set. In: MathWorld (englisch).
Robert P. Munafo: The Encyclopedia of the Mandelbrot Set. In: mrob.com. 1996; abgerufen am 1. November 2024 (englisch). 
Robert Devaney: Unveiling the Mandelbrot Set. In: plus.maths.org. 2006; abgerufen am 15. Juli 2021 (englisch). 
Wolfgang Beyer: Animationen zur Zoomfahrt im hiesigen Artikel. In: wolfgangbeyer.de. 2007; abgerufen am 15. Juli 2021. 
Anders Sandberg: Mercator Mandelbrot Maps. In: flickr.com. 2009; abgerufen am 15. Juli 2021 (englisch). 
Norbert Treitz: Mandelbrots Apfelmännchen. In: spektrum.de. 2017; abgerufen am 7. August 2021. 
Albert Lobo: Mandelbrot viewer. In: albertlobo.com. 2019; abgerufen am 1. November 2024 (englisch). 
Joao Maio: Mandelbrot Maps. In: github.com. 2021; abgerufen am 10. Juli 2023 (englisch). 
Zooming in a Mandelbrot Fractal by a Ratio of 10^227. In: youtube.com. 2016; abgerufen am 14. Oktober 2021. 
Sounds of the Mandelbrot Set. In: youtube.com. 2021; abgerufen am 15. Juli 2021 (englisch). 
Diverse Algorithmen zur Berechnung der Mandelbrot-Menge. In: Rosetta Code. Abgerufen am 15. Juli 2021 (englisch). 


## Einzelnachweise