# Ensemble de Mandelbrot

En mathématiques, l'ensemble de Mandelbrot est une fractale définie comme l'ensemble des points c du plan complexe pour lesquels la suite de nombres complexes définie par récurrence par :
    ${\begin{cases}z_{0}=0\\z_{n+1}=z_{n}^{2}+c\end{cases}}$
est bornée.
L'ensemble de Mandelbrot a été découvert par Gaston Julia et Pierre Fatou avant la Première Guerre mondiale. Sa définition et son nom actuel sont dus à Adrien Douady, en hommage aux représentations qu'en a réalisées Benoît Mandelbrot dans les années 1980. Cet ensemble permet d'indicer les ensembles de Julia associés à la suite : à chaque point du plan complexe correspond un ensemble de Julia différent. Les points de l'ensemble de Mandelbrot correspondent précisément aux ensembles de Julia connexes. Cet ensemble est donc intimement lié aux ensembles de Julia, ils produisent d'ailleurs des formes similairement complexes.
Les images de l'ensemble de Mandelbrot sont réalisées en parcourant les nombres complexes sur une région carrée du plan complexe et en déterminant pour chacun d'eux si le résultat tend vers l'infini ou pas lorsqu'on y itère une opération mathématique. On considère la partie réelle et imaginaire de chaque nombre complexe comme des coordonnées et chaque pixel est coloré selon la rapidité de divergence, ou si elle ne diverge pas.
Les images de l'ensemble de Mandelbrot exposent une limite élaborée qui révèle progressivement des détails récursifs toujours plus fins en augmentant le grossissement. La limite de l'ensemble est constituée de plus petites versions de la forme principale, donc la propriété fractale de l'autosimilarité s'applique à l'ensemble tout entier (et pas simplement à certaines parties).
L'ensemble de Mandelbrot est devenu populaire hors des mathématiques, comme inspiration artistique et comme exemple de structure complexe venant de l'application de règles simples. C'est l'un des exemples les plus connus de visualisation mathématique.

## Historique

L'ensemble de Mandelbrot tire ses origines de la dynamique complexe, un domaine défriché par les mathématiciens français Pierre Fatou et Gaston Julia au début du XXe siècle.
La première représentation de cet ensemble apparaît en 1978 dans un article de Robert W. Brooks et Peter Matelski. Le 1er mars 1980, au centre Thomas J. Watson de recherche IBM (dans l'État de New York), Benoît Mandelbrot obtient pour la première fois, une visualisation par ordinateur de cet ensemble. Mandelbrot étudie l'espace des paramètres des polynômes quadratiques complexes dans un article publié en 1980.
En 1984, l'étude de l'ensemble de Mandelbrot commence réellement avec les travaux d'Adrien Douady et de John H. Hubbard, qui établissent ses propriétés fondamentales et baptisent l'ensemble en l'honneur de Mandelbrot.
En 1985, les mathématiciens Heinz-Otto Peitgen (en) et Peter Richter popularisent l'ensemble de Mandelbrot par des images de qualité qui frappent les esprits$^,$$^,$. Dans le numéro d'août 1985 du magazine Scientific American, l'ensemble de Mandelbrot est présenté au grand public comme « l'objet mathématique le plus complexe jamais découvert » et présente l'algorithme qui permet de le tracer soi-même. La couverture de ce numéro reprend une image créée par Peitgen$^,$.
Les travaux de Douady et Hubbard coïncidaient avec un intérêt considérable pour la dynamique complexe et l'étude de l'ensemble de Mandelbrot a été le centre d'attention de ce domaine depuis. Parmi les mathématiciens qui apportèrent une contribution significative à l'étude de cet ensemble, il faut citer Tan Lei, Mikhail Lyubich, Curtis Tracy McMullen, John Milnor, Mitsuhiro Shishikura (en) et Jean-Christophe Yoccoz.

## Propriétés
### Définition

La théorie générale, développée par Pierre Fatou et Gaston Julia au début du XXe siècle, associe à toute fonction (suffisamment régulière) f(z, c) (à arguments et valeurs complexes) les ensembles de Julia J$_c$, définis (pour un c fixé) comme la frontière de l'ensemble des a complexes tels que la suite définie par z$_0$ = a et z$_{n+1}$ = f(z$_n$, c) reste bornée (en module) ; pour la fonction particulière f(z, c) = z$^2$ + c, on définit l'ensemble de Mandelbrot M comme l'ensemble des c pour lequel J$_c$ est connexe. Fatou et Julia ont démontré que cette définition est équivalente à celle donnée dans l'introduction, c'est-à-dire que c appartient à M si et seulement si la suite (z$_n$) définie par z$_0$ = 0 et z$_{n+1}$ = z$_n$$^2$ + c reste bornée (en module) ; restant dans les réels, il s'agit donc des points de coordonnées (a,b) tels que les deux suites (x$_n$) et (y$_n$) définies par la récurrence x$_0$ = y$_0$ = 0 et x$_{n+1}$ = x$_n$$^2$ – y$_n$$^2$ + a ; y$_{n+1}$ = 2x$_n$y$_n$ + b restent bornées.

### Barrière du module égal à 2

Si la suite des modules des z$_n$ est strictement supérieure à 2 pour un certain indice alors, cette suite est croissante à partir de cet indice, et elle tend vers l'infini.
Par inégalité triangulaire,
    |z$_{n+1}$| = |z$_n$$^2$ + c| ≥ |z$_n$$^2$| – |c| (*).
Soit α la racine positive de l'équation α = α$^2$ – |c|. Alors, α ↦ |c| est une bijection croissante de [1, +∞[ sur [0, +∞[, et (|c| > α ⇔ α > 2 ⇔ |c| > 2).
En posant
    x$_n$ = |z$_n$| – α,
(*) devient α + x$_{n+1}$ ≥ (α + x$_n$)$^2$ – |c| = α + 2αx$_n$ + x$_n$$^2$, d'où
    x$_{n+1}$ ≥ 2αx$_n$.
Par conséquent si, pour un certain indice k, x$_k$ > 0 alors, à partir de cet indice, la suite (x$_n$) croît plus vite qu'une suite géométrique.
Cette condition équivaut à |z$_k$| > α donc a lieu en particulier :
* pour |c| > 2, dès que k = 1 (car |c| > 2 ⇒ |z$_1$| = |c| > α) ;
* mais aussi, pour |c| ≤ 2, dès que |z$_k$| > 2 (car |c| ≤ 2 ⇒ α ≤ 2).
Pour que la suite (z$_n$) soit bornée, il suffit donc qu'elle ne tende pas vers l'infini, et il faut qu'elle reste bornée par 2 (en module). Ceci fournit deux définitions de l'ensemble de Mandelbrot équivalentes à celle donnée en introduction, et prouve de plus que cet ensemble est inclus dans le disque fermé de centre 0 et de rayon 2.

### Géométrie élémentaire

L'ensemble M de Mandelbrot est compact, symétrique par rapport à l'axe réel et contient le disque fermé de centre 0 et de rayon 1/4.
Son intersection avec l'axe réel est le segment [–2, 1/4].
Notons (z$_n$(c)) la suite associée à un complexe c et pour tout entier n > 0, M$_n$ l'ensemble des points c pour lesquels |z$_n$(c)| ≤ 2. D'après le paragraphe ci-dessus, M est l'intersection de la suite des M$_n$ (on peut, accessoirement, en déduire aussi que cette suite est décroissante). Or M$_1$ est un disque fermé donc un compact du plan, et chaque M$_n$ est fermé, comme image réciproque de M$_1$ par la fonction polynomiale c ↦ z$_n$(c). Il en résulte que M est compact.
Les fonctions c ↦ z$_n$(c) étant des polynômes à coefficients réels, tous les M$_n$ ainsi que M sont stables par conjugaison complexe, ce qui se traduit géométriquement par la symétrie des M$_n$ et de M par rapport à l'axe réel.
Si |c| ≤ 1/4 alors, par récurrence, |z$_n$(c)| ≤ 1/2 pour tout n donc c appartient à M.
Il reste à élucider la situation si c est un réel tel que 1/4 < |c| ≤ 2.
* si –2 ≤ c ≤ 0 (autrement dit si c$^2$ + c ≤ – c) alors, par récurrence, c ≤ z$_n$(c) ≤ – c pour tout n donc c appartient à M.
* pour tout réel c, z$_{n+1}$(c) – z$_n$(c) ≥ c – 1/4 donc si c > 1/4 alors c n'appartient pas à M.
Remarque. Toutes les racines des polynômes c ↦ z$_n$(c) sont dans M. En effet, si z$_r$(c) = 0 alors la suite des z$_n$(c) est de période r donc bornée.
Son aire est estimée autour de 1,50659177 ± 0,00000008.
Son centre de gravité est estimé à la position réelle –0,286 768 3 ± 0,000 000 1.
La principale structure remarquable est la cardioïde centrale, de centre $A\left({\frac {1}{4}},0\right)$ et d'équation polaire $\rho (\theta )={\frac {1}{2}}\left(1-\cos \theta \right)$. On peut aussi citer le disque de centre $B(-1,0)$ et de rayon $1/4$.

### Connexité

Douady et Hubbard ont montré, en 1985, que l'ensemble est simplement connexe. Ce résultat n'était pas évident à l'observation des premiers tracés de Mandelbrot, qui faisait apparaître des « îlots » semblant détachés du reste. Pour ce faire, ils ont montré que le complémentaire de l'ensemble de Mandelbrot est conformément isomorphe au complémentaire dans ℂ du disque unité.
On conjecture que l'ensemble de Mandelbrot est localement connexe$^,$$^,$.

### Autosimilarité

L'ensemble de Mandelbrot est autosimilaire dans le voisinage des points de Misiurewicz. Ces points sont denses sur toute la frontière de l'ensemble. On conjecture qu'il est également autosimilaire, à la limite, autour des points de Feigenbaum (ex. : −1,401 155 ou −0,152 8 + 1,039 7 i).
Des versions réduites de l'ensemble de Mandelbrot apparaissent sur toute sa frontière, jusqu'à des grossissements infinis, avec de légères différences, mais l'ensemble de Mandelbrot n'est pas, en général, strictement autosimilaire.

### Universalité

L'ensemble de Mandelbrot M a un caractère universel pour de nombreuses fonctions holomorphes. Des copies de M sont visibles sur les frontières de leurs bassins d'attraction, c'est-à-dire des ensembles des c pour lesquels les itérés f(f(…f(z)))) convergent vers un complexe donné. En voici quelques exemples, avec des fonctions transcendantes :
* $z\mapsto \cos(z)+c$ (fonction cosinus)
* $z\mapsto {\rm {e}}^{-z^{2}}+c$ (fonction de Gauss)
On retrouve également M lors de l'itération d'une famille de fonctions cubiques telles que $z\mapsto z^{3}+(\lambda -1)z-\lambda$, par la méthode de Newton. L'ensemble des points $\lambda$ ne convergeant pas vers une racine de ce polynôme prend la forme de M.
Là encore, l'autosimilarité n'est pas stricte (il n’y a pas de similitude exacte avec M), mais il existe une déformation de faible amplitude entre M et ces copies.
* L'ensemble des points ne convergeant pas, pour une fonction cubique, par la méthode de Newton.
* Points de non-convergence pour les itérées de f(z,c)=cos(z)+1/c, région [-0.09:-0.07]x[0.345:0.36].

### Dimension de Hausdorff

La dimension de Hausdorff de la frontière de l'ensemble de Mandelbrot vaut 2. Ce résultat a été démontré en 1990 par Mitsuhiro Shishikura. On ne sait pas si cette frontière a une aire (définie par la mesure de Lebesgue) positive.

### Lien avec l'équation logistique

Les paramètres de M, dans l'intervalle réel [–2, 1/4], peuvent être mis en correspondance bijective avec ceux de l'équation logistique : $z\mapsto \lambda z(1-z),\quad \lambda \in [1,4],$ la correspondance étant donnée par : $c={\frac {\lambda }{2}}\left(1-{\frac {\lambda }{2}}\right).$

### Lien avec les ensembles de Julia

L'ensemble de Mandelbrot M peut être défini comme l'ensemble des complexes c pour lesquels l'ensemble de Julia correspondant, J$_c$, est connexe. Mais, plus précisément, on a (à la limite) identité entre J$_c$ et M au voisinage de c, lorsque c appartient à la frontière de M. Ainsi, il y a une correspondance proche entre la géométrie de l'ensemble de Mandelbrot à un point donné et la structure de l'ensemble de Julia lui correspondant. Par exemple, un point est dans l'ensemble de Mandelbrot exactement là où l'ensemble de Julia correspondant est connexe.
Ce principe est exploité dans beaucoup de preuves et découvertes profondes sur l'ensemble de Mandelbrot. Par exemple, Mitsuhiro Shishikura a prouvé que la frontière de l'ensemble de Mandelbrot possède une dimension de Hausdorff égale à 2 et que l'ensemble de Julia correspondant à un point de cette frontière a aussi une dimension égale à 2.
Jean-Christophe Yoccoz a prouvé que l'ensemble de Julia est un espace localement connexe pour certains paramètres, avant de l'établir également pour l'ensemble de Mandelbrot aux points correspondant à ces paramètres (mais la connexité locale de M en tout point reste une conjecture ouverte en 2024).

### Bourgeons, antennes et périodicités

L'ensemble de Mandelbrot fait apparaître nombre de structures en forme de bourgeons entourant une structure principale en forme de cardioïde.
* La cardioïde est l'ensemble des points c qui convergent vers un point fixe (période 1). Ce sont les points de la forme $c={\frac {\mu }{2}}\left(1-{\frac {\mu }{2}}\right)$ pour tout $\mu$ appartenant au disque unité.
* Le bourgeon principal, à gauche de la cardioïde, lui est attaché au point c = –3/4. Il s'agit d'un disque centré en c = –1 et de rayon 1/4, formé de l'ensemble des points paramètres qui, à la limite, convergent vers un cycle de période 2.
* Les autres bourgeons tangents à la cardioïde sont les points admettant d'autres périodicités. Pour tout nombre rationnel p/q, avec p et q premiers entre eux, il existe un bourgeon presque circulaire tangent à la cardioïde au point $c_{p/q}={\frac {e^{2\pi ip/q}}{2}}\left(1-{\frac {e^{2\pi ip/q}}{2}}\right)$. Chacun de ces bourgeons est appelé « bourgeon p/q ». Il s'agit de l'ensemble des points paramètres convergeant vers un cycle de période q.
* Enfin, chaque bourgeon porte lui-même des bourgeons, représentatifs d'une périodicité différente, selon le même schéma. Par exemple, le bourgeon à gauche du grand bourgeon de période 2 a une périodicité de période 4. Le bourgeon immédiatement à sa gauche est de période 8, puis 16, etc., suivant ainsi le motif de doublement de période du diagramme de bifurcation de l'équation logistique.
Les bourgeons sont également surmontés de filaments en forme d'antenne. Le nombre d'antennes est directement lié à la périodicité du bourgeon. Ainsi, compter le nombre d'antennes permet de déterminer la périodicité du bourgeon.

## Dessiner l'ensemble

On a vu plus haut que dès que le module d'un z$_n$ est strictement plus grand que 2, la suite diverge vers l'infini, et donc c est en dehors de l'ensemble de Mandelbrot. Cela permet d'arrêter le calcul pour les points ayant un module strictement supérieur à 2 et qui sont donc en dehors de l'ensemble de Mandelbrot. Pour les points de l'ensemble de Mandelbrot, le calcul n'arrivera jamais à terme, donc il doit être arrêté après un certain nombre d'itérations déterminé par le programme. Il en résulte que l'image affichée n'est qu'une approximation du vrai ensemble.
Bien que cela n'ait aucune importance sur le plan mathématique, la plupart des programmes générant des fractales affichent les points en dehors de l'ensemble de Mandelbrot dans différentes couleurs. La couleur attribuée à un point n'appartenant pas à l'ensemble dépend du nombre d'itérations au bout desquelles la suite correspondante est déclarée divergente vers l'infini (par exemple quand le module est strictement supérieur à 2). Cela donne plusieurs zones concentriques, qui entourent l'ensemble de Mandelbrot. Les plus éloignées sont constituées des points c pour lesquels la suite (z$_n$) tend « plus rapidement » vers l'infini. Ces différentes zones délimitent d'une manière plus ou moins précise l'ensemble de Mandelbrot.
Le nombre maximum d'itérations influe donc fortement sur la représentation de l'ensemble de Mandelbrot. En effet si la forme de base de l'ensemble ne varie pas outre mesure à partir de 100 itérations il n'en est pas de même pour tous les centres et tous les zooms de cet ensemble.

### Optimisations
#### Cardioïde / bourgeon principal

Un moyen de limiter les calculs est de trouver au préalable les points appartenant à la cardioïde et au bourgeon principal, c'est-à-dire au disque de centre (–1, 0) et de rayon 1/4. Ces points coûtent particulièrement cher en calculs, puisqu'ils appartiennent à l'ensemble de Mandelbrot, et nécessitent donc d'aller jusqu'au maximum d'itérations. Avant de passer un point de coordonnées $x+iy$ dans l’algorithme itératif, il faut donc vérifier :
* si $x<p-2p^{2}+{\frac {1}{4}}$ avec $p={\sqrt {\left(x-{\frac {1}{4}}\right)^{2}+y^{2}}}$alors le point appartient à la cardioïde, donc à l'ensemble de Mandelbrot ;
* si $(x+1)^{2}+y^{2}<{\frac {1}{16}}$alors le point appartient au disque de centre (–1, 0) et de rayon 1/4, c'est-à-dire au bourgeon principal, et donc à l'ensemble de Mandelbrot.
Les autres bourgeons n'ont pas de test équivalent, puisqu'ils ne sont pas parfaitement circulaires. Il est néanmoins possible de vérifier si un point appartient à un cercle que l'on sait inscrit dans un de ces bourgeons, permettant à la plupart des points de ces derniers de ne pas lancer l'algorithme itératif.

#### Algorithme par perturbation

Il est possible de reformuler l'équation principale, de manière à raisonner sur la différence d'évolution par rapport à un point voisin. Cette méthode permet de ramener ces différences $\Delta _{n}$ à une faible précision (vingtaine de chiffres significatifs), dans le domaine des unités de calcul en virgule flottante matérielles de l'ordinateur, qui sont beaucoup plus rapides que les bibliothèques d'arithmétique multiprécision qui ne sont pas autant accélérées matériellement. L'ensemble est ainsi calculé à partir d'un (ou plusieurs) points de référence, calculés de manière absolue.
Si on raisonne en absolu, les nombres mis en œuvre pour le calcul de l'ensemble peuvent atteindre des milliers de chiffres significatifs pour des grands niveaux de zoom, et s'accroissant davantage avec le zoom. Avec cet algorithme, la vitesse de calcul reste relativement constante selon le niveau de zoom, la vitesse de calcul dépendant plus de la complexité de la scène.
Cet algorithme est à la base d'un des programmes de calcul le plus rapide de l'ensemble de Mandelbrot en 2017 : Kalles Fraktaler.
Principe :
$X_{n+1}=X_{n}^{2}+X_{0}$
$Y_{n+1}=Y_{n}^{2}+Y_{0}$
$\Delta _{n}=Y_{n}-X_{n}$
${\begin{aligned}\Delta _{n+1}=Y_{n+1}-X_{n+1}&=(Y_{n}^{2}+Y_{0})-(X_{n}^{2}+X_{0})\\&=((X_{n}+\Delta _{n})^{2}+X_{0}+\Delta _{0})-(X_{n}^{2}+X_{0})\\&=(X_{n}^{2}+2X_{n}\Delta _{n}+\Delta _{n}^{2}+X_{0}+\Delta _{0})-(X_{n}^{2}+X_{0})\\&=2X_{n}\Delta _{n}+\Delta _{n}^{2}+\Delta _{0}\end{aligned}}$
Tous les $\Delta _{n}$, calculés à partir des $X_{n}$ de référence mémorisés, ont une faible précision, permettant d'utiliser les unités de calcul en virgule flottante matérielles de l'ordinateur pour les opérations les plus nombreuses et coûteuses.
Le calcul peut encore être davantage optimisé en exprimant $\Delta _{n}$ en utilisant l'approximation :
$\Delta _{n}=A_{n}\Delta _{0}+B_{n}\Delta _{0}^{2}+C_{n}\Delta _{0}^{3}+o(\Delta _{0}^{4})$
avec $A_{n},B_{n},C_{n}$ ne dépendant pas de $\Delta _{n}$, mais seulement de $X_{n}$ (série du point de référence à partir duquel on calcule les perturbations).
Donc, une fois calculés et mémorisés tous les $A_{n},B_{n},C_{n}$ pour le point de référence ($\Delta _{0}=0$), on peut rechercher directement pour tout point voisin (en faisant varier $\Delta _{0}$), par recherche dichotomique, le $\Delta _{n}$ qui provoque la divergence de la suite (puisqu'on n'a pas besoin du $\Delta _{n-1}$ pour calculer $\Delta _{n}$) et le calcul d'un point se fait avec une complexité en temps $O(\log {N})$ au lieu de $O(N)$, $N$ étant le nombre maximum d'itérations.

## Un zoom commenté

L'ensemble de Mandelbrot doit beaucoup sa popularité à la variété et la beauté de ses structures et à la profondeur infinie de ses détails, mais aussi à la possibilité de l'explorer soi-même à l'aide des nombreux logiciels aujourd'hui disponibles.
La séquence d'exploration commentée ci-dessous est un zoom profond vers la valeur de c = –0,743643887037151 + 0,13182590420533i, à travers nombre de motifs caractéristiques. Le rapport de grossissement entre la dernière et la première image est d'environ 60 milliards (si la dernière image était en taille réelle, la première ferait une petite dizaine de fois la distance Terre-Lune). La séquence entière (et sa suite jusqu'à un grossissement d'environ 10$^{30}$) peut être visionnée sur l'animation ci-contre.
* Description
* L'ensemble de Mandelbrot initial.
* Nous allons zoomer sur la vallée située entre la cardioïde et le bourgeon principal.
* Cette vallée a été baptisée « vallée des hippocampes ».
* À gauche, des spirales doubles, à droite les « hippocampes ».
* Nous zoomons sur l'un d'eux.
* Un « hippocampe », tête en bas. Cet hippocampe est composé de 25 « antennes » consistant en 2 groupes de 12 et un filament relié à la cardioïde. Nous en déduisons que le bourgeon qui le porte a une périodicité de 25. Le point de rencontre de ces antennes est un « point de Misiurewicz ». Sur la plus longue antenne, celle qui mène à la « queue » de l'hippocampe, on reconnaît une copie réduite de l'ensemble de Mandelbrot, appelée aussi « satellite ».
* Nous zoomons maintenant sur la queue de l'hippocampe.
* L'extrémité de la queue, enroulée en spirale, est aussi un point de Misiurewicz.
* Nous zoomons sur le haut de l'image.
* Une section de la queue. Cette structure complexe est composée d'un seul et unique chemin qui mène jusqu'à l'extrémité de la queue. L'ensemble de Mandelbrot est un ensemble simplement connexe, ce qui signifie qu'il n'y a ni boucles, ni ilots.
* Zoomons vers le centre de l'image.
* Un deuxième « satellite » (ou « minibrot ») apparait, au cœur de ce carrefour. Un exemple d'autosimilarité : La frontière de l'ensemble de Mandelbrot contient une infinité de copies de lui-même. Quel que soit le lieu on l'on zoome on en trouvera toujours au moins un. Les deux spirales sont le début d'une série de couronnes concentriques, avec le satellite en son centre.
* Nous zoomons vers ce satellite.
* Chacune de ces couronnes est composée de spirales similaires. Leur nombre s'accroit en puissances de 2, un phénomène typique de l'environnement des satellites. Le chemin vers l'extrémité de la queue entre dans le satellite par le point d'inflexion de la cardioïde et ressort par l'extrémité son antenne.
* L'antenne du satellite. On peut distinguer plusieurs satellites du deuxième ordre.
* Zoomons en haut à droite de l'image.
* La « vallée des hippocampes » du satellite. Toutes les structures déjà rencontrées précédemment réapparaissent.
* Spirales doubles et hippocampes. Contrairement à la première vallée, celle-ci est peuplée, en plus, de légères structures spirales. Dans un satellite d'ordre n cohabitent n+1 types de structures différentes. Pour ce satellite d'ordre 1 coexistent donc deux types de structures différentes.
* Nous zoomons sur une double spirale, à gauche de la vallée.
* Double-spirale avec satellites du deuxième ordre.
* Nous zoomons, maintenant, sur une des structures dentelées blanches en haut à droite de l'image.
* Ces structures légères rappellent celles de certains ensembles de Julia. Là encore, ce motif très torturé n'est constitué que d'un seul filament.
* Nous zoomons sur le « double-crochet » que l'on distingue à droite de l'image.
* Ce double crochet rappelle, encore une fois, la forme en spirale de la queue d'un hippocampe.
* Nous zoomons vers le motif central.
* Des îlots apparaissent, à la manière d'une poussière de Cantor. La forme générale est celle d'un ensemble de Julia J$_c$. Toutefois, contrairement à un ensemble de Julia, ces points sont tous connectés, car nous sommes toujours dans l'ensemble de Mandelbrot. La forme générale de cet ensemble n'est pas celle de l'ensemble de Julia associé à cette position. Il s'agit de l'ensemble de Julia que nous aurions obtenu si nous avions sélectionné, au début de notre exploration, un point à proximité d'une double spirale au lieu d'un hippocampe.
* Ces structures sont elles-mêmes connectées à une structure centrale que nous pourrions découvrir si nous poussions encore plus loin le grossissement. En théorie, le grossissement pourrait ainsi être infini, et faire découvrir sans cesse de nouvelles structures. Pour en voir davantage, cliquer sur l'animation en haut à droite de cette section.

## Généralisations et variantes

L'ensemble de Mandelbrot peut être généralisé pour des puissances d réelles supérieures à 2 pour z ↦ z$^d$ + c. Ces généralisations sont parfois appelées « multibrot », mais pour certains auteurs (comme McMullen), le terme « ensemble de Mandelbrot » doit également désigner ces généralisations.
* Puissance 3
* Puissance 4
* Puissance 5
* Puissance 6
Une modification de la méthode de traçage, proposée par Melinda Green en 1993, mène à une variante baptisée « Buddhabrot ». On y voit la densité des points visités par les orbites correspondant aux valeurs de c qui divergent, donc choisies à l'extérieur de l'ensemble de Mandelbrot.
En 3 dimensions, il n'existe pas de structure de corps comparable à celle des nombres complexes, donc pas d'extension « naturelle ». Notons toutefois l'extension de Daniel White, réalisée en 2009, baptisée « Mandelbulb. »
En 4 dimensions, l'extension naturelle au corps des quaternions a été étudiée par John Holbrook en 1987. D'autres résultats théoriques ont été traités pour l'ensemble de Mandelbrot quaternionique$^,$ (voir vidéo) défini avec un polynôme quadratique de la forme $q^{2}+c$.
En 2000, Dominic Rochon a utilisé l’anneau commutatif des nombres bicomplexes pour donner une nouvelle version de l’ensemble de Mandelbrot en dimension trois et quatre. De ses travaux émerge le « Tétrabrot »$^,$.
Il est aussi possible de généraliser encore, et de considérer plutôt des itérations quelconques $z_{n+1}=f(z_{n},c)$.
* Une variante pour f(z,c)=cos(z/c)
* f(z,c)=sinh(z)+1/c²
* f(z,c)=cos(z)+1/c
* f(z,c)=c$^{z-1}$e$^{-c}$ (intégrande de la fonction Gamma)
* f(z,c)=sinh(z)+1/c$^k$ (k variant entre 0.1 et 10)
* f(z,c)=z²+sin(c$^3$)
* f(z,c)=z²+c$^3$-δ, avec δ≈-1.401155 le point de Feigenbaum
* f(z,c)=z²+(c$^3$+0.7-0.2i)$^3$+i (zoom en 0.36+0.7i)

## Références culturelles

L'ensemble de Mandelbrot est l'un des objets mathématiques les plus populaires et les plus connus. Il est évoqué par divers artistes, dans des livres ou des musiques notamment :
* La musique de Jonathan Coulton Mandelbrot Set est un hommage à l'ensemble de Mandelbrot lui-même et à son « père », Benoît Mandelbrot. Les paroles rappellent entre autres la définition par récurrence de l'ensemble$^,$$^,$.
* Le second livre de la série Mode (livres) (en) de Piers Anthony, "Fractal Mode", décrit un monde qui est un modèle 3D parfait de l'ensemble.
* Le roman de Arthur C. Clarke "Le Fantôme venu des profondeurs" contient un lac artificiel créé pour répliquer la forme de l'ensemble de Mandelbrot.
* Le chanteur de heavy metal sud-coréen Norazo (en)a fait un clip musical Ni pal za ya (ton destin), qui commence avec une vidéo hypnotique de l'ensemble de Mandelbrot.
* L'album Jupiter's Darling du groupe de rock américain Heart présente un ensemble de Mandelbrot sur sa couverture. L'ensemble est tourné de manière que le cardioïde soit en haut et que la figure ressemble à un cœur.

## Logiciels générateurs de fractales

* ChaosPro - (propriétaire - pour Windows)
* FractaNep (gratuit - pour Windows)
* Fractile Plus (gratuit pour iOS)
* Fractint (propriétaire - porté sur de nombreuses plates-formes)
* Gecif
* Gnofract 4D - (puissant et rapide - pour Linux, FreeBSD ou Mac OS X)
* Kalles Fraktaler (gratuit - pour Windows)
* Mandelbulb 3D (gratuit pour de nombreuses plateformes, exemple de réalisation : Vidéo)
* XaoS (Logiciel libre - Licence GPL, pour Windows, Mac OS X, Linux et systèmes type Unix)