suanfa.txt

用复合辛普森公式计算卫星轨道长度的程序.

#include <stdio.h>
#include <math.h>
#define PI 3.14159265358979323846 
void main()
{double i,m,n,s;
	double A=(439.0+2384+6371*2)/2,C;
double a,b;
C=6371.0+2384-A;
a=A*A;b=a-C*C;
m=0;n=0;
for(i=2*PI/1000;i<2*PI;i+=2*PI/1000)
m=m+sqrt(a*sin(i)*sin(i)+b*cos(i)*cos(i));
m=m*2;
for(i=2*PI/2000;i<2*PI;i+=2*PI/1000)
n=n+sqrt(a*sin(i)*sin(i)+b*cos(i)*cos(i));
n=n*4;
s=2*PI/6000*(sqrt(b)+sqrt(b)+m+n);
printf("%lf\n",s);
}结果：
C=48723.610382km