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