Orthogonal Polynomial Approximation in Higher Dimensions: Applications in Astrodynamics