Kuramoto-modelをExcelで解いてみる。

Kuramoto-Sivashinskyに続いて蔵本モデルを解いてみよう。時間発展にはRunge-Kutta5次のCash-Karpを使おう（8次Dormand-Princeはちょっと時間かかりすぎる）。

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

K=2の時の結果はこれ。

だんだん振動がそろってくるのがわかる。

この本を参考にしました。これでホタルの集団発光の話を初めて知った。でももう絶版みたいだ。