Simulación de 3 cuerpos y precisión de Lambert Interplanetary Solver

Descubrí que usar un solucionador Lambert junto con las ecuaciones de la aproximación cónica parcheada para un vuelo interplanetario de la Tierra a Marte es bastante inexacto, cuando se tiene en cuenta el campo gravitatorio de la Tierra. Supongo que esto es de esperar ya que es una aproximación , pero esto me hizo pensar: ¿Cómo se calcula con precisión una trayectoria interplanetaria entre dos planetas si se tienen en cuenta los efectos perturbadores del objetivo y, lo que es más importante, del planeta fuente? Supongo que una forma de arreglar la trayectoria una vez fuera de la esfera de influencia de la Tierra sería ejecutar el solucionador de Lambert nuevamente y aplicar el requeridodelta-va la nave espacial, pero esto parece algo muy ineficiente cuando el combustible es limitado. ¿Existen métodos más avanzados que el enfoque cónico parcheado de Lambert + y, de ser así, alguien tiene enlaces a artículos escritos sobre estos temas? Cualquier ayuda sería muy apreciada.

PD Este es el resultado que obtengo de una solución de Lambert al ignorar los efectos gravitacionales de la Tierra y Marte (el semicírculo interior es la órbita de la Tierra, el semicírculo exterior es la órbita de Marte y la línea que une los dos es la nave espacial) trayectoria):

Y este es el resultado que obtengo de una solución cónica parcheada de Lambert + cuando se toman en cuenta los efectos gravitacionales de la Tierra y Marte: