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

#include <math.h>
#include <iostream>
#include <stdio.h>
 
using namespace std;
 
int main()
{
    setlocale(0, "");
    double xn, xk, dx, eps, x, z=0, mathF;
    int n = 0;
    cout << " Введите xn, xk, dx ( шаг ), eps ( точность ): ";
    cin >> xn >> xk >> dx >> eps;
    printf("-----------------------------------------------\n");
    printf("|     x    |      f     |    n   |   math F   |\n");
    printf("-----------------------------------------------\n");
    x = xn;
    double ch =  pow((x - 1), (2 * n + 1)) / ((2 * n + 1) * pow((x + 1), (2 * n + 1))), y = ch;
    
    while (fabs(ch) > eps && round(x * 100) / 100 <= xk) {
        
            ch = (pow((x - 1), (2 * (n + 1) + 1)) / ((2 * (n+1) + 1) * pow((x + 1), (2 * (n+1) + 1)))) / (pow((x - 1), (2 * n + 1)) / ((2 * n + 1) * pow((x + 1), (2 * n + 1))));
            mathF = log(fabs(x));
            z += mathF;
            y += ch;
            
        
        printf("| %3.2lf | %10.5lf   |  %3d   | %10.5lf |\n", x, 2*y, n, z);
        x += dx;
        n++;
        
    }
    printf("-----------------------------------------------\n");
}

#include <iostream>
#include <iomanip>
#include <cmath>
 
using namespace std;
 
int main()
{    
    int n;
    double xn, xk, dx, u, an, F, Fmath, eps;
    
    cout << "xn="; cin >> xn;
    cout << "xk="; cin >> xk;
    cout << "dx="; cin >> dx;
    cout << "eps="; cin >> eps; 
    xk+=dx/2.;
    
    cout << "   x       F        Fmath\n";
    for(double x=xn; x<=xk; x+=dx)
    {
        u=(x-1.)/(x+1.);
        an=u; 
        F=an;
        n=1;
        while (an >= eps) 
        {            
            an*=u*u/(2.*n+1.);
            F+=an;
            n++;
            if(n==501) break;
        }
        
        F*=2.;
        
        Fmath=log(x);
        
        cout << setprecision(2) << fixed << setw(5) << x << setprecision(4) << setw(10) << F << setw(10) << Fmath << "\n";
        n++;        
    }
    
system("pause");
return 0;
}