We present the first numerical application of a method that we have recently proposed to solve the Non Perturbative Renormalization Group equations and obtain the n-point functions for arbitrary external momenta. This method leads to flow equations for the n-point functions which are also differential equations with respect to a constant background field. This makes them, a priori, difficult to solve. However, we demonstrate in this paper that, within a simple approximation which turns out to be quite accurate, the solution of these flow equations is not more complicated than that of the flow equations obtained in the derivative expansion. Thus, with a numerical effort comparable to that involved in the derivative expansion, we can get the full momentum dependence of the n-point functions. The method is applied, in its leading order, to the calculation of the self-energy in a 3-dimensional scalar field theory, at criticality. Accurate results are obtained over the entire range of momenta.