This paper is devoted to the description of computational algorithms for modeling quasi-one-dimensional non-steady flows of a multicomponent reacting gas. The particularity of the developed modeling technique is that the paths of strong and weak discontinuities are mobile computational nodes, and the parameters for them are calculated using special algorithms. A set of programs has been developed, which can be used to solve the problems of the reacting gas dynamics that are of applied importance, as well as serve as an illustrator for physical gas dynamics training courses. The paper provides the results of the numerical modeling of the supersonic flow in a flat channel simulating the operation of experimental facilities of the Institute for Problems in Mechanics and the Institute of Physics and Technology. A satisfactory correlation between the calculated and experimental data has been obtained.