Вычислить функцию через ряд Тейлора

#include <iostream>
#include <cstdio>
#include <math.h>
#include <clocale>
#include <float.h>
using namespace std;
int main()
{
    setlocale(LC_ALL, "Russian");
    double xn, xk, dx, x;
    double eps = DBL_EPSILON;
    int n=0;
    printf  ( "введите х начальное и конечное,|x|<1\n"  )  ;
    scanf  ( "%lf%lf", &xn  ,  &xk );
    printf  ( " введите шаг dx, точность e \n" ) ;
    scanf ( "%lf%lf", &dx  , &eps ) ;
    //-------------------------------------------------------------------
    printf (  "Вычисление функции через ряд Тейлора\n" ) ;
    printf (  "-------------------------------\n")  ;
    printf (  "|    x    |     F    |    n   |\n" ) ;
    printf (  "-------------------------------\n" ) ;
    //  Некоторый множитель T равен Cn+1 / Cn или (x*x)/3
    double arth = 0, y = arth;  //  первый член ряда и нач. значение суммы
    for ( x = xn; x <= xk; x += dx)
    {
        do {
            arth *= ( pow( x, (2*n+1) ) / ( 2*n + 1) ); // очередной член ряда
            y += arth;
            n++;
            printf ( "|%9.21f   |%9.21f    |%9.21f  |\n",  x,  y,  n );
            }
    while (fabs( arth ) > eps);
    }
 
 
 
return 0;
}

 arth *= ( pow( x, (2*n+1) ) / ( 2*n + 1) ); // очередной член ряда