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

#include <iostream>
#include <cmath>
using namespace std;
 
double fact(const int n) { 
    return n <= 0 ? 1.0 : n * fact(n - 1);
}
 
int main() {
    double k, sum, x, f, r = -1.0, eps;
    cout << "Enter epsilon value: ";
    cin >> eps;
    for (x = 1; x <= 15; x += 1) {
        k = 0.0;
        sum = 0.0;
        f = pow(-1.0, k) * pow(x, 2 * k) / fact(k);
        sum = sum + f;
        while (fabs(f) > eps) {
            k++;
            f = pow(-1.0, k) * pow(x, 2 * k) / fact(k);
            sum += f;
        }
        cout << "sum = " << sum << endl;
        cout << "exp(" << -pow(x, 2) << ") = " << exp(-pow(x, 2)) << endl;
    }
    return 0;
}

#include <iostream>
#include <cmath>
int main()
{
    double eps;
    std::cout << "Enter epsilon value: "; std::cin >> eps;
    std::cout << "x\tsum\t\te^(-x^2)" << std::endl;
    for (double x = 1; x <= 15; x += 1)
    {
        double u = 1.0, sum = 1.0;
        int k = 0;
        while (u < -eps || eps < u)
        {
            k++;
            u *= -x * x / k;
            sum += u;
        }
        std::cout << x << '\t' << sum << '\t' << exp(-x * x) << std::endl;
    }
    return 0;
}