已知两点经纬度求两点的大地方位角

代码说明：

应用背景物理大地测量学中，经常用到两点间的方位角，已知球面两点大地经纬度坐标，求解两点间的球面方位角。假设地球为球形进行的计算。关键技术通过两点的大地经纬度坐标，计算两点球面坐标方位角。需要输入文件INPUT，输入格式 纬度 经度 大地高（或正高），可以是一组数据也可以是连续的多个数据，排列顺序是按列三列数据。计算出的方位角是以第一个点到第二个点的方位角。例如，假如序号为1,2,3,4……。那么求出的方位角为α1-2，α2-3，α3-4……。输出文件为fangweijiao.dat