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

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
int main() {
    setlocale(LC_ALL, "Russian");
    int N = 100;
    float x, eps;
    cout << "x="; cin >> x;
    cout << "eps="; cin >> eps;
    float a = x/2, sum = a; //начальные условия
    cout << left << setw(3) << "i\t" << setw(10) << "a\t"; //заголовок
    cout << setw(10) << "sum" << endl;
    cout << left << setw(3) << 1 << "\t" << setw(10) << a; //1 шаг
    cout << "\t" << setw(10) << sum << endl;
    for (int i = 2; abs(a) > eps && i <= N; i++) { //i-ый шаг
        sum += a;
        a *= x / (2*i-1) ;
        cout << left << setw(3) << i << "\t" << setw(10) << a;
        cout << "\t" << setw(10) << sum << endl;        
    }
    return 0;
    
}

#include <iostream>
#include <cassert>
#define ABS(x) ((x) < 0 ? -(x) : (x))
 
using namespace std;
 
long foo(double x, double precision, double &sum)
{
  assert(x > 0 && precision > 0); // по условию задачи
  double m = 0.5 * x;
  size_t n = 1;
  for (sum = 0; ABS(m) > precision; n++)
    sum += m, m *= x / ((2 * n + 1) * (2 * n + 2));
  return n;
}
 
int main()
{
  double x[]   = {2, 2, 0.5, 0.5};
  double eps[] = {1e-3, 1e-5, 1e-3, 1e-5};  
  
  for (auto idx = 0; idx != sizeof(x) / sizeof(*x); idx++)
  {
    double sum;
    long steps = foo(x[idx], eps[idx], sum);
    printf("# %4g | %4g | %.6f | %ld\n", 
           x[idx], eps[idx], sum, steps - 1); 
  }  
 
  return 0;
}