Application of High Order Expansions of Two-Point Boundary Value Problems to Astrodynamics