Вывести таблицу значений функции и ее разложения в ряд Тейлора на заданном отрезке с заданным шагом

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{    
    int n;    
    double a=0.0, b=2.0, h=0.1, an, sum, delta, e, eps=1e-9;        
 
    do
    {
       cout << "0<e<1 e="; cin >> e;
       if (e<=0 || e>=1) cout <<"0<e<1, Repeat please...\n";
    }
    while (e<=0 || e>=1);
    
    cout <<" x      sum      exp(x)    delta\n";
    for (double x=a; x<=b+eps; x+=h)
    {
        n=1;
        an=1.0;
        sum=0.0;
        while(fabs(an)>=e)
        {
           sum+=an;           
           an*=x/n;
           n++;
        } 
    delta=sum-exp(x);    
    cout <<setprecision(2)<<fixed<<x<<"  "<<setprecision(6)<<sum<<"  "<<exp(x)<<"  "<<delta<<endl;
    }   
    
system("pause");
return 0;
}