Cette fractale a été découverte par Michel Hénon, de l'observatoire de Nice (France), en 1976, au cours de recherches sur les orbites d'objets astronomiques. |
x' = 1 + y - ax y' = bx |
Longtemps dans le passé (au siècle dernier), après avoir implémenté la formule ci-dessus, je me suis demandé ce qu'il se passait si on faisait varier la constante. En gros, nous allons balayer le plan x-y avec a et b, et en tirer les conclusions qui s'imposent. J'aime bien faire varier les constantes, on tombe parfois sur de petits joyaux tels que celui-ci:
Si j'ai bien retrouvé le code fortran77 qui calcule cette image, je n'ai pas retrouvé les explications qui vont avec. Il va donc falloir que je les reconstitues à partir du source... Pour les courageux, c'est map_henon.f dans ce tarball. Un de ces jours, je ferais le nettoyage dans tout ce fatras...