お題:二重振り子を数値計算でシミュレーションせよ