Cosmological Post-Newtonian Expansions to Arbitrary Order

We prove the existence of a large class of one parameter families of solutions to the Einstein-Euler equations that depend on the singular parameter epsilon = nu(T)/c (0 < epsilon < epsilon(0)), where c is the speed of light, and nu(T) is a typical speed of the gravitating fluid. These solutions are shown to exist on a common spacetime slab M congruent to [0, T) x T(3), and converge as epsilon SE arrow 0 to a solution of the cosmological Poisson-Euler equations of Newtonian gravity. Moreover, we establish that these solutions can be expanded in the parameter epsilon to any specified order with expansion coefficients that satisfy epsilon-independent (nonlocal) symmetric hyperbolic equations.