Дано разложение функции y (x) в ряд Тейлора. Составить программу, вычисляющую для значения x и точности Eps

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cmath>
using namespace std;
 
double Y(double x)
{
    return cosh(x);
}
 
double S(double x, double e)
{
    double ak = 1, sum = 0;
    long k = 0;
    while (e < fabs(ak))
    {
        sum += ak;
        ak *= pow(x, 1) / (1 * k + 1);
        k++;
    }
    return sum;
}
 
double D1 (double y, double ak, double d1)
{
    d1 = y-ak;
    return d1;
}
 
double D2(double d2, double d1, double y)
{
    d2 = d1 / y;
    return d2;
}
 
 
int main()
{
    double e, x, d1;
    cout << "Enter x : "; cin >> x;
    cout << "Enter e : "; cin >> e;
    cout << "                     Rezultati vichislenui " << "\n";
    cout << "Tochnoe znachenie funkcii Y(x) = " << Y(x) << endl;
 
 
 
 
 
    cout << "-------------------------------------------------\n";
    cout << "|  N  |  Slagaemoe |   S(x)   |    d1    |    d2    |\n";
    cout << "-------------------------------------------------\n";
    cout << "|  " << 1 << "   " << x << "   " << pow(x, 1) / (1 * 1 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 1) / (1 * 1 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 1) / (1 * 1 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "|  " << 2 << "   " << pow(x, 3) / 3 * 2 * 1 << "   " << pow(x, 3) / (1 * 3 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 3) / (1 * 3 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 3) / (1 * 3 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "|  " << 3 << "   " << pow(x, 5) / 5 * 4 * 3 * 2 * 1 << "   " << pow(x, 5) / (1 * 5 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 5) / (1 * 5 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 5) / (1 * 5 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "|  " << 4 << "   " << pow(x, 7) / 7 * 6 * 5 * 4 * 3 * 2 * 1 << "   " << pow(x, 7) / (1 * 7 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 7) / (1 * 7 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 7) / (1 * 7 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "|  " << 5 << "   " << pow(x, 9) / 9 * 8 * 7 * 6 * 5 * 4 * 3 * 2 * 1 << "   " << pow(x, 9) / (1 * 9 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 9) / (1 * 9 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 9) / (1 * 9 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "|  " << 6 << "   " << pow(x, 11) / 11 * 10 * 9 * 8 * 7 * 6 * 5 * 4 * 3 * 2 * 1 << "   " << pow(x, 11) / (1 * 11 + 1) << "  " << -1 * ((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 11) / (1 * 11 + 1) << "  " << -1 * (((pow(e, x)*pow(e, x) - 1) / (2 * pow(e, x))) - pow(x, 11) / (1 * 11 + 1) / pow(e, x)*pow(e, x) - 1) / 2 * pow(e, x) << " |""\n";
    cout << "--------------------------------------------------------------------------\n";
   
    system("pause");
    return 0;
}

double S(double x, double e)
{
    double ak,sum;
    ak=x;
    sum=ak;
    long k = 2;
    while (e < fabs(ak))
    {
        ak*=x*x/(2*k-2)/(2*k-1);
        sum += ak;
        k++;
    }
    return sum;
}