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

for (double x=xn; x<xk; x+=dx){
  double sum=M_PI/2.0;
  double power=-x; 
  int denominator=1;
  double step;
  do{
    step=power/denominator;
    power*=(-x*x);
    denominator+=2;
    step+=power/denominator;
    power*=(-x*x);
    denominator+=2;
    
    sum+=step;
  }while(fabs(step)>eps);
  std::cout<<"arcctg("<<x<<")="<<sum<<std::endl;
}

#include <conio.h>
#include <stdio.h>
#include <math.h>
 
using namespace std;
 
int main(int argc, char *argv[])
{
    double Xn, Xk, dx, e, x, q, s, y;
    
    while (true)
    {
        
        cout << "VVedite nachalnoe znachenie Xn" << endl;
        cin >> Xn;
        cout << "VVedite konechnoe znachenie Xk" << endl;
        cin >> Xk;
        cout << "VVedite shag dx" << endl;
        cin >> dx;
        cout << "VVedite pogreshnost' e" << endl;
        cin >> e;
        
        
        printf ("--------------------------------------------------------\n");
        printf ("       X     |    Levaya chast'   |   Pravaya chast'  | \n");
        printf ("--------------------------------------------------------\n");
        
        for (x = Xn; x <= Xk; x += dx)
        {
            s = 1;
            q = 1;
            
           y = exp(-x);
            
            for (int n = 0; fabs(q) > e; n++)
            {
                q *= (-1 * x) / (n + 1);
                s += q;
            }
            
            printf("|%9.1lf   |%9.5lf           |%9.5lf          |\n", x, y, s);
 
        } 
        
        printf ("-------------------------------------------------------\n");
    }
        
    getch();
    return 0;
}