Оптимизация разложения функции

double Razl(double x, int k)
{
    double x1 = pow(-1, k + 1);
    double x2 = pow(x - 1, k) / k;
    return x1 * x2;
}

double Razl(double x, int k)
{
    double x2 = pow(x - 1, k) / k;
    return x2 * (1.0 - ((k & 0x01) * 2));
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{    
    int n;
    double x, an, S, eps;    
    
    do 
    { 
       cout <<"0<x<2  x="; cin >>x;
       if (x <= 0 || x >= 2) cout <<"0<x<2 !!! Repeat please...\n";
    } 
    while (x <= 0 || x >= 2);
    
    do
    {
       cout <<"0<eps<1  eps="; cin >>eps;
       if (eps<=0 || eps>=1) cout <<"0<eps<1 !!! Repeat please...\n";
    }
    while (eps<=0 || eps>=1);
    
        n=2;
        an=x-1;
        S=0.0;
        do 
        {
            S+=an;            
            an*=-(x-1)*(n-1)/n;            
            n++;               
        }
        while(fabs(an)>=eps);
    cout << " x      S(x)     ln(x)  |S(x)-ln(x)|  N \n";     
    cout << x << "  " << S << "  " << log(x) << "  " << fabs(S-log(x)) << "   " << n-2 <<endl;      
    
system("pause");
return 0;
}