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

#include "iostream"
#include "iomanip"
#include "cmath"
using namespace std;
void main()
{
    setlocale(0,"rus");
    int n=0;
    float xn,xk,dx,e,sum=0,t=1,PI=3.1415;
    bool prov;
 
    cout<<"Введите Xn";
    cin>>xn;
    cout<<"Введите Xk";
    cin>>xk;
    cout<<"Введите шаг dx";
    cin>>dx;
    cout<<"Погрешность e=";
    cin>>e;
 
    if (xn<-1 ||xn>1 ||xk<-1 ||xk>1 ||e==0)
    {
        cout<<"В связи с ОДЗ: |x|<=1. Введенные значения не могут быть использованы"<<endl;
    }
    else
    {
        cout<<"-----------------------------------------------------------"<<endl;
        cout<<"|Значение X|Значение суммы|Кол-во элементов|"<<endl;
        cout<<"-----------------------------------------------------------"<<endl;
    }
    while (xn<=xk)
    {
        sum=0,t=1;
        n=0;
        prov=false;
        do
        {
            t=(pow(-1.,n+1)*pow(xn,2*n+1))/(2*n+1);
            sum=sum+t;
            n++;
            if (n>=50000) prov=true;
        }
        while (e<=abs(t));
        cout<<"|"<<setw(15)<<xn<<"|"<<setw(15)<<PI/2+sum<<"|"<<setw(15);
        if(prov==false)
        {
            cout<<n;
        }
        else
        {
            cout<<"Ряд расходится";
        }
        ;cout<<"|"<<endl;
        xn=xn+dx;
        cout<<"------------------------------------------------------------------------"<<endl;
        system("pause");
    }
}

if (n>=50000) prov=true;