#pragma once

#include <tuple>
#include <vector>

std::tuple< std::vector<double>, std::vector<double>, std::vector<T> >
  adaptive_dormand_prince( T f( double t, T y ),
                           double t0, T y0, double tf,
                           double h, double min_h,
                           double eps_abs, std::size_t max_iterations );

#include "dormand_prince.tpp"