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

#include <stdio.h>
#include <math.h>
using namespace std;
int main()
{
    setlocale(LC_ALL,".1251");
    double a,b,h,y,eps;
    cout<<"Введите начало отрезка a и конец отрезка b "<<endl<<"a=";
    cin>>a;
    cout<<"b=";
    cin>>b;
    if ((b<a)||(a<-1) || (a>1) || (b<-1) || (b>1));
    {
        cout<<"Неправильно задан промежуток"<<endl;
        return 0;
    }
    cout<<"Введите шаг h=";
    cin>>h;
    if (h>abs(a-b))
       {
           cout<<"Неправильно задан шаг"<<endl;  
           return 0;
       }
    cout<<"Введите точность вычислений k=";
    int k;
    cin>>k;
    eps=pow(10.0,-(k+2));
    cout.setf(ios::fixed | ios::left);
    cout.precision(k);
    cout.width(k+5);
    cout<<"x";
    cout.width(17);
    cout<<"Кол-во слагаемых";
    cout.width(k+5);
    cout<<"f(x)";
    cout.width(k+10);
    cout<<"Стандартная функция"<<endl;
    int n=1;
    for (double x=a ; x<b ; x+=h )
    { 
        double d=x;
        int slog =1;
    while (abs(d)>=eps)
    {
        n=n+2;
        d=(-1.0*d*x*x*(n-2))/n;
        slog++;
    }
        y=atan(x);
        cout.width(k+5);
        cout<<x;
        cout.width(17);
        cout<<slog;
        cout.width(k+5);
        cout.width(k+10);
        cout<<y<<endl;
        n=1;
    }
    cout<<endl;
    system("pause");
    cin.get();
    return 0;
}