  We reproduce the two-body gravitational conservative dynamics at third
post-Newtonian order for spin-less sources by using the effective field theory
methods for the gravitationally bound two-body system, proposed by Goldberger
and Rothstein. This result has been obtained by automatizing the computation of
Feynman amplitudes within a Mathematica algorithm, paving the way for
higher-order computations not yet performed by traditional methods.
