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

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

Euler法よりもマシなはずの中点法で計算してみる。 Euler法との違いは漸化式のの定義部分だけである。 ...snip double precision function nextx(t, x, h) double precision, intent(in) :: t, x, h double precision :: k1, k2 k1 = f(t, x) k2 = f(t + h /…