Sumatorias

//|Calcular Sumatorias|
#include <stdio.h>
#include <math.h>
using namespace std;
double const PI=3.14159265;

long factorial(int n){
    long factorial=1;
    for(int i=1;i<=n;i++){
        factorial*=i;
    }
    return factorial;
}

double sumatoria1(int k){
    double suma1=0;
    for(int n=1; n<=k; n++){
        suma1+=sin((n*(n+1))/factorial(2*n));
    }
    return suma1;
}

double sumatoria2(int k){
    double suma2=0;
    for(int n=0; n<=k; n++){
        suma2+=(PI/2)*((pow(-1, n))/(factorial((pow(n,2))+1)));
    }
    return suma2;
}

int main(){
    int j,k;
    printf("Ingrese El Limite Para la Primera Sumatoria:\n");
    scanf("%d", &j);
    printf("Ingrese El Limite Para la Segunda Sumatoria:\n");
    scanf("%d", &k);
    double P=sumatoria1(j);
    printf("Sumatoria P= %f\n", P);
    double Q=sumatoria2(k);
    printf("Sumatora Q= %f\n", Q);
    double S=P+Q;
    printf("Sumatoria S= %f\n", S);
}