Computing the Critical Points of the Distance Function Between Two Keplerian Orbits via Rigorous Global Optimization