我有一个实现该算法的sql Server函数。这应该很容易翻译成其他语言：

create function dbo.udf_Haversine(@lat1 float,@long1 float,@lat2 float,@long2 float) returns float begin
    declare @dlon float,@dlat float,@rlat1 float,@rlat2 float,@rlong1 float,@rlong2 float,@a float,@c float,@R float,@d float,@DtoR float

    select @DtoR = 0.017453293
    select @R = 3959      -- Earth radius

    select 
        @rlat1 = @lat1 * @DtoR,@rlong1 = @long1 * @DtoR,@rlat2 = @lat2 * @DtoR,@rlong2 = @long2 * @DtoR

    select 
        @dlon = @rlong1 - @rlong2,@dlat = @rlat1 - @rlat2

    select @a = power(sin(@dlat/2),2) + cos(@rlat1) * 
                     cos(@rlat2) * power(sin(@dlon/2),2)
    select @c = 2 * atn2(sqrt(@a),sqrt(1-@a))
    select @d = @R * @c

    return @d 
end

这返回了里程的距离。对于公里，更换地球半径与公里相当。

Here是一个更深入的解释。

编辑：该功能足够快，足够准确，可以使用邮政编码数据库进行半径搜索。在this site多年来一直在做得很好(但现在不再存在链接了)。