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

double a, b, h, x, y, s, p;
    int n, i;
    cout << "Vvedite a,b,n" << endl;
    cin >> a >> b >> n;
    x = a;
    h = (b - a)/10;
    do
    {
        p = s = 1;
        for (i = 1; i <= n; i++)
        {
            p*=2*(x/2);
            s += p;
        }
        y = (((x*2)/4)+(x/2)+1)*exp(x/2);
        cout << setw(15) << x << setw(15) << y << setw(15) << s << endl;
        x += h;
    } while (x <= b + h / 2);
    cout<<"END"<< endl;

#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
 
int main()
{
    double a, b, h, x, y, s, p;
    int n, i;
    cout << "Vvedite a,b,n" << endl;
    cin >> a >> b >> n;
    x = a;
    h = (b - a)/10.0;
    do
    {
        p = 1.0;
        double fact = 1.0;
        for (i = 1; i <= n; i++)
        {
            fact *= i;
            p += (i*i + 1) / fact * pow(x/2.0, i);            
        }
        y = (((x*x)/4.0)+(x/2.0)+1.0)*exp(x/2.0);
        cout << setw(15) << x << setw(15) << y << setw(15) << p << endl;
        x += h;
    } while (x <= b + h / 2.0);
    cout<<"END"<< endl;
}