A multi-phase pseudospectral method is formulated for aerospace trajectory optimization with interior point constraints (IPCs). A verification technique for the associated covector mapping theorem is then developed and the "initial guess" obtained by this mapping is also used to start the indirect method.