Вычисление функции c помощью разложения в ряд

#include <iostream.h>
#include <iomanip>
#include <math.h>
 
 
int _tmain()
{
    return 0;
}
 
 
void main()
{
    double xn, xk, dx, eps, s, sum, x, matf, s1,f;
    int n=1;    
    cout<<"Vvedite xn xk dx eps\n";
    cin>>xn>>xk>>dx>>eps;
    cout<<"-----------------------------------------------\n";
    cout<<"|     x    |      f     |    n   |   math F   |\n";
    cout<<"-----------------------------------------------\n"<<endl;
    x=xn;
    n=0;
    s1=(-3,14/2);
    f=s1;
    matf= atan(x);
 printf("%10.2lf | %8.2lf   |  %3d   | %10.2lf |\n",x,f,n,matf);    
    while(fabs(s)>eps && n <= 500 && (x<xk))
    {
        
            s=(((-1)*(2*n+1))/((2*n+3)*x*x));
            f=f+s;
            matf = atan(x);
            x=x+dx;
           printf("%10.2lf | %8.2lf   |  %3d   | %10.2lf |\n",x,f,n,matf);
                    n++;
    }
    cout<<"-------------------------------------------------\n";
}

#include <stdio.h>
#include <math.h>
 
#define M_PI 3.14159265358979323846
 
float farctg(int n, float xn, float xk, float eps, float dx)
{
if (n == 500)
 return 0.0;
 
float coef = pow( (float) -1, (float) n) / ( 2*n + 1 );
float step = 2.0*n+1.0;
 
float result = coef * pow(xn, step);
 
//printf("%d: %f %f\n", n, coef, step, result);
 
return (result+farctg(n+1, xn, 0,0,0));
}
 
int main()
{
float i;
for(i=0.0; i<1.0;i+=0.01)
{
printf("%f: %f ", i, farctg(0, i,0,0,0));
printf("%f\n", atan(i));
}
 
return 0;
}

#include <iostream.h>
#include <iomanip>
#include <math.h>
 
 
int _tmain()
{
    return 0;
}
 
 
void main()
{
    double xn, xk, dx, eps, s, sum, x, matf, s1;
    int n=1;    
    cout<<"Vvedite xn xk dx eps\n";
    cin>>xn>>xk>>dx>>eps;
    cout<<"-----------------------------------------------\n";
    cout<<"|     x    |      f     |    n   |   math F   |\n";
    cout<<"-----------------------------------------------\n"<<endl;
    for(x=xn; x<=xk; x+=dx)
    {n=0;
    sum=-3.14/2;
    s=-1/x;
    sum+=s;
 
    matf= atan(x);  
 
    while(fabs(s)>eps && n <= 500)
    {
        
            s*=(((-1)*(2*n+1))/((2*n+3)*x*x));
            sum+=s;
            n++;
    }
printf("%10.2lf | %8.2lf   |  %3d   | %10.2lf |\n",x,sum,n,matf);
    }
 
    cout<<"-------------------------------------------------\n";
}