2015-04-28から1日間の記事一覧

(23)振子の運動 4(rk23() で近似解を求める)

下の微分方程式の近似解を rk23() で求める。 1 階の連立に直してから解く。 θ' = vel vel' = θ'' = (-g/l)*sin(θ) {0, 2} は t の振り幅。 {1, 0} は、1 が θ の初期値、0 が vel の初期値。 0.5 は t の刻み幅。 1.E-10 は error tolerance。 出力された行…