  General relativity predicts the gravitational wave signatures of coalescing
binary black holes. Explicit waveform predictions for such systems, required
for optimal analysis of observational data, have so far been achieved using the
post-Newtonian (PN) approximation. The quality of this treatment is unclear,
however, for the important late-inspiral portion. We derive late-inspiral
waveforms via a complementary approach, direct numerical simulation of
Einstein's equations. We compare waveform phasing from simulations of the last
$\sim 14$ cycles of gravitational radiation from equal-mass, nonspinning black
holes with the corresponding 2.5PN, 3PN, and 3.5PN orbital phasing. We find
phasing agreement consistent with internal error estimates based on either
approach, suggesting that PN waveforms for this system are effective until the
last orbit prior to final merger.
