How do I update my integrated function array correctly?

Question

I have been trying to use the Runge-Kutta45 integration method to update a set of positions and velocities of particles in space to get the new state at some time step.

I