Background. Presentation of probability as an intrinsic property of the nature leads researchers to switch from deterministic to stochastic description of the phe- nomena. On the basis of the ideology of N. G. van Kampen and C. W. Gardiner the procedure of stochastization of one-step process was formulated. It allows to write down the master equation based on the type of of the kinetic equations (equations of interactions) and assumptions about the nature of the process (which may not neces- sarily process of birth–death). The kinetics of the interaction has recently attracted attention because it often occurs in the physical, chemical, technical, biological, envi- ronmental, economic, and sociological systems. However, there are no general methods for the direct study of this equation. The expansion of the equation in a formal Taylor series (so called Kramers–Moyal’s expansion) is used in the procedure of stochasti- zation of one-step processes. It is also possible to apply system size expansion (van Kampen’s expansion). Leaving in the expansion terms up to second order we can get the Fokker–Planck equation, and after the Langevin equation. It should be clearly understood that these equations are approximate recording of the master equation. Purpose. However, this does not remove the need for the study of the master equation. Moreover, that the power series produced during the decomposition of the master equation may be divergent (for example, in spatial models). This makes it impossible to apply the classical perturbation theory. Method. It is proposed to use quantum field perturbation theory for the statistical systems (so-called method Doi). The perturbation series are treated in the spirit of the Feynman path integral, where the Green’s functions of the perturbed Liouville operator of the master equation are propagators. For more convenience of selection of the perturbed and unperturbed parts of the Liouville operator and to obtain the explicit form of the Green function of the master equation we need to rewrite the equation in the occupation number representation (Fock state). Results. This work is a methodological material that describes the principles of master equation solution based on methods of quantum field perturbation theory. The feature presentation is that it is designed for non-specialists in quantum field theory. This example uses Verhulst model because of its simplicity and clarity (the first order equation, is independent of the spatial variables, however, contains non-linearity). Conclusions. The described method allows to solve directly the master equation and obtain the model parameters, such as the moments and the probability density distribution.