def METHOD_NAME(point_1, point_2):
    # lat lon to radians for haversine
    point_1 = [radians(_) for _ in point_1]
    point_2 = [radians(_) for _ in point_2]
    result = haversine_distances([point_1, point_2])
    # convert to km
    result *= 6371000 / 1000
    # result is a 2d distance matrix,
    #  0, dist
    #  dist, 0
    return result[0][1]