Сумма ряда с заданной точностью - C++

#include <stdio.h>
#include <math.h>
 
double series(double x, double epsilon, int *n, double *last);
 
int main()
{
    double x;
    double epsilon;
    double last;
    int n;
 
    while (1)
    {
        printf("Enter x (-25 < x < 25): ");
        scanf("%lf", &x);
 
        if (x < -25.0 || x > 25.0)
            printf("x must be in range of -25 to 25");
        else
            break;
    }
 
    while (1)
    {
        printf("Enter epsilon (0 < epsilon <= 1): ");
        scanf("%lf", &epsilon);
 
        if (epsilon <= 0.0 || epsilon > 1.0)
            printf("epsilon must be in range of 0 to 1");
        else
            break;
    }
 
    printf("Sum of series: %lf", series(x, epsilon, &n, &last));
    printf("\nNumber of terms: %d", n);
    printf("\nLast term: %lf", last);
 
    getch();
    return 0;
}
 
double series(double x, double epsilon, int *n, double *last)
{
    double sum;
    double an;
    int den;
    int k;
    int sign;
 
    sum = 0.0;
    k = 4;
    den = 1;
    sign = -1;
 
    do
    {
        an = sign * x * (k + x);
        an /= 2 * den;
        k += 4;
        den *= k / 4;
 
        if ((k / 4) % 2 == 0)
            sign *= -1;
 
        sum += an;
    }
    while (fabs(an) >= epsilon);
 
    *n = k / 4;
    *last = an;
 
    return sum;
}