Почему не работает программа? Ряды Тейлора

#include <iostream>
#include <clocale>
#include <iomanip>
#include <cmath>
 
using namespace std;
double taylor(double x,double e,int & n)
{
double sum=0,a,b;
do
        {
            a=pow(x,2*n+1)/(2*n+1);
            b=pow(x,2*n+2)/((2*n+2));
            n++;
            sum+=a;
        }
    while (abs(b-a)>e);
return sum;
 
}
 
int main()
{
    setlocale(LC_ALL,"RUS");
    //Ввод и проверка данных
    double x1,x2,dx,e;
    cout<<"Введите начальное значение промежутка: ";
    cin>>x1;
    while (cin.fail())
    {
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(),'\n');
        cout << "Введено некорректное значение. Введите начальное значение промежутка: ";
        cin >> x1;
    }
    cout<<"Введите конечное значение промежутка: ";
    cin>>x2;
    while (cin.fail())
    {
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(),'\n');
        cout << "Введено некорректное значение. Введите конечное значение промежутка: ";
        cin >> x2;
    }
    cout<<"Введите шаг: ";
    cin>>dx;
    while (cin.fail())
    {
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(),'\n');
        cout << "Введено некорректное значение. Введите шаг: ";
        cin >> dx;
    }
    cout<<"Введите точность: ";
    cin>>e;
    while ((e>1)||(e<0)||cin.fail())
    {
        cin.clear();
        cin.ignore(numeric_limits<streamsize>::max(),'\n');
        cout << "Введено некорректное значение. Введите точность (положительное значение меньше 1): ";
        cin >> e;
    }
    cout<<"     Таблица значений функции     "<<endl;
    cout<<"X       Сумма ряда      Количество элементов"<<endl;
        //Проверка на существование промежутка
    if(x1<x2)
    {
        //цикл для расчета значений по заданному промежутку
        double i,znach;
        int n;
        
        for(i=x1;i<=x2;i+=dx)
        {
            
                n=0;
                znach=taylor(i,e,n);
                cout<<fixed<<setprecision(2)<<i<<setw(19) << setprecision(10)<<znach<<setw(15)<<n<<endl;
 
            
        }
    }
    else 
        cout<<"Введенный интервал не существует"<<endl;
    system("Pause");
    return 0;
}

for(i = x1;  ; i += dx)
{
            
       n=0;
       znach=taylor(i,e,n);
       cout<<fixed<<setprecision(2)<<i<<setw(19) << setprecision(10)<<znach<<setw(15)<<n<<endl;
       if (i > x2) break;
}