Приближенные вычисления и рекуррентные соотношения

#include<iostream>
#include<cmath>
using namespace std;
 
//Рекурсивный расчет ф-ии
float rec_func(float x, int n) {
    int i;
    if (n < 1) {
        i = 2 * n - 1;
        x *= (x * x);
        return x / i  + rec_func(x, n - 1);
    }
    else {
        return x;
    }
 
}
 
//Итеративынй расчет ф-ии
float iter_func(float x, int n) {
    float y = 0, a = x;
    int b = 1;
    for (int i = 1; i <= n; i++) {
        y += (a / b);
        b += 2;
        a *= x * x;
    }
    y = 2 * y;
    return y;
}
 
int main() {
    setlocale(LC_ALL, "rus");
    float x, y1, y2, y3, a1, b1, h;
    int n;
 
    //Ввод данных для ф-ии
        cout << "Введите x ( |x| < 0 ) : " << endl;
        cin >> x;
        cout << "Введите n: " << endl;
        cin >> n;
        cout << "Введите промижуток [a, b]" << endl;
        cin >> a1 >> b1;
        cout << "Введите размер шага h: " << endl;
        cin >> h;
 
    y1 = iter_func(x, n); // Итеративная ф-ия
 
    y2 = log((1 + x) / (1 - x)); // Библиотечная ф-ия
 
    y3 = 2 * rec_func(x, n); //Рекурсивная ф-ия
        
    cout << endl << "x = " << x << "    " << "y1 = " << y1 << "    " << "y2 = " << y2 << "    " << "y3 = " << y3 << endl;
    cout << endl << "Табуляция ф-ии на [a, b]: " << endl;
 
    while (a1 <= b1) {
        y1 = iter_func(a1, n);
        y2 = log((1 + a1) / (1 - a1));
        y3 = 2 * rec_func(a1, n);
        cout << "x = " << a1 << "    " << "y1 = " << y1 << "    " << "y2 = " << y2 << "    " << "y3 = " << y3 << endl;
        a1 += h;
    }
    system("pause");
    return 0;
}

//Рекурсивный расчет ф-ии
double rec_func(double x, int i,int n) {
// i - текущий индекс 
    static double s=x,a=x,b=1.0;
    if (i == 1) { // начинать надо с i=1
        s=x;a=x;b=1.0;
    }
    if(i==n)
        return s;
    s+=a/b;
    a *= (x * x);
    b+=2.0;
    return rec_func(x, i+1,n);
}