Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その15) 蔵本モデルを計算してGIFアニメに。

先日は蔵本・シバシンスキー方程式でしたが、今日は蔵本モデル。

dφi/dt = ωi + K/N Σsin (φj - φi )          （和はj=1 ～ N)

位相がどんどんそろっていくのがわかる。

クリックでGIFアニメが始まります。