An optimal h6 scheme for solving TPBVP in astrodynamics