C语言习题5.17--求地球表面任意两…

Description

假设地球是一个标准的球体，其半径是R，忽略地形对距离的影响。球面上的点使用经纬度表示，规定东经为正、西经为负、北纬为正、南纬为负。设A点的经度是α1、纬度是β1；B点的经度是α2、纬度是β2，则A和B之间的距离可按下式近似计算：同时约定：东经为正,西经为负；南纬为90°+地理纬度值、北纬为90°-地理纬度值D=Rarccos(C)其中：C=sin(β1)sin(β2)cos(α1-α2)+cos(β1)cos(β2)；特别注意：arccosC的单位是弧度.假设R=6371km，编写函数计算地球表面任意两点之间的距离。

Input

两个点的经纬度

Output

两点之间的距离

Sample Input

121.457413 37.476
116.30729 39.98772

Sample Output

526.709160

#include <stdio.h>#include <stdlib.h>#include <math.h>#define pi 3.14159

int

main()

{

double

R=6371.0,a1,b1,a2,b2,d,c;

scanf

(

"%lf%lf%lf%lf"

,&a1,&b1,&a2,&b2);

a1=a1/180*pi;

b1=b1/180*pi;

a2=a2/180*pi;

b2=b2/180*pi;

c=

cos

(b1)*

cos

(b2)*

cos

(a1-a2)+

sin

(b1)*

sin

(b2);

d=

acos

(c)*R;

printf

(

"%.6lf"

,d);

return

0;

}