# Ensemble de Multibrot

En mathématiques, l'ensemble de multibrot est 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}^{d}+c\end{cases}}$
où d ≥ 2,
est bornée$^,$$^,$.
L'exposant d peut être généralisé à des valeurs négatives et fractionnelles.
Ce nom est la contraction de multiple et ensemble de Mandelbrot (ce dernier correspond au cas particulier où d=2).
Beaucoup de représentations graphiques sont disponibles, mais pour l'instant, aucune n'a réussi à afficher un empilement 3D de différentes images pour que l'évolution générale de la forme puisse être vue depuis un autre endroit que d'en haut.

## Exemples

Le cas de
    $d=2\,$
est l'ensemble de Mandelbrot classique depuis lequel le nom est dérivé.
Les ensembles pour les autres valeurs de d présentent aussi des images fractales lorsqu'ils sont tracés sur le plan complexe.
Chaque exemple des différentes puissances d ci-dessous est tracé à la même échelle. Les valeurs de c appartenant à l'ensemble sont noires. Les valeurs de c dont les valeurs ne sont pas bornées après itérations n'appartiennent pas à l'ensemble et sont représentées en différentes couleurs pour montrer leurs contours. Les couleurs dépendent du nombre d'itérations nécessaires pour que la valeur dépasse la magnitude définie.

### Puissances positives

L'exemple d = 2 est l'ensemble de Mandelbrot original. Les exemples pour d > 2 sont souvent appelés ensembles de multibrot. Ces ensembles incluent l'origine et ont un périmètre fractal, avec (d − 1) axes de symétrie rationnel(s).
| z $\mapsto$ z$^2$ + c | z $\mapsto$ z$^3$ + c | z $\mapsto$ z$^4$ + c |
| z $\mapsto$ z$^5$ + c | z $\mapsto$ z$^6$ + c | z $\mapsto$ z$^{96}$ + c |
* z $\mapsto$ z$^{96}$ + c détail x40

### Puissances négatives

Lorsque d est négatif, l'ensemble encercle l'origine mais ne l'inclut pas. Il y a un comportement complexe entre l'ensemble et l'origine, une zone en forme d'étoile avec (1 − d) axes de symétrie rationnel(s). L'ensemble semble avoir un périmètre circulaire, mais cette illusion est créée par le rayon maximal autorisé par l'algorithme et n'est pas la vraie limite de l'ensemble, qui s'étend en réalité dans toutes les directions jusqu'à l'infini.
| z $\mapsto$ z$^{−2}$ + c | z $\mapsto$ z$^{−3}$ + c | z $\mapsto$ z$^{−4}$ + c | z $\mapsto$ z$^{−5}$ + c | z $\mapsto$ z$^{−6}$ + c |

### Puissances fractionnelles

L'inclusion de puissances fractionnelles ne montre rien de particulier, si ce n'est l'apparition des axes de symétrie que l'on observe avec des puissances entières. La transition semble régulière, le nouvel axe "apparait" sans causer d'irrégularités par rapport aux formes précédemment obtenues.

## Tracer l'ensemble

Toutes les images ci-dessus sont tracées avec l'algorithme d'Évasion du Temps qui identifie des points hors de l'ensemble de manière simple. On peut révéler beaucoup plus de détails fractals en traçant l'exposant de Liapounov comme le montre l'exemple ci-dessous.
L'exposant de Liapounov est le taux de croissance de l'erreur d'une séquence donnée. On calcule d'abord la séquence d'itération avec N itérations, puis on calcule l'exposant avec
    $\lambda ={\frac {1}{N}}\ln |\mathbf {z} |$
et si l'exposant est négatif, la séquence est stable. Les pixels blancs sur l'image sont les paramètres c pour lesquels l'exposant est positif, c'est-à-dire instable.
Les couleurs montrent la période des cycles auxquels les orbites sont attirées. Tous les points bleu foncé (extérieur) sont attirés par un point fixe et tous les points au milieu, en bleu plus clair, sont attirés par un cycle à deux périodes, et ainsi de suite.

### Pseudocode

La valeur complexe z a les coordonnées (x, y) sur le plan complexe et est élevée à des puissances différents dans la boucle d'itération par les codes montrés ci-dessus. Les puissances non présentes dans ce code sont obtenues en concaténant des codes montrés.
| z$^{−1}$ | z$^2$ (pour l'ensemble de Mandelbrot) | z$^3$ | z$^5$ | z$^n$ |
| d=x^2+y^2 if d=0 then ESCAPE x = x/d + a y= -y/d + b | xtmp=x^2-y^2 + a y=2*x*y + b x=xtmp | xtmp=x^3-3*x*y^2 + a y=3*x^2*y-y^3 + b x=xtmp | xtmp=x^5-10*x^3*y^2+5*x*y^4 + a y=5*x^4*y-10*x^2*y^3+y^5 + b x=xtmp | xtmp=(x*x+y*y)^(n/2)*cos(n*atan2(y,x)) + a y=(x*x+y*y)^(n/2)*sin(n*atan2(y,x)) + b x=xtmp |