Floating point division by zero

[Nater]

Где slag*= при k=12 в знаменателе образуется ноль, хотя математически этого быть не должно. В чем прикол? не знаю как решить проблему.

//---------------------------------------------------------------------------

#pragma hdrstop

#include <tchar.h>
#include <conio.h>
#include <iostream.h>
#include <math.h>
//---------------------------------------------------------------------------

#pragma argsused
int u(double n)
{
	double c = 1;

	for (int i = 1; i <= n; i++)
	{
	c*=i*1.0;
	}
	return c;
}

int _tmain(int argc, _TCHAR* argv[])
{
	double x, E, S = 0, slag = 0;
	int k;

	cin >> x;
	cin >> E;
	k = 3;
	slag = x*x*x*x*x*x*7.0/24.0;
	do
	{
	S+=slag;
	k++;
	cout << k << "\n";
	cout << (k*1.0*u(2*k-2)*1.0*(u(2*k-4)*1.0-k*1.0+1.0)) << "\n";
	slag*=(-x*x*(((u(2*k-2)*1.0)-k*1.0)*(k*1.0-1)*u(2*k-4)*1.0))/(k*u(2*k-2)*(u(2*k-4)-k+1.0));
	//cout << slag << "\n";
	} while (fabs(slag)>E);

	cout << S << "\n";
	cout << log(1+x*x)-x*x*cos(x) << "\n";
	getch();
	return 0;
}
//---------------------------------------------------------------------------

[Базиля]

Сообщите значения остальных переменных в вашем выражении на шаге при k=12.
Также, хотелось бы увидеть формулу в более читаемом виде.