時間遅れのあるロジスティック方程式をExcel VBAで計算。

遅延微分方程式を計算するシリーズ。”常微分方程式の数値解法I”に載っていたルンゲクッタを基にした方法で計算している。そこにも載っていた人口力学の例、時間遅れのあるロジスティック方程式（ライトの方程式とも）はこんな感じの式。

dy(t)/dt = (a - y(t-1) )*y(t)

でaをパラメータにした計算結果がこちら。

振動してますな。