Subscribed unsubscribe Subscribe Subscribe

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

振子の運動

f:id:ti-nspire:20150423071437j:plain

 

下の微分方程式の近似解を rk23() で求める。

f:id:ti-nspire:20150428063955j:plain

 

1 階の連立に直してから解く。

        θ' = vel

        vel' = θ'' = (-g/l)*sin(θ)

f:id:ti-nspire:20150428065115j:plain

{0, 2} は t の振り幅。

{1, 0} は、1 が θ の初期値、0 が vel の初期値。

0.5 は t の刻み幅。

1.E-10 は error tolerance。

 

出力された行列は 1 行目が t、2 行目が θ、 3 行目が vel。

 

 

Remove all ads