A Sixth-Order Accurate Scheme for Solving Two-Point Boundary Value Problems in Astrodynamics