Вычисление суммы ряда с использованием метода рекуррентных соотношений

#include "stdafx.h"
#include <math.h>
#include <iostream>
#include <conio.h>
#include <iomanip>
 
using namespace std;
 
void main()
{
    float S, Xn, Xk, Y, X, Dx, n, a, p;
    int N = 1;
    cout << "\n  N                   X                   S                   Y" << "\n";
    Xn = 0.1;
    Xk = 1;
    X = Xn;
    Dx = (abs(Xk - Xn)) / 9;
    while (N <= 10)
    {       
        a = 1;
        S = 0;
        n = 1;
        while (n <=10)
        {
            S = S + a;
            n++;
            p = (2 * n + 1) * pow(X, 2 * n) / n * (n - 1);
            a = a * p;
        }
        Y = (1 + (2 * (X * X))) * exp(X * X);
        cout << setw(10) << N;
        cout << setw(20) << setprecision(5) << X;
        cout << setw(20) << setprecision(5) << S;
        cout << setw(20) << setprecision(5) << Y;
        cout << "\n";
        X = X + Dx; 
        N = N + 1;
    }
    _getch();
}

#include "stdafx.h"
#include <math.h>
#include <iostream>
#include <conio.h>
#include <iomanip>
 
using namespace std;
 
void main()
{
    float S, Xn, Xk, Y, X, Dx, n, a, p;
    int N = 1;
    cout << "\n  N                   X                   S                   Y" << "\n";
    Xn = 0.1;
    Xk = 1;
    X = Xn;
    Dx = (abs(Xk - Xn)) / 9;
    while (N <= 10)
    {       
        a = 1;
        S = 0;
        n = 1;
        float f=1;
        while (n <=10)
        {     
            f*=n;
            p = (2 * n + 1) * pow(X, 2 * n) / f;
            n++;
            S = S + a;
            a =  p;
        }
        Y = (1 + (2 * (X * X))) * exp(X * X);
        cout << setw(10) << N;
        cout << setw(20) << setprecision(5) << X;
        cout << setw(20) << setprecision(5) << S;
        cout << setw(20) << setprecision(5) << Y;
        cout << "\n";
        X = X + Dx;
        N = N + 1;
    }
    _getch();
}