#include <stdio.h>
#include <math.h>
double P(int N, double X)
{   
   if(N == 0)
       return 1;
   else
       if(N== 1)  
           return X;
       else
           return ((2*N-1)*X*P(N-1, X) - (N-1) * P(N-2, X))/N;                                    
}




double Sum(double a,double b,double c,double d)

{ int n;   double sum=0.00000000; double l,k,C1,D1;/* C1為r/r0，D1為d/r0 */

  for(n=0;n<3;n++)
	{C1=sqrt((0.5*a+c)*(0.5*a+c)+d*d);D1=sqrt((0.5*a-c)*(0.5*a-c)+d*d);
       k=(D1*D1+C1*C1-a*a)/(2*D1*C1);
       l=pow(D1,-(n+1))*pow(C1,-(n+1));
       sum=sum+n*(b-1)/(n+(n+1)*b)*l*P(n, k);
	}

  printf("%f\n",sum ); 
  return sum ; 
}

 void main()
{ double A=3.9 ;double B=0.0700000; double C=0.55;double D=1;double y;/* A=R/ro，B=P2/P1，C=x/r0，D=h0/r0；，其中，R為AM之間的距離，r0為球體半徑，h0為球體埋深，P1，P2分別為圍巖與球體的電阻率值。x為oo’的距離，o為求球心在地面上的投影點，o'為AM的中點*/
 
y= 1+2*A*Sum(A, B, C, D);
 printf("%f\n", y);
 }