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

#include <iostream>
#include <locale.h> 
#include <stdio.h> 
#include <math.h>
using namespace std;
const double e = 2.71828;
double function(double x, double d1)
{
    double res1 = 0, res2 = x, z = 1, res3 = 0;
    do
    {
        res1 = res1 + res3;
        res3 = res2;
        res2 = ((x * x) / ((2 * z + 1) * (2 * z + 2))) * res3;
        z = z + 1;
    } while (fabs(res2 - res3) > d1);
    return res1;
}
double function2(double x, double d2)
{
    double res1 = 0, res2 = x, z = 1, res3 = 0;
    
    do
    {
        res1 = res1 + res3;
        res3 = res2;
        res2 = ((x * x) / ((2 * z + 1) * (2 * z + 2))) * res3;
        z = z + 1;
    } while (fabs(res2 - res3) > d2);
    return res1;
}
int main()
{
    setlocale(LC_ALL, "RUS");
    double x1, x2, y;
    double d1 = 0.1;
    double d2 = 0.000001;
    double f;
    int k=1;
    cout << "Начальное значение х :";
    cin >> x1;
    cout << "Конечное значение х :";
    cin >> x2;
    cout << "Шаг табуляции :";
    cin >> y; 
    cout << "x   " << "f(x)   " << "first         " << "second"<<endl ;
 
    for (double x = x1; x <= x2; x += y)
    {
        f = (pow(e, (-x)) - pow(e, x)) / 2.;
        cout << x  <<"    " << f << "    " << 0.1 << "    " << function(x, d1) << "    " << 0.000001 << "    " << function2(x, d2) << endl;
        k++;
    }
    system("pause");
    return 0;
}

#include <iostream> 
#include <iomanip> 
#include <cmath>
using namespace std;
 
int main()
{ 
    int n, k;
    double xn, xk, dx, an, S1, S2, eps1=0.1, eps2=0.000001;
    
    cout << "xn="; cin >> xn;
    cout << "xk="; cin >> xk;
    cout << "dx="; cin >> dx;
    
    cout << "   x     eps1      f1(x)    eps2   f2(x)   sh(x)\n";  
    for (double x=xn; x<=xk+dx/2.; x+=dx) 
    { 
        n=1;
        S1=0.;
        an=x;       
        while(fabs(an)>eps1)    
        { 
            S1+=an;         
            an*=x*x/(2*n*(2*n+1));
            n++;
        }
        
        k=1;
        S2=0.;
        an=x;       
        while(fabs(an)>eps2)    
        { 
            S2+=an;         
            an*=x*x/(2*k*(2*k+1));
            k++;
        }      
    
    cout << x << "    " << eps1 << "    " << S1 << "    " << eps2 << "    " << S2 << "    " << sinh(x) << "\n";
    
    }
    
system("pause");
return 0;
}