已知地球的平均半径为6371.393千米,假设在地球的某一纬度上,有两个处于不同经度的点A、B,用C语言编程求出这两点之间的直线距离、这两点在该纬度剖面上切圆的劣弧长度,以及这两点之间的球面距离。
以下是我写的代码。
#include "stdio.h"
#include "math.h"
main()
{
double x,y,z,t,s,d,l,R,r=6378.137,pi=3.1415926;
scanf("%lf%lf%lf",&x,&y,&z);
if((t=fabs(z-y))>180)
t=360-t;
t=2*pi*t/360;
x=fabs(x);
x=2*pi*x/360;
R=r*cos(x);
l=2*R*sin(t/2);
s=R*t;
d=r*t;
printf("%f,%f,%f\n",l,s,d);
system("pause");
return 0;
}
老师给的参考数据是
北纬0度,东经0度,东经90度——直线距离=9010.510271,纬度切圆劣弧长度=10008.160550,球面距离=10008.160550
而我求出来的是,9020.047727,10018.754000,10018.754000
老师要求尽可能高的精确度
我是哪里出了问题