#pragma once

#include <tuple>
#include <vector>

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

#include "fehlberg.tpp"