Вычислить сумму ряда

сase 3: {
            int n, k, j, f;
            double x, y, p, s;
            cout << "Enter x: "; cin >> x;
            cout << "Enter y: "; cin >> y;
            cout << "Enter n: "; cin >> n;
            cout << endl;
            s = 0;
            for (k = 1; k <= n; k++) {
                p = 1;
                for (j = 1; j <= k; j++) {
                    double el = sin(y - j);
                    p = p * el;
            }
                f = 1;
                for (j = 1; j <= 2*k+1; j++) 
                     f = j * f;
                s = s + (pow(cos(x * k),2) / (f)) + p;
                cout << "p= " << p << endl;
                cout << "f= " << f << endl;
                cout << "s= " << s << endl;
                cout << endl;

сase 3: {
int n, k, j, f;
double x, y, p, s;
cout << "Enter x: "; cin >> x;
cout << "Enter y: "; cin >> y;
cout << "Enter n: "; cin >> n;
cout << endl;
s = 0;
for (k = 1; k <= n; k++) {
p = 1;//Произведение элементов ряда
for (j = 1; j <= k; j++) {
double el = sin(y - j);
p = p * el;/
}
f = 1;
for (j = 1; j <= 2*k+1; j++)
f = j * f;//Факториал
s = s + (pow(cos(x * k),2) / (f)) + p;
cout << "p= " << p << endl;
cout << "f= " << f << endl;
cout << "s= " << s << endl;
cout << endl;

#include <iostream>
#include <cmath> 
int main()
{
    int n;
    std::cout << "Enter n: "; std::cin >> n;
    double x, y;
    std::cout << "Enter x: "; std::cin >> x;
    std::cout << "Enter y: "; std::cin >> y;
    double f = 1.0, p = 1.0, s = 0.0;
    for (int k = 1; k <= n; k++)
    {
        p *= sin(y - k);
        f *= k;
        s += pow(cos(x * k), 2) / (2 * f + 1) + p;
        std::cout << "p = " << p << "; f = " << f << "; s = " << s << std::endl;
    }
    return 0;
}