请详细说明haversine函数的用法

`haversine`函数是用于计算地球上两个点之间的大圆距离的函数，其原理是根据两点的经纬度计算它们之间的弧长，再将弧长转换为实际距离。函数原型如下：

double haversine(double lat1, double lon1, double lat2, double lon2);

其中，`lat1`和`lon1`是第一个点的纬度和经度，`lat2`和`lon2`是第二个点的纬度和经度。这些参数的单位是弧度，而不是角度，因此如果你已经有了以角度为单位的经纬度坐标，需要先将其转换为弧度。可以使用下面的公式进行转换：

double deg2rad(double deg) {
    return deg * M_PI / 180.0;
}

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

`haversine`函数返回的是两个点之间的距离，单位通常是千米或英里，具体取决于你输入的经纬度是使用公制单位还是英制单位。下面是一个示例代码，用于计算纽约和旧金山之间的距离：

#include <cmath>
#include <iostream>

double deg2rad(double deg) {
    return deg * M_PI / 180.0;
}

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

double haversine(double lat1, double lon1, double lat2, double lon2) {
    double dlat = deg2rad(lat2 - lat1);
    double dlon = deg2rad(lon2 - lon1);
    double a = sin(dlat / 2) * sin(dlat / 2) + cos(deg2rad(lat1)) * cos(deg2rad(lat2)) * sin(dlon / 2) * sin(dlon / 2);
    double c = 2 * atan2(sqrt(a), sqrt(1 - a));
    double d = 6371 * c;  // 地球半径为6371公里
    return d;
}

int main() {
    double lat1 = 40.7128;  // 纽约的纬度
    double lon1 = -74.0060;  // 纽约的经度
    double lat2 = 37.7749;  // 旧金山的纬度
    double lon2 = -122.4194;  // 旧金山的经度
    double distance = haversine(deg2rad(lat1), deg2rad(lon1), deg2rad(lat2), deg2rad(lon2));
    std::cout << "The distance between New York City and San Francisco is " << distance << " km." << std::endl;
    return 0;
}

输出结果为：

The distance between New York City and San Francisco is 4133.11 km.

需要注意的是，由于地球并不是完全的球形，因此使用`haversine`函数计算的距离只是一个近似值，可能与实际距离存在一定的误差。