由经纬度计算距离

    前些天接触到一个由经纬度计算距离的算法，非常简明了哦。
    怎奈进行操作的时候，却报出来类型错误的结果。后来经过分析，发现问题出在下面代码中的13行。
     /* 计算两个GPS点直线距离 结果单位：M

 3     */
 4     public static double distance(double jingdu1, double weidu1,
 5             double jingdu2, double weidu2) {
 6         double theta = jingdu1 - jingdu2;
 7         double dist = Math.sin(deg2rad(weidu1)) * Math.sin(deg2rad(weidu2))
 8                 + Math.cos(deg2rad(weidu1)) * Math.cos(deg2rad(weidu2))
 9                 * Math.cos(deg2rad(theta));
10
11
12
13         dist = Math.acos(dist);
14         dist = rad2deg(dist);
15         return dist * 60 * 1.1515 * 1.609344 * 1000;
16     }
17     // 将角度转换为弧度
18     static double deg2rad(double degree) {
19         return degree / 180 * Math.PI;
20     }
21     
22     // 将弧度转换为角度
23     static double rad2deg(double radian) {
24         return radian * 180 / Math.PI;
25     }

        第7行得到的dist，要在第13行进行反余弦计算。但是在某个具体的经纬度时(概率特别特别小，很遗憾，没有记下来)，会得到一个1.00000000000000000001的值。  我们知道，余弦值最大只能为1啦，所以反余弦的时候就会处理为错误喽。所以在10、11、12行加上如下代码就ok啦！

 if(dist>1){ 
      dist = 1.0;
  }