  In this paper, we propose a discretization of the multi-dimensional
stationary compressible Navier-Stokes equations combining finite element and
finite volume techniques. As the mesh size tends to 0, the numerical solutions
are shown to converge (up to a subsequence) towards a weak solution of the
continuous problem for ideal gas pressure laws p($\rho$) = a$\rho$ $\gamma$ ,
with $\gamma$ > 3/2 in the three-dimensional case. It is the first convergence
result for a numerical method with adiabatic exponents $\gamma$ less than 3
when the space dimension is three. The present convergence result can be seen
as a discrete counterpart of the construction of weak solutions established by
P.-L. Lions and by S. Novo, A. Novotn{\'y}.
