Рассмотрена разработка компьютерной программы для квантово-механических расчетов вариационных квантовых систем нескольких частиц. Создано и в настоящее время поддерживается хранилище нового кода. Применен новый вычислительный метод, что позволяет точно рассчитать среднее значение некоторых сингулярных операторов, которые появляются в расширении релятивистской энергии малых атомов в условиях постоянной тонкой структуры. Произведены подробные численные расчеты на низком спектре атомов бора и лития. Такой результат был достигнут впервые для пяти-электронной системы.