In a companion paper (hep-th/0512317), we have presented an approximation
scheme to solve the Non Perturbative Renormalization Group equations that
allows the calculation of the n-point functions for arbitrary values of the
external momenta. The method was applied in its leading order to the
calculation of the self-energy of the O(N) model in the critical regime. The
purpose of the present paper is to extend this study to the next-to-leading
order of the approximation scheme. This involves the calculation of the 4-point
function at leading order, where new features arise, related to the occurrence
of exceptional configurations of momenta in the flow equations. These require a
special treatment, inviting us to improve the straightforward iteration scheme
that we originally proposed. The final result for the self-energy at
next-to-leading order exhibits a remarkable improvement as compared to the
leading order calculation. This is demonstrated by the calculation of the shift
ΔTc, caused by weak interactions, in the temperature of Bose-Einstein
condensation. This quantity depends on the self-energy at all momentum scales
and can be used as a benchmark of the approximation. The improved
next-to-leading order calculation of the self-energy presented in this paper
leads to excellent agreement with lattice data and is within 4% of the exact
large N result.Comment: 35 pages, 11 figure