A numerical method for the analysis of the fast axial flow glow discharge CO2 laser has been developed. The method is based on the self-consistent solution to the 1-D steady-state glow discharge equations, the gas dynamic equations, and the vibrational relaxation equations. The discharge equations include the continuity ones for the electrons, the positive and negative ions, and Poisson's equation for the electric field. The three-mode relaxation model for the vibrational kinetics and the plane-parallel optical resonator model have been used. This approach does not require previous assignment of the discharge power distribution and enables obtaining the discharge structure including the near-electrode regions in addition to the laser characteristics.