Вычислить значения функций заданных разложением в ряд Тейлора

#include <iostream>
#include <iomanip>
#include <cmath>
 
// Функция для вычисления факториала
double factorial(int n) {
    double result = 1.0; // Используем double для избежания переполнения
    for (int i = 1; i <= n; ++i) {
        result *= i;
    }
    return result;
}
 
// Функция для вычисления экспоненты
double exponential(double x, int terms) {
    double result = 1.0; // Начальное значение для e^x
    double term = 1.0;   // Текущий член ряда
    for (int n = 1; n <= terms; ++n) {
        term *= x / n;   // Вычисляем n-й член
        result += term;  // Добавляем к результату
    }
    return result;
}
 
// Функция для вычисления y(x, k)
double y(double x, int k) {
    double sum = 0.0;
    double e_neg_x = exponential(x, k); // e^(x)
    for (int n = 0; n <= k; ++n) {
        sum += (pow(2, n) * pow(e_neg_x, n)) / factorial(n);
    }
    return sum;
}
 
int main() {
    double a = 0.1;
    double b = 1.0;
    int k = 120;
    double h = (b - a) / 10;
 
    std::cout << std::fixed << std::setprecision(6);
    std::cout << "x\t\t y(x, " << k << ")\n";
    std::cout << "-------------------------\n";
 
    for (double x = a; x <= b; x += h) {
        std::cout << x << "\t " << y(x, k) << "\n";
    }
 
    return 0;
}

double func(double x, int k = 120) {
    double tmp = -calculate_exp(x); // Используем нашу реализацию exp(x)
    double memb = 1.0, sum = memb;
 
    for (int n = 1; n <= k; ++n) {
        memb *= 2 * tmp / n;
        sum += memb;
    }
 
    return sum;
}