単振動の微分方程式を刻みhについてルンゲクッタで求める。

質問

m*d^2x/dt^2=-kx x(0)=1 dx/dt=0
というのが与えられて二階微分だから一階微分にするために
dx(t)/dt=v(t) dv(t)/dt=-k*x(t)/mという式を立てました。オイラー法ではできたのですが2次、4次のルンゲクッタだとできません。どなたか回答お願いします。

回答No.2

ojisan7

ひとつ補足をするのを忘れていました。
なにか、ルンゲクッタ法を使ったプログラムのサンプルに目を通すことも、使い方を学ぶ上で参考になると思います。そして、自分で、簡単なプログラムを作り、何か簡単な微分方程式で実行してみることです。

回答No.1

ojisan7

ルンゲクッタ法は２次にしても、４次にしてもステップ数はそんなに多くはないので（数ステップ程度）簡単にできると思います。教科書にある公式通りに入力すればできると思います。がんばってください。