#pragma once

#include <tuple>
#include <vector>

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

#include "dormand_prince.tpp"