2011-08-21から1日間の記事一覧

微分方程式の数値解析 #6

落体の運動方程式をRK4により数値的に解くFORTRANプログラム。 Euler法によるプログラムのうちy(t)とv(t)からy(t+h)とv(t+h)を求める漸化式部分をRK4に変更しただけ。 同じ計算をnextyとnextvで行っているのでおいしくないプログラムになっている。 もっとま…