This procedure for computing analytically the orbital phase corresponds to what is called in the jargon the “Taylor T2
approximant”. We refer to Ref. [98] for discussions on the usefulness of defining several types of approximants for computing
(in general numerically) the orbital phase.