功能推导指南

我有一个用于积分微分方程系统的代码。它是一个带有参数的模板函数——系统的维度，类型为int。它看起来像这样：

template<int n> std::array<double, n> operator+(const std::array<double, n>& a, 
    const std::array<double, n>& b)
{       // Вспомогательная функция
    std::array<double, n> ret;
    for (int i = 0; i < n; i++) ret[i] = a[i] + b[i];
    return ret; // NRVO
}

template<int n> std::array<double, n> operator*(double h,const std::array<double, n>& b)
{        // Вспомогательная функция
    std::array<double, n> ret;
    for (int i = 0; i < n; i++) ret[i] = h * b[i];
    return ret; // NRVO
}

template<int n> std::vector<std::array<double, n>>
    rk4(std::function<std::array<double, n>(double, const std::array<double, n>&)> f,
        double t0, const std::array<double, n>& x0, double tend, double h)
{        // Функция для решения системы
    if (tend < t0) return {};
    size_t points = static_cast<size_t>(1 + floor((tend - t0) / h));
    std::vector<std::array<double, n>> ret{ points };
    ret[0] = x0;
    double t = 0;
    for (size_t i = 1; i < points; i++)
    {       // Метод Рунге-Кутты 4 порядка с пост. шагом
        auto k1 = f(t, ret[i-1]);
        auto k2 = f(t + h / 2, ret[i - 1] + (h / 2 ) * k1);
        auto k3 = f(t + h / 2, ret[i - 1] + (h / 2 ) * k2);
        auto k4 = f(t + h, ret[i - 1] + h * k3);
        ret[i] = ret[i - 1] + (h / 6) * (k1 + 2*k2 + 2*k3 + k4);
    }
    return ret;
}

这是我的使用方法：

std::array<double, 2> test1(double t, const std::array<double, 2>& x)
{    // функция, считающая правую часть системы размерности 2
     return { x[1], -x[0] };  
}

代码中的某处：

auto result = rk4<2>(test1, 0, { 1,0 }, 6.28*10, 0.05);

您不能删除尖括号中的两个，编译器将无法理解什么n等于两个。尽管这里显然是多余的，因为只有 2n才能相同。我试图为rk4演绎指南编写，以便从第三个参数的大小推断维度：

template <typename F,typename T0, typename X0, typename Tend, typename H> 
    rk4(F, T0, X0, Tend, H)->rk4<std::tuple_size<X0>::value>;

但它不起作用。我的问题是：是否可以在这里写一个演绎指南，如果可以，怎么写？