The principal approximation methods used to compute the inspiral of compact binary systems are the post-Newtonian (pN) expansion, in which an orbital angular velocity MΩ serves as the expansion parameter; and the self-force or extreme-mass-ratio-inspiral approach, in which the small parameter is the mass ratio m/M of the binary’s two components. We work in an overlapping regime where both approximations are valid and find numerical values of pN coefficients at orders beyond the reach of current analytical work. In this talk we present a novel analytic extraction of high-order pN parameters that govern quasi-circular binary systems using ultra-high accuracy numerical computations.