Solving the Lambert Problem

This function solves the Lambert orbital boundary value problem by determining the velocity vectors at two specified position vectors given the transfer time between them. The implementation handles both short-way and long-way trajectories with numerical iteration methods to compute orbital parameters, though the approach may produce high delta-v requirements.