この英文の和訳をお願いします。
Before a detailed description of our numerical procedures, we comment on above simplifications. As seen in an example of e=1.0 and i=0.5, n-recurrent (n≧3) collision orbits contribute by only 1% or less to the collision rate. On the one hand, from calculations with other e and i, the degree of contribution by n-recurrent (n≧3) collision orbits is found to be largest in the case e≒1. Hence, the error in <P(e, i)> introduced by simplification (i) is of the order of 1% or less. As for the applicability of the two-body approximation, we have confirmed in PaperII that the orbit are well described by the two-body formula inside the two-body sphere, whose radius is given by Eq. (13). No appreciable error in <P(e, i)> comes from simplification (ii). Simplification (iii) follows the discussion in the last section.
Using above simplifications, we have developed numerical procedures for obtaining <P(e, i)> efficiently; their flow chart is illustrated in Fig. 10. Choosing initial values of orbital elements (e, i, b, τ_s, ω_s), we start to compute numerically Hill’s equations (6) by an ordinary fourth-order Runge-Kutta method from a starting point given by Eqs. (21) and (22). The distance r is checked at every time step of the numerical integration. If the particle flies off to a sufficient distance from the protoplanet after approaching it, i.e., if
|y|>y_0+2e, (26)
then, the orbital computation is stopped. If a particle approaches the protoplanet and crosses the two-body sphere surface, i.e., if r≦r_cr, the two-body formula is employed to predict whether or not a collision occurs. When no collision occurs at the first encounter, the numerical integration of Hill’s equations is continued. Since a particle which enters the two-body sphere inevitably escapes from the sphere (see Paper II), the particle follows alternatives: one is that it departs to such a distance that Eq. (26) is satisfied, and the other is that it crosses the two-body sphere surface again. In the former case, we stop the computation, considering that the orbit is non-collisional. In the latter case, the occurrence of collision is checked in the same way as earlier by means of the two-body formula, and the orbital calculation is terminated.
Using the numerical procedures developed in this way, we obtain <P(e, i)> in many sets of (e, i); the results are presented in the preceeding sections.
Fig. 10. Flow chart of orbital calculation for finding collision orbits.
Fig. 10. 拡大画像↓
http://www.fastpic.jp/images.php?file=2113240192.jpg
長文になりますが、よろしくお願いします。
お礼
ありがとうございます。 (#^.^#)