Вывести таблицу значений функции

/*
Для х , которое изменяется от а до b с шагом ( b - а ) / ( k - 1),
где k - количество точек, вычислить функцию у (х), используя ее
расписание в ступенчатый ряд для двух случаев : для заданного п
и для заданной точности ε = 0.0001 . Для сравнения найти точное
значение функции .
*/
 
#include <cstdlib>
#include <cmath>
#include <iostream>
 
unsigned int factorial(unsigned int n)
{
    return (n <= 1) ? 1 : n * factorial(n - 1);
}
 
double func_by_n(double x, unsigned int n)
{
    return (n == 0) ? 1 : pow(x, 2.0 * n) / factorial(2 * n) + func_by_n(x, n - 1);
}
 
double func_by_e(double x, double e)
{
    double res = 0.0;
    double prev_res = 0.0;
    unsigned int n = 0;
    do {
        prev_res = res;
        res += pow(x, 2.0 * n) / factorial(2 * n);
        ++n;
    } while (abs(res - prev_res) > e);
    return res;
}
 
double func(double x)
{
    return (exp(x) + exp(-x)) / 2.0;
}
 
int main()
{
    setlocale(0, "rus");
 
    double x0 = 0.0, x1 = 1.0;
    double x = x0;
    unsigned int k = 11;
    double h = (x1 - x0) / (k - 1);
 
    std::cout << "x\ty(x, n)\ty(x, e)\ty(x)\n\n";
    do {
        std::cout << x << '\t' << func_by_n(x, 5) << '\t' << func_by_e(x, 0.0001) << '\t' << func(x) << '\n';
        x = x + h;
    } while (!(x > x1));
 
    system("pause");
    return 0;
}