Actualmente estoy trabajando en un algoritmo para calcular las posiciones planetarias utilizando la ley de Kepler. Mientras lo probaba con el servicio de efemérides suministrado por https://ssd.jpl.nasa.gov/horizons.cgi La anomalía real es correcta, pero la distancia se desvía en un 1%, así que me pregunto si hay algún problema con la distancia suministrada por el servicio o si he calculado la distancia de forma incorrecta.
Los valores se refieren a Mercurio el 1 de enero de 2000, a las 00:00. Para obtener los elementos orbitales, se introducen los siguientes parámetros:
Tipo de efemérides: ELEMENTOS
Cuerpo objetivo: Mercurio
Origen de coordenadas: 10 (es el ID del cuerpo del sol, que será el centro de coordenadas)
Rango de tiempo: Start=2000-01-01 00:00, Stop=2000-01-01 01:00, Step=1
Configuración de la tabla: por defecto
Visualización/salida: por defecto
El resultado muestra que a la hora indicada (1 de enero de 2000 00:00), los elementos eran los siguientes:
verdadera anomalía: 1.751155303115542E+02 (degrees)
eje semimayor: 3,870982252717257E-01 (AU)
excentricidad: 2.056302512089075E-01
Si se utilizan estos tres valores para calcular la distancia $r$ utilizando la fórmula $$r=\frac{a(1-e^2)}{1 + e\cdot\cos(v)}$$ donde $a$ es el semieje mayor, $e$ es la excentricidad y $v$ es la verdadera anomalía, se obtiene un resultado de aproximadamente 0,466 UA.
Pero si en cambio se solicita la distancia del lugar cambiando el tipo de efemérides de ELEMENTOS a VECTORES, se obtiene la distancia (indicada por el valor RG) a 0,47 UA.
Estoy muy confundido con esto, y espero que alguien pueda arrojar algo de luz sobre este misterio.