Функция вычисления суммы ряда ln(x)

#include <iostream>
#include <iomanip>
#include <stdio.h>
#include <conio.h>
#include <time.h>
#include <stdlib.h>
#include <locale.h>
#include <math.h>
using namespace std;
double sum (double x, double eps)
{
double s=0,sn=1.,z;
int n=1;
z=(x-1)/(x+1);
for (s=0., sn=1,n=1;fabs(sn)>eps;n++)
{
    sn=2/(2*n-1)*pow (z,(2*n-1));
    s+=sn;   
}
 
return s;
}
int main () 
{
    double x;
    for (x=0.1;x<=1.;x+=0.1)
        cout<<"x="<<x<<" sum="<<setprecision(4)<<sum(x,0.0001)<<" Ln(x)= "<<setprecision(4)<<log(x)<<endl;
    
    getch();
    return 0;
}

2.0/(2*n-1.0)

double my_ln (double x)
{
    const double eps=1e-6;
    const double lnln=M_LN2; //ln2
    int k=0;
    while(x>2.0)
    { 
        x/=2.0;
        k++;
    }
    x-=1.;
    double s=0;
    int n=1;
    double an=x;
    while (fabs(an)>eps)
    {
        s+=an;
        ++n;
        an*=-x*(n-1)/n;
    }
    s+=k*lnln;
    return s;
}

double my_ln (double x)
{
     if(x<1)
        return 0;
    const double eps=1e-6;
    double s=0;
    int n=0;
    double an=2.*x;
    double z =(x-1.)/(x+1.);
    while (fabs(an)>eps)
    {
        s+=an;
        ++n;
        an*=z*z*(2.*n-1.)/(2.*n+1);
    }
    return s;
}