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

#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
 
int main()
{
    const double ESP = 0.0001;
    int pow_num = 0;
    int lit = 0;
    int n = 0;
    double x = -3;
    double result = 1;
    double sum = 0;
    double sin;
 
    cout << "X      |Sum         |Function         |Literals" << endl;
    cout << "-----------------------------------------------" << endl;
    do
    {
        do
        {
 
             while (pow_num != 2 * n);
            {
                if (x == 0)
                {
                    result = 0;
                    pow_num = 2 * n;
                }
                else
                {
                    result *= x;
                    pow_num++;
 
                }
            }
            sum += result /(2n-1)!; //  - здесь нужно поделить результат на факториал (2n-1)!
            sin = sinh(x) / x;
            pow_num = 0;
            result = 1;
            n++;
            lit++;
        }  while (sum - sin > ESP || sin - sum > ESP);
        cout << fixed;
        cout << setprecision(2) << x << "\t|";
        cout << setprecision(6) << sum << "\t|";
        cout << sinh(x) / x << "\t|";
        cout << lit << endl;
        sum = 0;
        lit = 0;
        n = 0;
        x += 0.05;
    } while (x < 3);
      cout << "-----------------------------------------------" << endl;
 
    return 0;
 
}

#include <iostream>
#include <iomanip> 
#include <cmath>
using namespace std;
 
int main()
{ 
    int n;    
    double xstart=-3., xend=3., h=0.05, an, S, Y, eps=0.0001; 
    
    cout <<"   x     S(x)     Y(x)\n";
    for (double x=xstart; x<=xend+h/2.; x+=h)
    {   
        n=0;
        an=1.; 
        S=0.;
        while(fabs(an)>eps)    
        { 
            S+=an;
            n++;
            an*=x*x/(2*n*(2*n+1));            
        }    
    Y=sinh(x)/x;
    cout << setprecision(2) << fixed<< x << "  " << setprecision(5) << S << "  " << Y << "\n";
    }   
          
system("pause");
return 0;
}