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
00022
00023 int main() {
00024 double a = 6378137, f = 1/298.257223563;
00025 double lat, lon, A, P;
00026 int n;
00027 struct geod_geodesic g;
00028 struct geod_polygon p;
00029 geod_init(&g, a, f);
00030 geod_polygon_init(&p, 0);
00031
00032 while (scanf("%lf %lf\n", &lat, &lon) == 2)
00033 geod_polygon_addpoint(&g, &p, lat, lon);
00034 n = geod_polygon_compute(&g, &p, 0, 1, &A, &P);
00035 printf("%d %.8f %.3f\n", n, P, A);
00036 return 0;
00037 }