code: def _integrate_odeint(self, *args, **kwargs):
          
          import pyodeint  # Python interface to boost's odeint integrators
          kwargs['with_jacobian'] = kwargs.get(
              'method', 'rosenbrock4') in pyodeint.requires_jac
          return self._integrate(pyodeint.integrate_adaptive,
                                 pyodeint.integrate_predefined,
                                 *args, **kwargs)
docs: Do not use directly (use ``integrate(..., integrator='odeint')``).
  
          Uses `Boost.Numeric.Odeint <http://www.odeint.com>`_
          (via `pyodeint <https://pypi.python.org/pypi/pyodeint>`_) to integrate
          the ODE system.