Интерполяционный полином Лагранжа

#include <iostream>
#include <array>
#include <algorithm>
 
template <size_t degree>
class LagrangeInterpolation {
    std::array<double, degree + 1> _x;
    std::array<double, degree + 1> _y;
    std::array<double, degree + 1> _denominators;
    double get_product(int index, double x);
public:
    LagrangeInterpolation(double a, double b);
    double operator()(double t);
};
 
template<size_t degree>
double LagrangeInterpolation<degree>::get_product(int index, double x) {
    int i = 0;
    return std::accumulate(_x.begin(), _x.end(), 1.0, [&](double prev, double xi) {
                           return i++ != index ? prev * (xi - x) : prev;
                           });
}
 
template <size_t degree>
LagrangeInterpolation<degree>::LagrangeInterpolation(double a, double b) {
    double h = std::abs(b - a) / degree;
    int i = 0;
    std::generate(_x.begin(), _x.end(), [&](){ return a + h * i++; });
    std::transform(_x.begin(), _x.end(), _y.begin(), [](double x) { return x * (std::log(x) + 2.0 * std::cos(x)); });
    i = 0;
    std::generate(_denominators.begin(), _denominators.end(), [&]()
                  {
                      double x = _x[i];
                      return get_product(i++, x);
                  });
}
 
template <size_t degree>
double LagrangeInterpolation<degree>::operator()(double t) {
   int i = 0;
   return std::accumulate(_y.begin(), _y.end(), 0.0, [&](double prev, double y) {
                            double d = _denominators[i];
                            return prev + y * get_product(i++, t) / d;
                          });
}
 
int main()
{
    LagrangeInterpolation<5> lagrange(1.0, 2.0);
    double x = 0.0;
    std::cout << "Введите x: " << std::endl;
    std::cin >> x;
    std::cout << "Значение интерполяционного полинома Лагранжа: " << lagrange(x) << std::endl;
}