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

#include<stdio.h>
#include<iostream>
#include<cmath>
#include<math.h>
using namespace std;
int main() {
    setlocale(LC_ALL, "rus");
    long double f, Sx;
    int N,  k;
    float x, Yx, xn, xk,n1, h;
    Yx = 0;
    Sx = 0;
    f = 1;
    printf("Введите xn= "); scanf_s("%f", &xn);
    printf("Введите xk= "); scanf_s("%f", &xk);
    printf("Введите N="); scanf_s("%d", &N);
    h = (xk - xn) / 10;
    for (x = xn; x <= xk; x = x + h)
    {
        Yx = cos(x);
        for (n1 = 0; n1 <= N; n1++)
        {
            for (k = 1; k <= n1; k++)
            {
                f = f * k;
            }
            float w=-1;
            int q=pow(w,n1);
            Sx = Sx + ((q)*pow(x,(2*n1))) /(2*f); 
        }
        printf(" При x=%f, Yx=%f, Sx=%f, \n", x, Yx, Sx);
        
        Yx = 0;
        Sx = 0;
        f = 1;
    }
    system("pause");
}

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{    
    int N;    
    double xn, xk, h, an, S, eps=1e-9; 
    cout <<"xn="; cin >>xn;
    cout <<"xk="; cin >>xk;    
    cout <<"N="; cin >>N;
    h=(xk-xn)/10;    
    
    cout <<"   x    Y(x)    S(x)\n";
    for (double x=xn; x<=xk+eps; x+=h)
    {
        an=1.0;
        S=1.0;
        for (int k=1; k<N; k++)
        {
           an*=-x*x/(2*k*(2*k-1));
           S+=an;           
        }    
    cout <<"x="<<x<<"   "<<cos(x)<<"   "<<S<<"\n";
    }   
    
system("pause");
return 0;
}

S*= an;
S = S * an;

#include<stdio.h>
#include<iostream>
#include<cmath>
#include<math.h>
using namespace std;
int main()
{   setlocale(LC_ALL, "rus");
    double f,Sx,eps=1e-9;
    int N, k;
    float x, Yx, xn, xk,n1, h,g;
    Yx=0;
    Sx=0;
    f=1;
    printf("Введите xn= ");
    scanf_s("%f", &xn);
    printf("Введите xk= "); 
    scanf_s("%f", &xk);
    printf("Введите N="); 
    scanf_s("%d", &N);
    h=(xk-xn)/10;
    for (x=xn; x<=xk+eps; x=x+h)
    {  
        Yx=cos(x);
        for (n1=0; n1<=N; n1++)
        {
            for (k=1; k<=n1; k++)//проблема с этим факториалом
            {
                f=2*f*k;
                 
            }
            Sx=Sx+((pow((-1),n1)*pow(x,(2*n1)))/g);     
        }
        printf("При x=%f Yx=%f Sx=%f\n", x, Yx, Sx);
        Yx=0;
        Sx=0;
        f=1;
    }
    system("pause");
}

Sx=Sx+((pow((-1),n1)*pow(x,(2*n1)))/f);