pub fn direct_gravity_step(sim: &mut NBodySimulation)
Direct O(n²) gravity calculation and Euler step.