Probablemente hay muchas formas de hacer esto.
Ten en cuenta que la ecuación paramétrica de la esfera es muy similar a la ecuación paramétrica de un esferoide. Si puedes calcular la geodésica de una esfera a partir de su representación paramétrica, es posible que puedas aplicar el mismo procedimiento al esferoide.
http://www.cs.iastate.edu/~cs577/handouts/geodesics.pdf
También podrías utilizar las herramientas del cálculo tensorial calculando la Matriz de Transformación de Coordenadas Esféricas a Oblatas $[\frac{dx^a}{d\bar{x}^u}]$, y luego calcular la métrica oblata usando $\bar{g}_{uv} = \frac{dx^a}{d\bar{x}^u} \frac{dx^b}{d\bar{x}^u}g_{ab}$ y luego aplicar la ecuación geodésica a la métrica de las Coordenadas Oblatas y resolverla.
Creo que también es posible calcular la geodésica en Coordenadas Esféricas y transformar tu respuesta a Coordenadas Oblatas debido a que la solución será una ecuación vectorial. Obviamente, sería mucho más fácil.
Debido a que las representaciones paramétricas son muy similares, las matrices de transformación son fáciles de calcular.
También está este método directo: http://mathworld.wolfram.com/OblateSpheroidGeodesic.html