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

#define _USE_MATH_DEFINES
#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
 
double foo (int, int, int)
 
int main () {
    int dx, X1, X2, eps;
    cout << "Enter values for X1, X2, dx and eps" << endl;
    cin >> X1; X2; dx; eps;
    cout << setprecision(eps);
    cout << "Table of results" << endl;
    cout << "X1" << "\t" << "Result" << endl;
    foo(X1, X2, dx);
}
 
double foo(int X1, int X2, int dx) {
    
}

#define _USE_MATH_DEFINES
#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
 
double foo(int, int, int);
 
int main () {
    int dx, X1, X2, eps;
    cout << "Enter values for X1, X2, dx and eps" << endl;
    cin >> X1; X2; dx; eps;
    cout << setprecision(eps);
    cout << "Table of results" << endl;
    cout << "X1" << "\t" << "Result" << endl;
    foo(X1, X2, dx);
}
 
double foo(int X1, int X2, int dx) {
    double x = 1, f1 = 1;
    for (X1; X1 <= X2; X1 += dx) {
        f1 += (2 * static_cast<_int64>(X1) - 1) * pow(x, (2 * X1 + 1)) /  ;  // тут вообще не знаю что делать
        x = f1;
        cout << X1 << "\t" << f1 << endl;
    }
    return f1;
}

#include <iostream>
#include <iomanip>
#include <cmath>
 
using namespace std;
 
void foo(double, double, double, int);
double bar(double, double);
 
int main()
{
    double dx, X1, X2;
    int prec, space;
    cout << "Enter values for X1, X2, dx and precision" << endl;
    cin >> X1 >> X2 >> dx >> prec;
    space = prec + 4;
    cout << left << fixed << setprecision(prec);
    cout << "Table of results" << endl;
    cout << setw(space) << "X1" << setw(space) << "Result" << setw(space) << "Real result" << endl;
    foo(X1, X2, dx, prec);
}
 
void foo(double X1, double X2, double dx, int prec)
{
    double eps = pow(0.1, prec); // можно повысить prec для повышения точности
    int space = prec + 4;
    for (X1; X1 <= X2; X1 += dx)
        cout << setw(space) << X1 << setw(space) << bar(X1, eps) << setw(space) << asin(X1) << endl;
}
 
double bar(double x, double eps)
{
    const size_t maxn = 1e7; // спасаемся от бесконечного цикла при слишком большой точности
    double f = x, add = x, padd = -1;
    x = x * x;
    for (size_t n = 2; abs(add - padd) >= eps && n < maxn; n += 2) // на самом деле тут хранится 2n, а не n
    {
        padd = add;
        add *= x * (n - 1) / n;
        f += add / (n + 1);
    }
    return f;
}