We demonstrate the power of a recently-proposed approximation scheme for the non-perturbative renormalization group that gives access to correlation functions over their full momentum range. We solve numerically the leading-order flow equations obtai
ned within this scheme, and compute the two-point functions of the O(N) theories at criticality, in two and three dimensions. Excellent results are obtained for both universal and non-universal quantities at modest numerical cost.
