Dormand-Prince(ルンゲクッタ8次)DOP853 VBA移植版でローレンツ方程式を計算

さて、ようやく計算できるようになった。まずはおなじみローレンツ方程式。

dx/dt = - σ*(x - y)

dy/dt = -y - x * z + r * x

dz/dt = x * y - b * z

を計算。一発で

が計算できた。これ使っていろいろ計算していこう（続く）。

（オリジナルのDOP853は　http://www.unige.ch/~hairer/software.html　です）

ソースコードはこちら。

「DOP853.bas」をダウンロード

「Lorenz.bas」をダウンロード

ライセンスは

「licence.txt」をダウンロード

に従います。