Initialize and use the propagator of the multisystem class
This should be relatively simple to implement, but there is one question that is not so clear to me: what time step should this propagator use? The smallest time step of the subsystems or the largest one? Probably the smallest one makes more sense.