Помогите решить задачу...

[molodzo]

Задача:

k изменяется от 0 до бесконечности ((4k+1)! / [(2k)!]^2) * x^2k = =1/2*[(1-4x)^(-3/2) + (1+4x)^(-3/2)]

X принадлежит интервалу [-1/4 ; 1/4)

Найти мин. k, которое будет удовлетворять условию:
abs(правая часть - левая часть) < 10^(-15)


Вот код(тока у меня не работает):

#include <iostream.h>
#include <conio.h>
#include <math.h>

long double k,x2,x,p,p2,n,znam,znam1,znam2,p3,my_res,p4,p5,p6,p7,p8,sum1,sum2,ch,sravn;
long double factorial(long double);

main()
{
k=-1;
x2=pow(1/10,15);
do
{
      my_res=1;
      x=-0.25;
      k=k+1;
      p=4*k+1;
      ch=factorial(p);
      p2=2*k;
      znam=factorial(p2);
      znam1=znam*znam;

         do
         {
            x=x+x2;
            if (2*k>=0) p3=pow(x,2*k);
            else p3=(1/(pow(x,abs(2*k))));
            sum1=ch/znam1*p3;

            p4=1-4*x;
            p5=(1/(pow(p4,3/2)));
            p6=1+4*x;
            p7=(1/(pow(p6,3/2)));
            p8=p5+p7;
            sum2=0.5*p8;
            sravn=pow(1/10,15);
         }
         while ((x<=0.25) && (abs(sum1-sum2)>=sravn));

         if(abs(sum1-sum2)<sravn)
            {
               cout<<"Sum1="<<sum1;
               cout<<"Sum2="<<sum2;
               cout<<"K="<<k;
               getch();
               break;
            }
         else continue;
      }
while(abs(sum1-sum2)>=sravn);
getch();
}

long double factorial(long double n)
{
   if (n==1) {exit(1);}
   else
   {
      my_res=my_res*n;
      factorial(n-1);
   }
}

Помогите, пожалуйста, ее решить!

[Svetlag]

Мне не очень понятна задача.
Но факториал я бы считала так

long double factorial(const long double n)
{
 	if(n<0) return 0.0;
	if(n==1 || n==0) return 1.0;
	else return n*factorial(n-1);
}

Вместо

abs(sum1-sum2)

поставьте

fabs(sum1-sum2)

И еще мне кажется вы задаете слишком маленький шаг по x.