написать программу вычисления функции представленной бесконечным рядом чисел с точностью е=0.0001

#include "stdafx.h"
#include "iostream"
 
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{ 
    double x = 0, eps = 0, expf = 0, p = 0, xn = 0, xk = 0, dx = 0;
    int n = 0;
    cout << "xn, xk, dx, eps \n";
    cin >> xn >> xk >> dx >> eps;
    x = xn;
 
    cout << "+-----------+-----------+-------------+\n";
    cout << "|     x     |znac. f-cii|     n       |\n";
    cout << "+-----------+-----------+-------------+\n";
    while (x <= xk)
    {
    expf= 1;
    n= 1;
    p = 1;
    while (abs(p) >= eps)
    {
        p = p*x / n;
        expf = expf + p;
        n = n + 1;
    }
    cout << "  "<< x << "          " << expf << "          " << n << "\n";
     x = x + dx;
    }
    cin >> n;
 
 
    return 0;
}

double sum=1,an=-x/3.;
int n=2;
while(fabs(an)>1e-5)
{
   sum+=an;
   an*=-(n+3.)/(n*3.)*x;
   n++;
}
cout<<sum;