Вычисление ln(1 + x) разложением в ряд

#include <iostream>
#include <cmath> 
 
double ln(double x, int threshold){
        double result=0;
        for (int i=1;i<threshold;i++)
                result += (i%2==0?-1:1) * pow(x,i)/i;
        return result;
}
int main (){
 
    double x;
    int n;
 
    std::cout << "Input x: ";
    std::cin >> x;
 
    std::cout << "\nInput n: ";
    std::cin >> n;
 
    std::cout << "\n\nln = " << ln(x,n);
   
    return 0;
}

double eps=0.0000000001;
double ln(double a)
{
    double t = 0;
    double ans = 0, st = a, j = 1, zn = 1;
    double e = exp(1);
    while(a > sqrt(e))
    {
        a /= sqrt(e);
        t+=0.5;
    }
    a--;
    st = a;
    while(abs(st / j) > eps)
    {
        ans += zn * st / j;
        st *= a;
        zn *= -1;
        j++;
    }
    ans += t;
    return ans;
}

#include <iostream>
#include <cmath>
 
using namespace std;
double eps=0.0000000001;
double ln(double a)
{
    double t = 0;
    double ans = 0, st = a, j = 1, zn = 1;
    double e = exp(1);
    while(a > sqrt(e))
    {
        a /= sqrt(e);
        t+=0.5;
    }
    a--;
    st = a;
    while(abs(st / j) > eps)
    {
        ans += zn * st / j;
        st *= a;
        zn *= -1;
        j++;
    }
    ans += t;
    return ans;
}
int main()
{
    double y;
    cin>>y;
    cout<<ln(y);
}