namespace boost {
  namespace numeric {
    namespace odeint {
      template<typename Stepper, typename System, typename State, 
               typename TimeIterator, typename Time, typename Observer> 
        size_t integrate_times(Stepper, System, State &, TimeIterator, 
                               TimeIterator, Time, Observer);

      // Solves the forwarding problem, can be called with Boost.Range as start_state. 
      template<typename Stepper, typename System, typename State, 
               typename TimeIterator, typename Time, typename Observer> 
        size_t integrate_times(Stepper stepper, System system, 
                               const State & start_state, 
                               TimeIterator times_start, 
                               TimeIterator times_end, Time dt, 
                               Observer observer);

      // The same function as above, but without observer calls. 
      template<typename Stepper, typename System, typename State, 
               typename TimeRange, typename Time, typename Observer> 
        size_t integrate_times(Stepper stepper, System system, 
                               State & start_state, const TimeRange & times, 
                               Time dt, Observer observer);

      // Solves the forwarding problem, can be called with Boost.Range as start_state. 
      template<typename Stepper, typename System, typename State, 
               typename TimeRange, typename Time, typename Observer> 
        size_t integrate_times(Stepper stepper, System system, 
                               const State & start_state, 
                               const TimeRange & times, Time dt, 
                               Observer observer);
    }
  }
}