The present paper deals with the numerical computation of the viscous fingering phenomenon occurring when gas is pushed into a thin layer of liquid between two parallel flat plates. The phenomenon is formulated mathematically as a gas-liquid two-phase flow problem. The moving interface between gas and liquid is not treated as a part of boundaries of a solution domain, but is recognized as the surface of density discontinuity in the domain. The density is treated as one of primary field variables of the problem. The three-dimensional flow equations, namely the momentum, continuity and energy equations, are integrated and averaged in the direction of thickness of the layer to yield the two-dimensional governing equations expressed in terms of velocity, pressure and density. The equations thus derived are discretized in space by the finite element method and in time by the finite difference method. Three characteristic growth patterns called spreading, splitting and shielding and their combined phenomenon have been calculated. Encouraging results have been obtained.