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

#include <iostream>
#include <cmath>
#include <cstdlib>
#include <iomanip>
using namespace std;
 
int main(){
    double a, b, an, h, Sx, Yx;
    int n;
    a = 0.1;
    b = 1;
    h = (b - a) / 10;
    n = 80;
    cout << " x     Y(x)    S(x)\n";
    for (double x = a; x <= b; x += h){
        an = x;
        Sx = 0;
        for (int i = 1; i < n; i++){
            an *= pow(1, i) * (pow(x, 2 * i) / (i * 2));
            Sx += an;
        }
        Yx = cos(x);
        cout << setprecision(2) << fixed << x << setprecision(4) << "  " << Yx << "  " << Sx << "\n";
    }
 
}

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{    
    int n=80;    
    double a=0.1, b=1., h=(b-a)/10., an, Y, S;   
    b+=h/2.;
    
    cout << " x       Y(x)         S(x)\n";
    for (double x = a; x <= b; x+=h)
    {
        Y=cos(x);
        an=S=1.;        
        for (int k = 1; k < n; k++)
        {
           an*=-x*x/(2.*k*(2.*k-1.));
           S+=an;           
        } 
        
    cout << setprecision(2) << fixed << x << "  " << setprecision(9) << Y << "  " << S << "\n";
    }   
    
system("pause");
return 0;
}