Borland C++, какая ошибка в задаче? (выбивает floating point divide by zero at 0x401197)

Задание
Для Х с данного диапазона (шаг выбрать самостоятельно) вычислить значение функции с помощью:
1)Стандартной функции –Y=f(X);
2)суммы ряда с точностью F (вибрати самостійно) .
Результат нужно вывести в виде таблицы:
іХY=f(X)S(X)fabs(s-y)

Код:

# include
# include

void main()
{
const F=0.1;
long int i,j,P2,P;
double z,x,P1,Y,S;
//scanf ("%g",&z);
z=0.15;
for (x=-1; x<=1; x+=z)
{
S=0;
P2=1;
i=0;
do
{
i++;
P2*=(-1);
P=1;
for (j=1; j<=i; j++)
P*=j;
P1=pow(2*x,2*i);
S+=P2*P1/P;
printf ("h=%g\n", P);
}
while (fabs(P2*P1/P)>=F);

Y=1+S/2;

printf ("Y=%g\n", Y);
}
getchar ();
getchar ();

}