#pragma once

#include <tuple>
#include <vector>

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

#include "euler_heun.tpp"