Abstract
A novel method to compute all critical points of the distance function between two Keplerian orbits (either bounded or unbounded) with a common focus is presented. The problem is attacked as a global optimization problem, solved by a rigorous global optimizer based on Taylor models. Thus, thigh enclosures of the stationary points are obtained. The embedded capability of the method of delivering high-order Taylor expansions is then used to analyze how uncertain orbital parameters affect the position of the stationary points and the associated distance values. Sample orbital sets and Apophis asteroid are used as test cases.