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

#include <iostream>
#include <cmath>
using namespace std;
int main()
{
    setlocale(LC_ALL, "Russian");
    double dx, x, x2, f, E, esp; int n;
    cout << "Введите Xнач. и Xкон. \n"; cin >> x >> x2;
    cout << "Введите dX\n"; cin >> dx;
    cout << "Введите e\n"; cin >> esp;
    while (x < x2)
    {
        n = 1;
        f = pow(x, n) / n;
        while (((pow(x, n + 1) / (n + 1)) - ((pow(x, n + 2) / (n + 2)) > E)))
        {
            f = ((pow(x, n + 1)) / (n + 1) - ((pow(x, n + 2)) / (n + 2)));
            n = n + 1;
            
            cout << "   " << x << "   " << x2 << "   " << E << endl;
 
            //E = round(E * 100000) / 100000; // 0.01
        }
        cout << "f= " << f << " n= " << n << "      x= " << x << endl;
        x = round(x * 100) / 100; // 0.01
        x = x2 + dx;
    }
}

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;
}