Henon-HeilesをExcelで解いてみる（シンプレクティック法）。

三体問題、FPU問題に続いてHenon-Heilesの微分方程式をシンプレクティック法で解いてみた。ポテンシャルは

U(x, y) = 1/2 (x² + y² + 2 x²y - 2/3 y³)

で、E=0.125になるようにいろんな初期値で計算した結果10条件を重ね合わせた結果がこれ。

前からこの図は人の顔（叫びみたいな）を横向けたように見えると思ってたんだが、こうみるとバケタンにも似てるな（https://sci.tea-nifty.com/blog/2008/06/post_535a.html）。