下の微分方程式の近似解を rk23() で求める。 1 階の連立に直してから解く。 θ' = vel vel' = θ'' = (-g/l)*sin(θ) {0, 2} は t の振り幅。 {1, 0} は、1 が θ の初期値、0 が vel の初期値。 0.5 は t の刻み幅。 1.E-10 は error tolerance。 出力された行…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。