Given two bodies on Keplerian orbits, what is the procedure to find the next closest approach of the two bodies?  The bodies would be effectively massless objects like two spaceships orbiting around a uniform central body.  The orbits can be given by state vectors or the six keplerian elements or any representation of elements that leads to the easiest solution.

I expect it is a numerical root finding problem followed by a second derivative check of some sort.