International 75 타원체의 총 면적과 반경을 계산하는 C 언어
1. 타원체의 면적을 얻기 위해 코드를 실행합니다.
International 75 타원체의 총 면적 : 510066100693931.625000 평방 미터 (510661 백만 평방 킬로미터)
지구 반경 Re = 6371010.172219 미터
2. 소스 코드의 일부
#include<stdio.h>
#include<math.h>
#define guoji75_a 6378140.00
#define guoji75_b 6356755.2881575287
#define PI 3.1415926535897932
#define keshi_a 6378245.0
#define keshi_b 6356863.0187730473
double guoji75quanmianji()
{
double P,L1,L2,B1,B2,a,b,e1;
L1 = 0.0;
L2 = PI;
B1 = 0.0;
B2 = PI/2;
a = guoji75_a;
b = guoji75_b;
e1 = (sqrt((a * a) - (b * b))) / a;
P=b*b*(L2-L1)*(((sin(B2)/(1-(e1*e1*sin(B2)*sin(B2))))+((1/(2*e1))*log((1+e1*sin(B2))/(1-e1*sin(B2)))))-((sin(B1)/(1-(e1*e1*sin(B1)*sin(B1))))+((1/(2*e1))*log((1+e1*sin(B1))/(1-e1*sin(B1))))));
return P;
}
main()
{
double a, e1,e2,M,N,R,_a,_b,A,B,X,P,x1,x2;
int weidu[3],fangweijiao[3],L1[3],L2[3],B1[3],B2[3],weidu1[3],weidu2[3];
int i,fg;
printf("下面计算国际75椭球全面积和与地球椭球面积相等的地球圆球的半径Re\n");
P=guoji75quanmianji();
printf("国际75椭球全面积为:%f平方米(%f亿平方公里)\n",2*P,2*P/100000000000000);
printf("地球圆球半径Re=%f米\n",sqrt(P/2/PI));
scanf_s("%d",&fg);
}