c++实现根据地图上两点经纬度计算两点间实际距离

#include <iostream>
#include <math.h>


using namespace std;


const double EARTH_RADIUS = 6371.004;//地球半径


double rad(double d)
{
   return d * M_PI / 180.0;
}


double GetDistance(double lat1, double lng1, double lat2, double lng2)
{
   double radLat1 = rad(lat1);
   double radLat2 = rad(lat2);
   double a = radLat1 - radLat2;
   double b = rad(lng1) - rad(lng2);


   double s = 2 * asin(sqrt(pow(sin(a/2),2) +
    cos(radLat1)*cos(radLat2)*pow(sin(b/2),2)));
   s = s * EARTH_RADIUS;
   //s = (s * 10000) / 10000;
   return s;
}


int _tmain(int argc, _TCHAR* argv[])
{
double dist = GetDistance(39.97274, 116.38649, 39.99274,116.39649);