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

#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
 
double f1(double x) 
{
    return atanh(x);
}
 
double f2(double x, double e) 
{
    double result = 0;
    double summa = x; // первое слагаемое равно x
    int N = 0;
 
    // вычисляем следующее слагаемое
    double cheslit = pow(x, 2 * N + 1); 
    double znamen = 2 * N + 1; 
    summa += cheslit / znamen;
    return summa;
}
 
 
double f3(double x, double e) {
    double sum = 0; // инициализируем сумму нулем
    double term = x; // первое слагаемое равно x
    int N = 1; // начинаем с первого слагаемого
    double x_squared = x * x; // вычисляем x^2
 
    while (abs(term) >= e) { // проверяем условие сходимости ряда
        sum += term; // добавляем текущее слагаемое к сумме
        term *= (x_squared) / (2 * N + 1); // вычисляем следующее слагаемое
        N++; // переходим к следующему слагаемому
        x_squared *= x * x; // обновляем значение x^2
    }
 
    return sum;
}
 
int main() 
{
    int N; // количество точек
    double A, B; // начало и конец отрезка
    int k; // точность
 
    cin >> N >> A >> B >> k;
    double step = (B - A) / (N - 1);
 
 
    cout << fixed << setprecision(k); // установка точности вывода
 
    cout << "x\tf1(x)\tf2(x)\tf3(x)" << endl; // заголовок таблицы
 
    for (int i = 0; i < N; i++)
    {
        double x = A + i*step;
        // проверка, принадлежит ли точка отрезку [A, B]
        if (x > A || x < B) 
        {
            cout << x << "\t" << f1(x) << "\t" << f2(x, pow(10, -k)) << "\t" << f3(x, pow(10, -k)) << endl; 
        }
        else {
            cout << x << "\t-\t-\t-" << endl;
        }
    }
 
    return 0;
}

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
double f1(double x)
{
    return log((1.+x)/(1.-x))*0.5;// atanh(x)
}
double f3(double x, double e) {
    double sum = 0; 
    double term = x; 
    int N = 0; 
    double x_squared = x * x;
    while (fabs(term) >= e) { 
        sum += term; 
        N++; 
        term *= x_squared / (2. * N + 1.); 
        x_squared *= x * x; 
    }
    return sum;
}
int main() 
{
    int N; // количество точек
    double A, B; // начало и конец отрезка
    double k=1e-6; // точность 1e-6
    cout<<"points:";
    cin >> N;
    cout<<"A(-1<A<1):";
    cin >> A;
    cout<<"B(A<B<1):";
    cin >> B;
    cin.get();
    double step = (B - A) / (N - 1);
    cout << fixed << setprecision(6); // установка точности вывода
    cout << "    x\t\ty=f3(x)\t\tatanh(x)" << endl; // заголовок таблицы
    for (int i = 0; i < N; i++)
    {
        double x = A + i*step;
        double y3=f3(x,k);
        cout << x << "\t" << y3<<'\t'<<f1(x) << endl; 
    }
    cin.get();
    return 0;
}