Binary black hole systems in the pre-coalescence stage are numerically constructed by demanding that the associated spacetime admits a helical Killing vector. Comparison with third order post-Newtonian calculations indicates a rather good agreement until the innermost stable circular orbit.