レポートで、微惑星集積のプログラムを組むことに。 dm/dt = nmvR^2(1+2θ) (θ: Safronov 数) 学部２年の演習以来２年半ぶりに、ルンゲクッタで書いてみようと思い立ったが、３０分で断念。面倒。 だいたい、精度いらないしね。 計算すると、地球サイズの惑星…