double  hypot( double  x,  double  y)  ... {
      return sqrt(x * x + y * y);
}

double  distance( double  wd1,  double  jd1,  double  wd2,  double  jd2)  ... {// 根据经纬度坐标计算实际距离
      double x, y, out;
      double PI = 3.1415926535898;
      double R = 6.371229 * 1e6;
       x = (jd2 - jd1) * PI * R * cos( ( (wd1 + wd2) / 2) * PI / 180) / 180;
       y = (wd2 - wd1) * PI * R / 180;
       out = hypot(x, y);
       return out;
  }