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

#include<iostream>
#include<cmath>
 
long double factorial(int x) {
    int res = 1;
    for (int i = 1; i <= x; ++i) {
        res *= i;
    }
    return res;
}
 
long double okruglenie(long double x, long double e) {
    int count = 0;
    long double r = e;
 
    for (int i = 1; r < 1; ++i) {
        r *= 10;
        count = i;
    }
 
    return (static_cast<long long>(x * pow(10, count))) / pow(10, count);
}
 
long double func(long double x, long double e) {
 
    long double res = 1.0;
 
    for (double i = 1.0; fabs(pow(-1.0, i) * pow(2.0, (2.0 * i - 1.0)) * pow(x, 2.0 * i) / factorial(2.0 * i)) > e ; ++i) {
 
        res += pow(-1, i) * pow(2.0, (2 * i - 1)) * pow(x, 2 * i) / factorial(2 * i);
    }
 
    return res;
}
 
long double func1(long double x, long double e) {
 
    long double res = 1.0;
    long double tmp = 0.5;
 
    for (long double i = 1; abs(tmp) > e; ++i) {
 
        tmp *= (-1) * 4 * pow(x, 2) / ((2 * i + 1) * (2 * i + 2));
 
        res += tmp;
    }
 
    return res;
}
 
int main() {
    long double a, b, n, e;
    std::cin >> a >> b >> n >> e;
    for (long double i = a; i < b; i = i + (abs(a - b) / n)) {
        double x = cos(i);
 
        std::cout << okruglenie(pow(x, 2), e) << '\t';
        std::cout << okruglenie(func(i, e),e) << '\t';
    /*  std::cout << okruglenie(pow(x, 2), e) << '\t';*/
        std::cout << okruglenie(func1(i, e), e) << '\t';
        if (okruglenie(func(i, e), e) == okruglenie(pow(x, 2), e)) {
            std::cout << "YES" << std::endl;
        }
        else {
            std::cout << "NO" << std::endl;
        }
    }
}

double my_func (double x, double eps = 1e-6)
{
  double sum = 0, memb = 1, tmp = 0.5;
  unsigned long long i = 0;
  while (fabs(memb) >= eps)
  {
    sum += memb;
    i += 2;
    tmp *= -4 * x * x / i / (i - 1);
    memb = tmp;
  }
 
  return sum;
}