Найти ошибку (Вычисление функции с помощью разложения в ряд)

#include <iostream>
#include <conio.h>
#include <clocale>
#include <cmath>
#include <iomanip>
 
 
using namespace std;
 
 
inline double my_abs(double val)
{
        double ret = val;
        if(val < 0)
                ret *= -1;
        return ret;
}
 
inline double Taylor(double x, double e, long &n)
{
        double sum = 0;
        double an = 1;
        n = 1;
        do
        {
                sum += an;
                n   += 1;
                an *= (-1)*(x*x/(2.0*n*2*(n - 1)));
        }
        while(e <= my_abs(an));
        return sum;
        
}
 
 
int main()
{
    setlocale(LC_ALL, "Russian");
 
    long n = 0;
    double xn, xk, dx, e;
 
    cout << "Программа для вычисления значения функции\ny = cos(x) с помощью ряда Тейлора\n";
 
    cout<<" xn = ";
    cin>>xn;
    cout<<"Конечное значение xk = ";
    cin>>xk;
    cout<<"Шаг dx = ";
    cin>>dx;
    cout<<"Погрешность e = ";
    cin>>e;
 
    cout << "-----------------------------------------------------------\n";
    cout << "|   x  | f(x) | f(x) в ряде Тейлора  |   n   | точность e | \n";
    cout << "-----------------------------------------------------------\n";
 
    while(xn <= xk)
    {
        
        cout<< "|" << setw(10)<< setprecision(8) << xn << "  |    "
            << setw(10) << setprecision(8) << cos(xn) <<"| "
            << setw(10) << setprecision(8) << Taylor(xn, e, n) <<"  | "
            << setw(10) << n << "|" 
            << setw(10) << e << "|" << endl;
        xn = xn + dx;
    }
    cout << "--------------------------------------------------------------------------\n";
 
    _getch();
    return 0;
}

inline double Taylor(double x, double e, long &n)
{
        double sum = 0;
        double an = 1;
        n = 1;
        do
        {
                sum += an;
                n   += 1;
                an *= (-1)*(x*x/(2.0*n*2*(n - 1)));
        }
        while(e <= my_abs(an));
        return sum;
        
}

#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
 
double Radian(double x);
double Period(double x);
double Taylor(double x, double e, long &n, double &an);
 
 
int main()
{
    long n;
    double x, y, xn, xk, dx, e, err;
 
    system("chcp 1251");//это нестандартная функция потому заменяйте на свою ф-цию руссификации
    cout<<"Программа для вычисления значения функции"<<endl
        <<"y = cos(x) с помощью ряда Тейлора"<<endl;
    cout<<"\tВвод значений углов в градусах"<<endl;
    cout<<"Начальное значение xn = ";cin>>xn;
    cout<<"Конечное значение  xk = ";cin>>xk;
    cout<<"Шаг аргумента      dx = ";cin>>dx;
    cout<<"Погрешность e = ";cin>>e;
    cout<<"Радианные меры углов : "<<endl;
    cout<<(xn = Radian(xn))<<endl;
    cout<<(xk = Radian(xk))<<endl;
    cout<<(dx = Radian(dx))<<endl;
    cout<<"Удаление целых периодов : "<<endl;
    cout<<(xn = Period(xn))<<endl;
    cout<<(xk = Period(xk))<<endl;
    cout<<(dx = Period(dx))<<endl;
 
    cout<<"|    x   |   N   |   Sum   |   cos   |   fabs(an)   |     e     |"<<endl;
    for(x = xn; x <= xk; x += dx)
    {
        y = Taylor(x, e, n, err);
        cout<<((x < 0) ? "|-" : "| ")
            <<setw(5)<<setprecision(3)<<fabs(x)<<" | "
            <<setw(5)<<n<<" | "
            <<((y < 0) ? "-" : " ")
            <<setw(6)<<setprecision(4)<<fabs(y)<<" | "
            <<((cos(x) < 0) ? "-" : " ")
            <<setw(6)<<setprecision(4)<<fabs(cos(x))<<" | "
            <<setw(8)<<setprecision(6)<<err<<" | "
            <<setw(8)<<setprecision(6)
            <<(cos(x) ? fabs((y - cos(x)) / cos(x)) : 0)<<" | "
            <<endl;
    }
    system("pause");//это нестандартная функция потому заменяйте на свою ф-цию стоп-точки
    return 0;
}
 
double Radian(double x)
{
    return x*asin(1.0) / 90;
}
 
double Period(double x)
{
    double period = 4*asin(1.0);
    while( x - period > 0)
        x -= period;
    return x;
}
 
double Taylor(double x, double e, long &n, double &err)
{
    double sum = 0;
    double an  = 1;
    n = 0;
    for(n = 0; e <= (err = fabs(an)); n++)
    {
        sum += an;
         an *= (-1)*x*x/((2*n + 1)*(2*n + 2));
    }
    return sum;
    
}