рекуррентная формула

//Подключение библиотеки стандартного ввода и вывода
#include<stdio.h>
//Подключение библиотеки математических функций
#include<math.h>
 
int main()
{
    //Объявление и инициализация
    //переменных программы
    double e, eps=0.0001, x, sum = 0, fact=1;
    int j, i, xk, xn, dx, count=0;
    bool flag = false;
 
    do
    {
        //Приглашение к вводу x нач
        printf("xn = ");
        //Ввод значения 
        scanf("%d",&xn);
 
        //Приглашение к вводу x кон
        printf("xk = ");
        //Ввод значения
        scanf("%d",&xk);
 
        if(xn > xk)
        {
            printf("xn doljno bit' < xk!\n");
        }else
        {
            flag = true;
        }
    }while (!flag);
 
    //Приглашение к вводу dx
    printf("dx = ");
    //Ввод значения
    scanf("%d",&dx);
 
    //Приглашение к вводу x
    printf("x = ");
    //Ввод значения
    scanf("%lf",&x);
 
    //Выводим шапку на экран
    printf("+-----------------------------------------------------------------------+\n");
    printf("| znachenie argumenta\t| znachenie funkcii\t| kol-vo prosumm el-tov\t|\n");
    printf("+-----------------------------------------------------------------------+\n");
 
    //Рассчитываем значения функции
    for(i = xn; i <= xk > eps; i+=dx)
    {
        fact = 1;
        //Цикл вычисления значения факториала числа
        for(j=2;j<=i;j++) 
            fact*=(double)j;
 
        e = (pow((-1), i) * pow(x, 2*i)) / fact;
        
        sum += e;
        count++;
 
        printf("|\t   %d \t\t|\t%lf \t|\t  %d  \t\t|\n", i, e, count);
 
    }
 
    printf("+-----------------------------------------------------------------------+\n\n");
    printf("e^-%.1lf^2 = %.3lf\n", x, sum);
 
    return 0;
}