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