Costate Estimation for a Multiple-Interval Pseudospectral Method Using Collocation at the Flipped Legendre-Gauss-Radau Points