Взятие производной любого порядка

//Обёртка для функции
    function<double(double)> f;
//Пакет точек функции
    vector<pair<double, double>> Data_;
 
//Задача функции
    f = [](double x)
    {
        double y = cos(x) * cos(x) - 0.3 * cos(x) - 0.4;
        return y;
    };
 
    function_point_cloud(f, Data_, -1000., 1000., 2500000);
    wcout << get_derivative_n_in_point(Data_, 2., 4u);
    Data_.clear();

//Облако точек функции
void function_point_cloud
(
    const function<double(double)> &f,
    vector<pair<double, double>> &v,
    const double &min,
    const double &max,
    const double &count_points_split
)
{
    double step = (max - min) / count_points_split;
    double left = min;
    for (size_t u = 0u; u <= count_points_split; ++u, left += step)
    {
        v.push_back(make_pair(left, f(left)));
    }
}

//Взятие любой производной
double get_derivative_n_in_point(const vector<pair<double, double>> &points, const double &x, const size_t &n)
{
    vector<pair<double, double>> box_old(points);
    vector<pair<double, double>> box_new;
    for (size_t p = 0u; p < n; ++p)
    {
        box_new = vector<pair<double, double>>(box_old.size() - 2);
        for (long long u = 1u; u < box_old.size() - 1; ++u)
        {
            double y = (box_old[u + 1].second - box_old[u - 1].second) / (box_old[u + 1].first - box_old[u - 1].first);
            double x = box_old[u].first;
            box_new[u - 1] = make_pair(x, y);
        }
        box_old = box_new;
    }
    size_t ind = 0u;
    for (const auto &value : box_new)
    {
        if (x < value.first)
        {
            break;
        }
        ++ind;
    }
    wcout << ind << endl;
    double res = abs(box_new[ind].first - x) < abs(box_new[ind - 1].first - x) ? box_new[ind].second : box_new[ind - 1].second;
    return res;
}