調和振動子の運動方程式は以下の通りです。ここで、pは運動量、qは位置、mは質量、kはバネ係数です。 これは解析的に解くことができます。 初期条件をp(0)=0, q(0)=1として、Maximaに解かせてみましょう。 assume(m > 0, k > 0)$ atvalue(q(t), t = 0, 1)$ a…