Разложение arcctg(x) в ряд Тейлора (переделать программу)

#include<iostream>
#include<math.h>
#include<iomanip>
#define _USE_MATH_DEFINES
using namespace std;
void main(){
float xn,xk,dx,eps,x,s,f;
int i,n;
cout<<"\n xn= ";cin>>xn;
cout<<"\n xk= ";cin>>xk;
cout<<"\n dx= ";cin>>dx;
cout<<"\n eps= ";cin>>eps;
cout<<"\n     x      s(x)          arctgx        n\n";
n=0;
for(x=xn;x<=xk;x+=dx)
{
s=x;
for(f=x,i=3; fabs(f)>eps ; i+=2)
{
f=f*(-x*x/(i-1)/i);
s+=f;n=i;
}
cout<<"\n"<<setprecision(4)<<setw(10)<<x
<<setprecision(4)<<setw(10)<<s
<<setprecision(4)<<setw(10)<<arctg(x)
<<setw(10)<<n;
}
cout<<"\n";
system("PAUSE")
}

double a, b, eps;
cin >> a >> b >> eps;
 for (int i = 0; i < 10; ++i){
        double x = a + i*(b - a)/9;
        
        double dx = x;
        double y = dx;
        int n = 1;
        while(fabs(dx) > eps){
            ++n;
            dx *= -x*x*(2*n - 3)/(2*n - 1.);
            y += dx;
        }
        
        cout << setw(7) << fixed << setprecision(2) << x;
        cout << setw(15) << fixed << setprecision(5) << y;
        cout << setw(15) << fixed << setprecision(5) << atan(x) << endl;