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

#include <stdio.h>
#include <math.h>
#include <iostream>
#include <fstream>
#include <cmath>
#include <iomanip>
 
using namespace std;
int main() {
    setlocale(LC_ALL, "Russian");
    unsigned int N;
    double a, b;
    double znam = 1;
    double e = 2;
    cout << " Лабораторная работа №3. Вариант №9. Выполнена студентом 1 курса ДИПРБ Матвеевым Алексеем. \n";
    cout << "Введите а\n";
    cin >> a;
    while (a <= 0) {
        cout << " Введите а \n";
        cin >> a;
    }
    cout << " Введите b\n";
    cin >> b;
    while (b <= 0) {
        cout << " Введите b\n";
        cin >> b;
    }
    while (a>b) {
        cout << " Введите а <b \n";
        cin >> a;
        cout << " Введите  b\n";
        cin >> b;
    }
    cout << " Введите кол-во слагаемых \n";
    cin >> N;
    cout << "Введите точность \n";
    cin >> e;
    while (e>0.1 && e <= 0) {
        cout << "Введите точность (e>0.1 и e<=0)\n";
        cin >> e;
    }
    double h = (b - a) / 10;
    double x;
    cout << "x" << setw(10) << "S" << setw(10) << "Sn" << setw(10) << "Se" << endl;
    for (x = a; x <= b; x = x + h)
    {
        double x4;
        cout << x;
        x4 = x * x*x*x;
        double S = 0;
        S = (0.25*log((1 + x) / (1 - x)) + 0.5*atan(x));
        cout << setw(10) << S;
        double Sn = x, m = x;
        int i = 1;
        while (i <= N) {
            m = m * x4;
            Sn = Sn + (m) / (4 * i + 1);
            i++;
        }
        cout << setw(10) << Sn;
        double Se = x;
        i = 1;
        m = x;
        while (((m) / (4 * i + 1))>e) {
            m = m * x4;
            Se = Se + (m) / (4 * i + 1);
            i++;
        }
        cout << setw(10) << Se;
 
        cout << "\n";
    }
}

Se = Se + (m) / (4 * i + 1);