We describe a hybrid algorithm to calculate the formation of planets from an
initial ensemble of planetesimals. The algorithm uses a coagulation code to
treat the growth of planetesimals into oligarchs and explicit N-body
calculations to follow the evolution of oligarchs into planets. To validate the
N-body portion of the algorithm, we use a battery of tests in planetary
dynamics. Several complete calculations of terrestrial planet formation with
the hybrid code yield good agreement with previously published calculations.
These results demonstrate that the hybrid code provides an accurate treatment
of the evolution of planetesimals into planets.Comment: Astronomical Journal, accepted; 33 pages + 11 figure