Estou montando um pequeno programa para obter duas coordenadas geográficas de um usuário e depois calcular a distância entre elas (levando em consideração a curvatura da Terra). Então, procurei na Wikipedia qual é a fórmulaaqu.

Basicamente, configurei minha função python com base nisso e é isso que eu criei:

def geocalc(start_lat, start_long, end_lat, end_long):
    start_lat = math.radians(start_lat)
    start_long = math.radians(start_long)
    end_lat = math.radians(end_long)
    end_long = math.radians(end_long)

    d_lat = start_lat - end_lat
    d_long = start_long - end_long

    EARTH_R = 6372.8

    c = math.atan((math.sqrt( (math.cos(end_lat)*d_long)**2 +( (math.cos(start_lat)*math.sin(end_lat)) - (math.sin(start_lat)*math.cos(end_lat)*math.cos(d_long)))**2)) / ((math.sin(start_lat)*math.sin(end_lat)) + (math.cos(start_lat)*math.cos(end_lat)*math.cos(d_long))) )

    return EARTH_R*c

O problema é que os resultados são realmente imprecisos. Eu sou novo em python, então alguma ajuda ou conselho seria muito apreciado!