haversine.c

/*
Recibe dos pares de coordenadas en grados (así como las da Google maps)
Regresa la distancia que hay entre esos dos puntos en Kilómetros
*/
double calcularDistanciaEntreDosCoordenadas(double lat1, double lon1, double lat2, double lon2)
{
    // Convertir todas las coordenadas a radianes
    lat1 = gradosARadianes(lat1);
    lon1 = gradosARadianes(lon1);
    lat2 = gradosARadianes(lat2);
    lon2 = gradosARadianes(lon2);
    // Aplicar fórmula
    double RADIO_TIERRA_EN_KILOMETROS = 6371;
    double diferenciaEntreLongitudes = (lon2 - lon1);
    double diferenciaEntreLatitudes = (lat2 - lat1);
    double a = pow(sin(diferenciaEntreLatitudes / 2.0), 2) + cos(lat1) * cos(lat2) * pow(sin(diferenciaEntreLongitudes / 2.0), 2);
    double c = 2 * atan2(sqrt(a), sqrt(1 - a));
    double d = RADIO_TIERRA_EN_KILOMETROS * c;
    return d;
}