Excelで戸田格子をシンプレクティック法で解いてみる。

前にFermi-Pasta-Ulamの問題を計算したけど、ちょっと変えたらすぐToda Latticeも計算できる。計算はシンプレクティック8次だ。

式は

dQn/dt = Pn

dPn/dt = exp(Qn-1 - Qn) - exp(Qn - Qn+1)

だ。

N=128で、周期境界条件、初期条件は正弦波として計算した結果はこれ。