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

#include <iostream>
#include <cmath>
 
using namespace std;
 
int rec (double x, int i, int eps)
{
    double sum=0;
    while (fabs(sum)>eps)
    {
        sum += pow(x,4*i-3)*(4*i-x)/((4*i)*fact(4*i-1));
        i++;
    }
    return sum;
 
}
 
 int main()
 {
     double x, i(1), e, Rez1, Rez2, raz;
     cout << "Vvedite zadannuyu tochnost e: ";
     cin >> e;
     cout << "Vvedite peremennuyu x: ";
     cin >> x;
     Rez1 = rec(x, i,e);
     Rez2 = (2-sin(x)-cos(x)-(pow(exp,-x))/2*pow(x,2);
     cout << "V ryad: "<< Rez1;
     cout << "Formula: "<< Rez2;
     return 0;
 }

int main() {
    auto task = [](double x, double e)->double {
        const double x4 = x*x*x*x;
        auto abs = [](double x)->double {return x<0 ? -x : x;};
        auto go = [=](int j, double u, double d, auto& yl)->double {
            double d_ = d*j*(j-1)*(j-2)*(j-3);
            double a  = u * (j-x) / d_;
            return abs(a)<e ? 0 : a + yl(j+4, -u*x4, d_, yl);};
        return go(4, x, 1, go);};
 
    cout<<task(5, 1e-6)<<endl;
    return 0;
}