Runge-Kutta

レポートで、微惑星集積のプログラムを組むことに。



dm/dt = nmvR^2(1+2θ)　(θ: Safronov 数)



学部２年の演習以来２年半ぶりに、ルンゲクッタで書いてみようと思い立ったが、３０分で断念。面倒。

だいたい、精度いらないしね。



計算すると、地球サイズの惑星に成長するまで約 5000 億年！



うーん。

パラメータの見直しはまた今度。