For inspiration, see ´Precise_orbit(path,t_0)´ in https://github.com/lupemba/syntese/blob/master/ToolBox/src/Load.jl
-- I think we should drop the t_0 input.
-- Make a ´struct´ named OrbitState with the fields; time, position, and velocity.
-- Make the function return an array of OrbitState, ´Array{OrbitState}´