Нахождение погрешности при вычислении ln

#include <iostream>
#include <cmath>
#include <iomanip>
#include <vector>
 
const double U = pow(2.0, -53);
 
 
 
double z = 1.0;
const double x = (z - 1.0) / (z + 1.0);
const int N = 100;
 
 
double GM(int j)
{
    return (U * j) / (1 - j * U);
}
 
 
 
std::vector<double> Array(std::vector<double>& ep)
{
    std::vector<double> Result;
 
    
 
    double x_ = x;
 
    double error_ = 0.0;
 
    for (int i = 1; i < N + 1; ++i)
    {
 
        for (int k = 1; k < 2 * i - 1; ++k)
        {
            x_ *= x;
        }
 
        Result.push_back(x_ / (2 * i - 1));
 
        error_ = (Result[i - 1]) * GM(2 * (2 * i - 1));
 
        ep.push_back(error_);
 
        x_ = x;
    }
 
    return Result;
}
double SumArray(std::vector<double> arr, std::vector<double>& ep, double& error2, std::vector<double> sum_error3, double& error1)
{
    if (arr.size() > 1)
    {
        if (arr.size() % 2 != 0)
        {
            arr.push_back(0.0);
 
            ep.push_back(0.0);
        }
 
        std::vector<double> res; 
 
        std::vector<double> eps; 
 
        std::vector<double> error_sum;
 
        for (int i = 1; i < arr.size(); i += 2)
        {
            res.push_back(arr[i] + arr[i - 1]);
 
            if (arr[i] != 0.0)
            {
                eps.push_back((ep[i] + ep[i - 1]) / (1 - 2 * U));
 
                error_sum.push_back(res.back() / (1 - 2 * U));
 
            }
            else {
                eps.push_back(ep[i - 1]);
 
                error_sum.push_back(res.back());
            }
 
        }
 
 
        return SumArray(res, eps, error2, error_sum, error1);
    }
    else {
 
        error1 = ep[0];
        error2 = sum_error3[0];
        return arr[0];
    }
}
 
int main()
{
    setlocale(LC_ALL, "rus");
    std::vector<double> ep;
 
    std::vector<double> sum_error3;
 
    std::cout.precision(20);
    std::vector<double> res = Array(ep);
 
 
    double error1 = 0.0;
 
 
 
 
 
    double error2 = 0.0; 
    double error3 = 0.0;
    double sum = 0.0;
    sum = SumArray(res, ep, error2, sum_error3, error1);
    std::cout << 10 * sum << "\n";
    error3 = (x / (1 - x * x));
 
 
 
 
    error3 = error3 / (2 * N + 1) * U;
    
    std::cout << error2 << " Ошибка суммирования\n";
    std::cout << error1 << " Погрешность\n";
   
    std::cout << error3 << " Ошибка отбрасывания остаточного члена ряда\n";
    std::cout << error1 + error2 + error3 << " Суммарная ошибка округления\n";
 
 
    return 0;
};

std::vector<double> getResultArray(double _x, double _n)   
{
    std::vector<double> result;
    result.reserve(_n);
 
    double mult_x = _x;
    int odd = 1;
    for(int i=0; i<_n; ++i) {
        result.push_back(mult_x / odd); //заполняем ряд: x ; (x^3)/3 ; (x^5)/5 ; (x^7)/7 ... (x^(2n-1)) / (2n-1)
        odd += 2;
        mult_x *= _x * _x;
    }
    return result;
}
 
 
std::vector<double> getEpArray(const std::vector<double> &_arr, double _u) 
{
    int n = _arr.size();
    std::vector<double> ep;
    ep.reserve(n);
 
    int odd_2 = 2;      //2, 6, 10, 14 ... 2(2n-1)
    for(double val : _arr) {
        double odd_2_u = _u * odd_2;
        ep.push_back(val * odd_2_u / (1 - odd_2_u));  
        odd_2 += 4;
    }
    return ep;
}