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

#include<iostream>
#include<cmath>
 
 
double factorial(double k)
{
 
    double j = 1;
    for (double i = 1; i <= k; i++)
    {
        j *= i;
    }
    return j;
}
 
using namespace std;
int main() {
    setlocale(LC_CTYPE, "rus");
 
        cout << "введите точность вычислений" << endl;
        double e;
        cin >> e;
 
    int N;
    cout << "N: ";
    cin >> N;
 
        int x;
        cout << "x: ";
        cin >> x;
    
 
 
    for (double n = 0; n <= N; n++) {
 
        float y0 = x*pow(4*x, 2*n) / factorial(2*n) + 1.2*sqrt(1 + 2*x);
        float k = k + pow(4*x, 2*n) / factorial(2*n);
        float y = k*x + 1.2*sqrt(1 + 2 * x);
 
        if (fabs(y - y0) < e) {
            cout << "y: " << y;
            break;
        }
    }
    return 0;
}

int fact(int x) {
    int result = 1;
    
    for(int i(1); i <= x; ++i)
        result *= i;
        
    return result;
}