二重ばね振り子の運動方程式をオイラーラグランジュ方程式から導出し、それを基に数値計算を行います。まずラグランジアンを計算します。図を用意するのが面倒だったので、二重振り子の場合と同様の変数設定をしていると考えてください。ただし質点はばねで…

結果 コード 太陽、地球、月と書いてありますが上の結果は太陽でも地球でも月でもありません。 初期条件と書かれたstateの部分には太陽、地球、月のx, y座標およびx, y方向の初期速度が入ります。 import numpy as np import matplotlib.pyplot as plt from …