Вывести на экран таблицу значений функции Y(x) и ее разложения в ряд S(x) для x, изменяющегося от a до b с шагом

#include
#include
#include
int main()
{
double a,b,h,x,y,s,k;
int n,i;
cout << "Vvedite a,b,n\n";
cin >> a >> b >> n;
x=a;
h=(b-a)/10;
while (x<=b)
{
s=k=1;
y=(1-pow(x,2)/2)*cos(x)-x/2*sin(x);
for (i=1; i<=n; i++)
{
k/=2*i;
s += pow(-1,i)*(2*pow(i,2)+1)*pow(x,2*i)*k;
}
cout << setw(15) << x << setw(15) << y << setw(15) << s << endl;
x += h;
}
return 0;
}