Relevance.Forthe time being, when calculating the flow rate of a vertical gas well, a trinomial formula is used, which discounts the fact that at cost-efficient flow rates there are two simultaneous zones around a vertical gas well, where each zone has its own Darcy’s law. The extents of these zones depend both on the formation’s permeability/porosity parameters and on formation’s draw-down pressure. The aim of the research is todetermine a relationship between the flow rate of a vertical gas well and the formation’s draw-down pressure, as well as the formation’s permeability/porosity parameters, taking into account the extents of areas of linear and nonlinear gas flows. Objects: gas formations and perfect gas well, laws of linear and non-linear gas flow. Methods: analytical methods of subsurface hydrodynamics. Results. Currently, to calculate a flow rate of a vertical gas well, the well-known trinomial formula is used. This formula implies that gas flow is governed by nonlinear law throughout the entire flow area starting from the well through to the external reservoir boundary. Whereas in reality, the law of nonlinear flow only occurs in a narrow area around the well.The extent of this area depends on the flow rate of the well, formation parameters, gas properties, etc.The use of analytical subsurface hydrodynamics methods enables to find formulas and an algorithm for calculating the flow rate of a gas well, as well as to determine flow coefficients, accounting for the extents of the non-linear flow areas.The radius of the nonlinear flow area can be determined from the derived equation taking into account the critical flow rate calculated based on the formation’s permeability/porosity parameters. Equations of flow coefficients have been derived for the trinomial formula of the gas inflow to the well. The derived formula was used to calculate gas well production rate with various values ​​of formation’s draw-down pressure.The difference in calculation results of flow rates obtained using traditional versus the derived formula decreases after the draw-down pressure growth due to the increase in the extent of the non-linear flow area.Indicator charts are constructed according to the obtained and traditional formulas.Inflow performance relationship curves were plotted based on the derived and traditional formulas.Compared to the traditional trinomial formula, the equations for the A and B flow coefficients are changed, which now also account for the extents of the non-linear flow area.