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

#include <iostream>
#include <stdlib.h>
#include <clocale>
#include <math.h>
//Наша рекурсивная функция
int fact(int n)
{
 
  if (n<2) return 1;
else return n*fact(n-1);
} //тут функция вызывает саму себя
using namespace std;
int main() {
    setlocale( LC_ALL, "Russian" );
    float x, Xn, Xk, dX, eps, y, fd, ft;  //Xn,Xk,dX-нач.конеч.знач.и шаг, eps-точность, у-знач.элемента ряда, fd-знач.элем.ряда,ft-знач.функции с помощью ряда.
    cout<< "\nВведите интервал функции, шаг функции и точность:\n";
    cout<<"Xn=";                  //Начальное значение интервала
    cin>>Xn;
    cout<<"Xk=";                 //Конечное значение интервала
    cin>>Xk;
    cout<<"dX=";                 //шаг аргумента
    cin>>dX;
    cout<<"eps=";               //точность
    cin>>eps;
    x=Xn;
    float e=exp(pow(-x,2));    //Вычисляем значение функции
 
    int n=0;
    ft=0;
    while(fd>eps) {                     //пока разность больше точности вып. цикл
    y=(pow(-1,n)*pow(x,2*n))/fact(n);  //Нный элемент ряда
    ft+=y;                             //Значение функции с помощью ряда
    fd=abs(e-ft);                      //Вычисляем модуль разности между значением функции и её значением
    n=n+1;                             // найденным с помощью ряда Тейлора
    };
    cout<<"x="<<x<<endl;
    cout<<"e="<<e<<endl;
    cout<<"fd="<<fd<<endl;
    cout<<"ft="<<y<<endl;
    cout<<"n="<<n<<endl;
 
 
    return 0;
 
};