Derivation of N-burn analytic solutions for propellant-optimal transfer trajectories of a vehicle in a vacuum between arbitrary boundary conditions. Variational changes in the desired boundary conditions are expressed, in general, in terms of variational changes in the control vector and in the initial state vector. All coefficient matrices are computed recursively in terms of the analytic matrices established from the subarcs of the N-burn solution. The solution is applicable to shuttle ascent (exoatmospheric), rendezvous, and deorbit problems. Consideration is also given to state-variable and control-variable inequality constraints.