A further development of the Kostyuk-Pikelner's model is presented. The response of the chromosphere heated by non-thermal electrons of the power-law energy spectrum has been studied on the basis of the numerical solution of the one-dimensional time-dependent equations of gravitational gas dynamics. The ionization and energy loss for the emissions in the Lyman and Balmer lines have been determined separately for the optically thin and thick Lα-line layers. Due to the initial heating, a higher-pressure region is formed. From this region, disturbances propagate upwards (a shock wave with a velocity of more than 1000 km s-1) and downwards. A temperature jump propagates downwards, and a shock is formed in front of the thermal wave. During a period of several seconds after the beginning of this process, the temperature jump intensifies the downward shock wave and the large radiative loss gives rise to the high density jump (ϱ2/ϱ1 ∼ 100). The numerical solution has been analyzed in detail for the case heating of the ionized and neutral plasma, and a value of this heating is close to the upper limit of the admissible values. In this case, the condensation located between the temperature jump and the shock wave front, may emit in the observed optical continuum.