В данной методической работе в высокотемпературном пределе систематически выведено выражение для матрицы плотности для системы частиц, взаимодействующих посредством кулоновского потенциала. Выкладки соответствуют оригинальной работе Кельбга \cite{Kelbg:Annalen:1963}, в которой решается уравнение Блоха в первом порядке теории возмущений. Особое внимание было уделено объяснению всех преобразований при выводе, чтобы упростить понимание этой нетривиальной теории. Решение Кельбга широко используется в моделировании методом Монте-Карло с интегралами по траекториям.

In this paper the density matrix for a system of particles interacting via the Coulomb potential is systematically derived in the high–temperature limit. The derivation corresponds to the original work by Kelbg [\href{https://doi.org/10.1002/andp.19634670308}{Kelbg G 1963 Annalen der Physik \textbf{467} 219-224}], which solves the Bl\"och equation in the first order of perturbation theory. Special attention is paid to explain all the transformations in the derivation in order to simplify the understanding of this non-trivial theory. The solution of Kelbg is widely used in path integral simulations of Coulomb systems.