The problem of trajectory design requires the determination of spacecraft position and velocity as a function of time that satisfy design constraints. The constraints that must be satisfied are supplied to the trajectory designer as parameters that are generally functions of the Cartesian state. Thus, the main interest in developing solutions of the equations of motion for navigation is to enable computation of parameters that satisfy mission constraints and state vectors that may be used to initialize numerical integration for further refinement of the trajectory design.