В работе исследуется динамика квантовых систем, взаимодействующих с электромагнитным полем. Вне рамок резонансного приближения и приближения вращающейся волны статистическая матрица плотности и вероятности квантовых переходов представлены интегралами по траекториям в пространстве энергетических состояний исследуемой системы.На базе полученных выражений разработан алгоритм численного моделирования динамики квантовых систем, взаимодействующих с электромагнитным полем.