code: def _initialize_integrator(self, t0, y0, integrator, **kwargs):
          
          # set the initial condition
          self.set_initial_value(y0, t0)
  
          # select the integrator
          self.set_integrator(integrator, **kwargs)
docs: Initializes the integrator prior to integration.