2008-03-13から1日間の記事一覧

ルンゲクッタ法で1階微分方程式を解く

微分方程式を解くルンゲクッタ法を試してみます。 今回は という微分方程式を解いて、空気抵抗がある自由落下の速度を求めてみます。g=9.8、k=0.01としてます。 空気抵抗により、速度が一定になることがわかります。 プログラム中、funcメソッドが求める微分…