Программирование вложенных циклов

#include <iostream>
#include <cmath>
using namespace std;
int main() 
{
    double a, b, h, s = 0, y, x1, i, per, ans;
    cout << "a= ";
    cin >> a;
    cout << "b= ";
    cin >> b;
    cout << "h= ";
    cin >> h;
    for (double x1 = a; x1 < b; x1 += h)
    {
        y = exp(2*x1);
        i = 0;
        s = 0;
        ans = 0;
        {
            //int fact=1, n;
            //cin>>n;
           // for (n=1; n<=i; n++)
            //{
               // fact=fact*n;
           // }
           // cout << fact;
            while (fabs(s) > 0.0001 || i == 0) {
                s = (pow(2*x1, i)) / fact;
                ans += s;
                i++;
                
            }
            cout << "x= " << x1 << " y= " << y << " sum s= " << ans << endl;
 
            
        }
        
    }
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{    
    int n;    
    double a=0.1, b=1., h=0.05, e=0.0001, an, S, Y;
    b+=h/2.;
    
    cout <<" x       S(x)     Y(x)\n";
    for (double x = a; x <= b; x+=h)
    {        
        n=1;
        an=S=1.;        
        while(fabs(an)>=e)
        {           
           an*=2.*x/n;
           S+=an;
           n++;
        } 
        
        Y=exp(2.*x);
        
        cout << setprecision(2) << fixed << x << "  " << setprecision(6) << S << "  " << Y <<"\n";
    }   
    
system("pause");
return 0;
}