00001
00002
00003
00004
00005
00006 #include <stdio.h>
00007 #include "geodesic.h"
00008
00009 #if defined(_MSC_VER)
00010
00011 # pragma warning (disable: 4996)
00012 #endif
00013
00014
00015
00016
00017
00018
00019
00020
00021 int main() {
00022 double a = 6378137, f = 1/298.257223563;
00023 double lat1, lon1, azi1, lat2, lon2, azi2, s12;
00024 struct geod_geodesic g;
00025
00026 geod_init(&g, a, f);
00027 while (scanf("%lf %lf %lf %lf\n", &lat1, &lon1, &azi1, &s12) == 4) {
00028 geod_direct(&g, lat1, lon1, azi1, s12, &lat2, &lon2, &azi2);
00029 printf("%.15f %.15f %.15f\n", lat2, lon2, azi2);
00030 }
00031 return 0;
00032 }