- The program uses predictor-corrector methods to advance the planets’ positions. The forces come from Newton’s law of gravitation and F=ma. The new position is first guessed using the present position, velocity, and acceleration only. This “prediction” is then “corrected” by incorporating the knowledge of the force at the newly guessed position. Without this correction step, the numerical scheme is unstable. This scheme is really simple.