Dg.rar_DG_台阶重力异常_重力异常

#include<stdio.h> #include<math.h> #define h 100 #define H 200 #define e 4 #define pi 3.1415926 #define G 0.0067 float fun(float x,float a) { float r1,r2,f1,f2,w,dg,m,n,k; r1=sqrt(h*h+x*x); {if (a==pi/2) r2=sqrt(H*H+x*x); else r2=sqrt(H*H+pow(x-(H-h)/tan(a),2)); } {if (x<0) f1=asin(h/r1); else if (x==0) f1=pi/2; else f1=pi-asin(h/r1); } {if (a==pi/2) w=x; else w=x-(H-h)/tan(a); } {if (w<0) f2=asin(H/r2); else if (w==0) f2=pi/2; else f2=pi-asin(H/r2); } m=(x+h*tanh(a))*sin(a); n=sin(a)*log(r2/r1)-(f2-f1)*cos(a); k=H*f2-h*f1; dg=2*G*e*(m*n+k); return (dg); //printf("%f\n",dg); } void main() {int j; float a,c,x[21]; FILE *fp; fp=fopen("D:\\数据.txt","w+"); scanf("%f",&c); a=c*pi/180; printf("当θ为%f时的计算结果\n",a); for(j=0;j<21;j++) {x[j]=100*j-1000; printf("%f\n",fun(x[j],a)); fprintf(fp,"%f\n",fun(x[j],a)); } }