Solve Lambert's Problem in Two-Body Dynamics
Solve the boundary value problem in two-body dynamics (known as Lambert's problem) using the universal variable formulation as presented in Bate, Mueller, and White. The function has been tested for elliptic, parabolic, and hyperbolic orbits.
The function does not currently handle the case where the secondary completes multiple revolutions about the central body within the specified time.
Inputs:
: Standard gravitational parameter of central body
,: Time scalar and column position vector for initial position in orbit
,: Time scalar and column position vector for final position in orbit
: If equal to 1, solve for the long path (delta_nu >= pi), else for the short path (delta_nu < pi)
Outputs:
,: Velocity vectors at initial and final times
: Number of iterations needed for convergence
Example:<