复件 (2) 跨孔同步法程序.txt

求二极法低阻球体上视电阻率值的程序

#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);
 }