In this paper, we examine the energy distribution function of electrons in the case of a very weakly ionized
argon plasma at sub-atmospheric pressure and external electric field using Boltzmann kinetic equation. We are
interested in the behavior of the collisional part of the equation, thus spatially uniform model is considered. The
goal of the research is to compare two different numerical approaches: a deterministic one (using a two-term local
non-stationary approximation) and a stochastic approach (using the Monte Carlo method) over a wide range of
reduced electric fields: from several Td to kTd. We present the comparison for steady-state and time-dependent
solutions, isotropic and anisotropic parts of the electron energy distribution function, and reaction constants.
The research will also help to identify any limitations and challenges of these methods.