The influence of electron-electron collisions on the formation of a nonlocal electron distribution function (EDF) and other characteristics of the electron gas are analyzed. Correct expressions for the coefficients in the integral of electron-electron collisions are obtained that are suitable for substitution into the kinetic Boltzmann equation averaged over the volume. This method of EDF calculation is implemented in COMSOL Multiphysics software. The example of positive column plasma in argon shows that accounting for electron-electron collisions leads both to differences between nonlocal and local EDFs and to dependence on the ionization degree.