Fractales autour de Hénon

{diagramme de Hénon}

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

Le mapping de tTh-Hénon

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:

{mapping de Hénon}
La même, mais en plus grand (1Mo)

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...