Fermi-Pasta-Ulam(FPU)問題をExcel VBA-シンプレクティック8次で計算

だいぶ前にシンプレクティック8次の公式をExcel VBAに移植して8の字解を計算したりフェルミ・パスタ・ウラム問題を計算したりした。

https://sci.tea-nifty.com/blog/2009/04/8excel-vba-7680.html

PCが新しくなったので規模を大きくしようと思って忘れていたので、まずはフェルミ・パスタ・ウラム問題（FPU問題）のN数を増やそう。

使った式はこれ。

 

α=1/4、N=64（オリジナルはN=32）としてシンプレクティック8次でモード毎のエネルギーを計算した結果がこれ。

N=32より格段にシャープな波形が得られた。なるほど。

で、前から気になっていた、このエネルギーを与える波形はどんなのかを次に見てみる（続く）。