c語言----編程

#include "math.h"

#include "malloc.h"

#define G 6.67 /*萬有引力常數（G：m^3/(kg*s^2))*/

main ()

{

FILE *fp;

int i;

double ra,rc,fia,fic;

double *Vxz,*Vzz;

double p,b,h,q; //double變量用%lf輸入可以用%f或%lf輸出

printf("\n***********A program of Gravitytabular***********\n");

printf("Remain Dentisity is p\n");

printf("Breadth is b\n");

printf("Burial Depth is h\n");

printf("Dip is q\n");

printf("Please enter the parameters:\n");

scanf("%lf%lf%lf%lf",&p,&b,&h,&q); //變量的順序和上面提示的順序不吻合

Vxz=(double *)calloc(50,sizeof(double));

Vzz=(double *)calloc(50,sizeof(double));

fp=fopen("answer2.txt","w");

for(i=0;i<=49;i++)

{

ra=sqrt(pow((i+b),2)+pow(h,2));

rc=sqrt(pow((i-b),2)+pow(h,2));

fia=atan((b+i)/h);

fic=atan((i-b)/h);

Vxz[i]=2*G*p*sin(q)*(sin(q)*log(rc/ra)+cos(q)*(fia-fic));

Vzz[i]=2*G*p*sin(q)*(sin(q)*(fia-fic)-cos(q)*log(rc/ra));

Vxz[i]=10*Vxz[i];

Vzz[i]=10000*Vzz[i];

fprintf(fp,"Vxz=%10.4f Vzz=%10.4f\n",Vxz[i],Vzz[i]);

}

fclose(fp);

printf("\n end of transform.\n");

}

//如果還有疑問我們可以交流交流

//祝妳學習進步