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

[Андрей979]

Доброго времени суток!
Помогите вычислить сумму ряда на Си. Условие такое:
f(x)=∑ от i=0 до бесконечности ((-1)^i(4i+1))/(4i+1)(2i)!
значок "^" - степень

[<- DelpHIFanat ->]

Ну, насчёт, рассчитать до бесконечности, не знаю как, а, вот, до определённого числа n - знаю (если n >= 5 результат получается один и тот же). Если компилятор не тормозит программу после выполнения, нужно раскомментировать закомментированные строки.

#include <stdio.h>
//#include <conio.h>

unsigned long fakt(unsigned short);
long IntPower(int, int);

int main()
{
    unsigned short n;
    printf("n = ");
    scanf("%u", &n);

    unsigned short i;
    double Sum = 0;
    for (i = 0; i <= n; i++)
        Sum += (double)(IntPower(-1, i) * (4 * i + 1)) / ((4 * i + 1) * fakt(2 * i));

    printf("Sum = %f", Sum);

	/*puts("\n======================\n\nPress ANY KEY!!!");
	getch();*/
	return 0;
};

unsigned long fakt(unsigned short n)
{
    if (n < 2)
        return 1;
    else
    {
        unsigned long f = 1;
        unsigned short i;
        for (i = 1; i <= n; i++)
            f *= i;
        return f;
    }
}

long IntPower(int x, int n)
{
    if (n == 0)
        return 1;
    unsigned int i = 1;
    long Pow = 1;
    for (i = 1; i <= n; i++)
        Pow *= x;
    return Pow;
}

Для компилятора Microsoft Visual Studio (студия любит тупить) код должен выглядеть так:

#include <stdio.h>
#include <conio.h>

unsigned long fakt(unsigned short);
long IntPower(int, int);

int main()
{
    unsigned short n;
	unsigned short i;
	double Sum = 0;
    printf("n = ");
    scanf("%u", &n);

    for (i = 0; i <= n; i++)
        Sum += (double)(IntPower(-1, i) * (4 * i + 1)) / ((4 * i + 1) * fakt(2 * i));

    printf("Sum = %f", Sum);

	puts("\n======================\n\nPress ANY KEY!!!");
	_getch();
	return 0;
};

unsigned long fakt(unsigned short n)
{
    if (n < 2)
        return 1;
    else
    {
        unsigned long f = 1;
        unsigned short i;
        for (i = 1; i <= n; i++)
            f *= i;
        return f;
    }
}

long IntPower(int x, int n)
{
    unsigned int i = 1;
    long Pow = 1;
	if (n == 0)
        return 1;
    for (i = 1; i <= n; i++)
        Pow *= x;
    return Pow;
}