yumulog | 社会人博士の日記

はてなブログに引っ越しました / 石川県に引っ越しました

Runge-Kutta

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



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



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

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



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



うーん。

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