Вычислить значение функции разложением в ряд Тейлора - C (СИ)

#include <stdio.h>
#include <math.h>
 
double factorial(double x)
{
    return tgamma(x + 1.0);
}
 
int main()
{
    int n, i;
    double x, res = 0.0, cur;
    
    scanf("%lf %d", &x, &n);
    
    for(i = 0; i <= n; i++)
    {
        cur = pow(x, 2 * i + 1) / factorial(2 * i + 1);
        if(i & 1) cur *= (-1);
        res += cur;
    }
    
    printf("%lf", res);
    
    return 0;
}