El JPL, que calcula las órbitas de los cuerpos celestes con gran precisión, utiliza una expresión para la aceleración de un cuerpo de masa despreciable debida a la fuerza gravitatoria de otro cuerpo que se parece a la siguiente:
$\frac{d\bar{v}}{dt}=-\frac{GM}{r^2}{(1-\frac{4GM}{rc^2}}+\frac{v^2}{c^2})\hat{r}+ \frac{4GM}{r^2}(\bar{v}\cdot{\hat{r}})\frac{v^2}{c^2}\hat{r}$
Ignorando las partes dependientes de la velocidad tenemos: $\frac{d\bar{v}}{dt}=-\frac{GM}{r^2}\hat{r}+\frac{4(GM)^2}{c^2r^3}\hat{r}$
Así que en realidad mantener $a = 2$ pero añadir una pequeña parte "cúbica inversa" se hace en realidad para ajustarse mejor a los valores experimentales, aunque el término cúbico inverso no es inventado sino que proviene de intentar aproximar los efectos relativistas generales utilizando lo que se conoce como "la expansión post-newtoniana".
La razón de esto es un poco complicada, pero básicamente añadiendo una pequeña parte cúbica inversa así como partes dependientes de la velocidad se explica lo que se conoce como "precesión anómala del perihelio".
Véase, por ejemplo, la expresión 4-61 en este documento titulada "Formulación de los valores observados y calculados de los tipos de datos de la red de espacio profundo para la navegación"