連続系課題。重力２体問題を単純化した常微分方程式を、刻み幅制御を行う線形多段階法によって解けという課題。
id:yoiraに助けてもらいながら実装。多項式の加算、乗算を使ってNewton補間をする必要がありました。多項式を実現するためにはやっぱりリストを使えるような言語の方がいいなーと思い、練習も兼ねてPythonで書いてみました。yoira君はSchemeで実装していたみたいです。
授業ではScilabという環境を薦められていたのですが、さすがにこの言語を覚える元気は無かったです。
全然バグが取れなくて悩んでいたら、そもそも積分の仕方が間違っていた...。もう微分積分とか久しく触れてないのでヤバイ。高校生に余裕で負ける自信があります。
comments