Вычислить функции с использованием их разложения в степенной ряд. Функции и Перегрузка функций

#include <cstdlib>
#include <iostream>
#include <math.h>
#include <iomanip>
using namespace std;
 
unsigned long long factorial(unsigned long long i) //вычисление факториала
{
    if (i == 0) 
        return 1;
    else 
        return i * factorial (i - 1);
}
 
int main()
{ 
setlocale(LC_ALL, "Russian"); // функция локализации вызывается только один раз
double  a, b, e;
int N;
    cout<<" Данные по варианту:"<< endl;
    cout<<" Разложение в ряд функций f(x): S=1+(x/1!)+(x^2/2!)+...+(xn/n!)"<< endl;
    cout<<" Диапазон изменения аргумента: a=1 b=2"<< endl;
    cout<<" N = 15: "<< endl;
    cout<<" Функций f(x) - сумма ряда (для проверки): S=e^x"<< endl;
    cout<<" eps близкое к нулю положительное число, например, 0.1, 0.01, 0.001, 0.0001"<< endl;
    cout<<" Введите Данные:"<< endl;
do
    {
        cout<<" Введитеи интервал а:"<< endl; cin>>a;// интервал
        cout<<" Введите интервал b:"<< endl; cin>>b;// интервал
    }while(a>=b);//проверка условий 
do
    {
        cout<<" Введите N:"<< endl; cin>>N;// диапазон изменения аргумента
 
    }while(N<1);//проверка условий 
do
    {
        cout<<"Введите eps: "; //точность 
        cin>>e;
    } while(e>0.1 || e<0.0001); //0.0001 < e < 0.1
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<"x"; // таблица
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<"Sn";
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<"Se";
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<"S"<<endl;
float x,h=(b-a)/10,sn=0,se=0,t,y; //вычисление шага
int i,j;
for(x=a;x<=b;x+=h)//вычисление x
{
    sn=x;se=x; //вычисление присваиваем значение х
    for(i=1;i<=N;i++)
{
    sn+=pow(x,i)/factorial(i);//правельное ли вычисление по формуле и факториала?
};
    j=1;
do 
{
    t=pow(x,j)/factorial(j);
    j++;se+=t;
} while(t>e);
    y=exp(x);
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<x; // заполнение таблицы
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<sn;
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<se;
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<y<<endl;
}
system("PAUSE"); 
return 0;
}

float x=a,h=(b-a)/10; //вычисление шага
for(int i = 1;i<=10;i++)//вычисление x
{
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<x; // заполнение таблицы
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<Myexp(x, N);
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<Myexp(x, e);
    cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<exp(x)<<endl;
}

h=(b-a)/10; //вычисление шага
for(x=a;x<=b;x+=h)//вычисление x
    {
        cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<x; // заполнение таблицы
        cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<My_exp(x, N);
        cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<My_exp(x, e);
        cout<<setw(15)<<setprecision(6)<<setiosflags(ios:: fixed)<<exp(x)<<endl;
    }