4変数の常微分方程式をルンゲクッタ法（4段4次）で計算する自作式をカシオの高精度計算サイトkeisan.casio.jpにUP!

ちょっと必要になって4変数の常微分方程式を計算しようとしたが、お手軽にweb上で計算できるものが欲しいな、と思ったのでカシオの高精度計算サイトに自作式として作ってみた。

以前は3変数のものをUPしていたが、あまり需要がなかったので消してた。

今度は4変数として復活。あんまりweb版では見かけないパターンじゃないかと思って（Wolfram Alphaなら何でも計算できますが、、、）

リンクはこちら：

4変数の常微分方程式をルンゲクッタ法（4段4次）で計算

画面イメージはこんな感じ。

説明文：

dx/dt=f1(x,y,z,w) dy/dt=f2(x,y,z,w) dz/dt=f3(x,y,z,w) dw/dt=f4(x,y,z,w) という4変数の常微分方程式をルンゲクッタ法（4段4次）で計算します。

例題として、x,yが角度が小さいときの振り子、z,wが任意の角度の振り子を計算します。 3変数の時はf4=0、2変数の時はさらにf3=0とでもしてください。