Among methods for finding numerical solutions of equations describing stationary two-dimensional flows of a viscous gas in a bounded domain, the majority deal with flows which are incompressible in the sense that the density can be considered constant in the whole domain; i.e., the velocity field is diver-genceless. For many problems of practical importance this simplification does not hold. In particular, there are convective problems in which a strong heat source generates nonnegligible temperature gradients throughout the flow while the pressure gradients are considerably smaller. From the equation of state it follows that the gas density has then to be treated as variable in the domain [1,2].