求纬度线上任意两点的球面距离

已知地球的平均半径为6371.393千米，并已知位于同一纬度上两点的纬度、经度值（单位为度），编写程序计算这两点之间的弦长和球面距离分别是多少千米。要求所有的实数使用double类型，圆周率的近似值取3.141592653589793。 提示：可以用正负来区分北纬与南纬、东经与西经。即北纬为正，南纬为负；东经为正，西经为负。

输入格式:

在第一行中输入一个位于区间[-90,90]内代表纬度的实数；在第二行中输入两个位于区间[-180,180]内代表经度的实数，以空格分隔。

输出格式:

输出弦长和球面距离的值，保留3位小数，以空格分隔，最后换行。

输入样例:

32.6
-56.1 163.8

输出样例:

10090.980 11645.760

#include <stdio.h>
#include <math.h>
#define pi 3.141592653589793
#define R 6371.393
int main()
{
    double c1,c2;
    double w,j1,j2,angle,r,a,y;
    scanf("%lf%lf%lf",&w,&j1,&j2);
    angle=j1-j2;
    if(angle<0)
    {
        angle=-angle;
    }
    if(angle>180)
    {
        angle=360-angle;
    }
    if(w<0)
    {
        w=-w;
    }
    w=(w/180)*pi;
    angle=(angle/180)*pi;
    r=R*cos(w);
    if(r<0)
    {
        r=-r;
    }
    c1=sqrt(2*r*r*(1-cos(angle)));
    y=(2*R*R-c1*c1)/(2*R*R);
    a=acos(y);
    c2=a*R;
    printf("%.3lf %.3lf\n",c1,c2);
    return 0;
}