The paper presents a mathematical model, and its solution, of non-isothermal radial gas flow in porous media. The system of partial differential equations has been converted into two ordinary differential equations and solved numerically. In the case of the gas condensate system the properties of the liquid and gas phases were evaluated using the Peng-Robinson equation of state. The condensation process depends on the simultaneous decrease of pressure and temperature. These two parameters are intrinsically conjugated, and it is not possible to separate each from other. The equilibrium model was solved by the Quasi-Newton Successive Substitution (QNSS) and Dominant Eigenvalue Method (DEM). A new function describing the saturation of drop out condensate for the non-equilibrium process has been proposed.