Computation of numerical method is inherently iterative: at each time step, the solution depends on the immediately preceding one. |
At each time step, Grav is computed:
|
We will use 4th order Runge Kutta to integrate in time and the program is designed an overall routine looping over time with parallelism hidden in Grav routine |