A method for the calculation of the nonlocal electron distribution function (EDF) with programs commonly used for solving the one-dimensional Boltzmann kinetic equation (e.g., COMSOL Multiphysics) was proposed and implemented. The capabilities of the proposed method were illustrated using the example of the positive column plasma in argon. Significant differences between the local and nonlocal EDFs were observed, especially at the plasma periphery. This can result in significant differences in the electron transport coefficients and reaction rate constants for electron-induced processes. The proposed approach increases the precision of the numerical modeling of gas discharge devices and equipment.