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

#include "stdafx.h"
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <iomanip>
 
 
int main()
{   setlocale( LC_ALL,"Russian" );
    using namespace std;
    const float eps = 0.001f;      
    float x, chisl = 1, znam = 1, step=1, y=1, y0, y1;    
    int n=0;
    cout << "Введите x:  " << endl; 
    cin>> x;
    do   
    {
    chisl*=(x-1);
    znam*=(x+1);
    y0=y;
    y +=((chisl/znam)/n+1);       
    n++;  
    }   
    while (abs(y - y0) <= eps); 
    y1 = log(x);
    cout<< "Результат итеративного калькулятора: " << setw(9) << setprecision(5) << y <<endl;
    cout<< "Результат встроеного калькулятора: " << setw(9) << setprecision(5) << y1 <<endl;
    cout<< "Число иттераций: " << setw(9) << setprecision(5) << n <<endl;
    system ("PAUSE");
}

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
 
int main()
{   
    setlocale( LC_ALL,"Russian" );
    const double eps = 1e-6;      
    cout << "Введите x:  " << endl; 
    double x;
    cin>> x;
    double y=0;
    double an=(x-1.)/x;
    int n=0;
    do   
    {
        y+=an;
        an*=(x-1.)/x*(n+1.)/(n+2.);
        n++;  
    }   
    while (abs(an) > eps); 
    double y1 = log(x);
    cout<< "Результат итеративного калькулятора: " << setw(9) << setprecision(5) << y <<endl;
    cout<< "Результат встроеного калькулятора: " << setw(9) << setprecision(5) << y1 <<endl;
    cout<< "Число иттераций: " << setw(9) << setprecision(5) << n <<endl;
    system ("PAUSE");
    return 0;
}