2015-06-13 (29)二重振子 1(rk23() が適用できる形に運動方程式を変形する) 振子の運動 運動方程式は下のとおりである。m は質量、v は角速度、a は角加速度である。(1) Nspire の組込函数 solve() を使って加速度 a について解く。(2) 下のように解ける。(3) 今度は θ' = v θ'' = v' = a と置き換えて一階の連立に直す。(4) これで rk23() が適用できる形になった。(5) 行列が返ってくる。1 行目が時間 t、2 ~ 5 行目がそれぞれ θ1、v1、θ2、v2。 参考文献:工業力学 (機械工学基礎講座)作者:敏博, 入江,元, 山田発売日: 1980/09/01メディア: 単行本(pp. 229 - 231)