Вообще-то нужно Вас отослать к файлу, в котором все написано. Но так уж и быть повторю выкладки. Во первых нельзя переходить в систему, связанную с центром инерции тела, так как решение должно получиться комплексное. Кроме того, решение для двух тел, которое решают второкурсники не распространяется на множество тел, а предлагаемый метод распространяется. Решение в комплексной плоскости отличается от решения, которое реализуют второкурсники и оно определяет дискретные уровни энергии.
Имеем уравнение движения для каждого тела
Используя комплексные полюса знаменателя которые непрерывны, получим уравнение относительно каждой переменной.
где величина
определяется по формуле
Причем эта величина в ноль и бесконечность не обращается и следовательно является монотонной функцией времени. Если эту величину подставить в формулу (1), то получим уравнение движения.
Формулу (1) можно представить в виде
как это сделать, см. в прилагаемом файле. Уравнение (2) можно представить в виде
Умножаем (3) на величину
и интегрируем по величине
, получим
Далее извлекаем корень из левой части и интегрируем это дифференциальное уравнение ., получим
Получим решение уравнения движения, причем координата зависит от времени и от целого числа p. Считаем интеграл энергии и получаем зависимость энергии от целого числа p. Т.е. получен дискретный спектр энергии.