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

Используя комплексные полюса знаменателя которые непрерывны, получим уравнение относительно каждой переменной.
![$\frac{d^2 x_l}{dt^2}=\exp[H_l(t)]\frac{\exp[\int \ln(x_l-a(\alpha))d\alpha]}{\exp[\int \ln(x_l-c(\alpha,\beta))^2d\alpha d\beta]}\eqno(1)$ $\frac{d^2 x_l}{dt^2}=\exp[H_l(t)]\frac{\exp[\int \ln(x_l-a(\alpha))d\alpha]}{\exp[\int \ln(x_l-c(\alpha,\beta))^2d\alpha d\beta]}\eqno(1)$](https://dxdy-02.korotkov.co.uk/f/d/1/a/d1a3b6789d3a1627f1b78c9b9255efb882.png)
где величина
![$\exp[H_l(t)]$ $\exp[H_l(t)]$](https://dxdy-04.korotkov.co.uk/f/3/1/9/31936e4dc5fcafa59705876daf2f01c382.png)
определяется по формуле
![$\exp[H_l(t)]=-\frac{\alpha(y_l-x_l)}{|\vec y-\vec x|^3}\frac{\exp[\int \ln(x_l-c(\alpha,\beta))^2d\alpha d\beta]}{\exp[\int \ln(x_l-a(\alpha))d\alpha]}$ $\exp[H_l(t)]=-\frac{\alpha(y_l-x_l)}{|\vec y-\vec x|^3}\frac{\exp[\int \ln(x_l-c(\alpha,\beta))^2d\alpha d\beta]}{\exp[\int \ln(x_l-a(\alpha))d\alpha]}$](https://dxdy-02.korotkov.co.uk/f/5/7/2/572733f9d59fbe13e1cef5ac957904e782.png)
Причем эта величина в ноль и бесконечность не обращается и следовательно является монотонной функцией времени. Если эту величину подставить в формулу (1), то получим уравнение движения.
Формулу (1) можно представить в виде
![$\frac{d^2 x_l}{dh_l^2}=\frac{\exp[\int \ln(x_l-a(\alpha))d\alpha]}{\exp[\int 2\ln(x_l-c(\alpha,\beta))d\alpha d\beta]}\eqno(2)$ $\frac{d^2 x_l}{dh_l^2}=\frac{\exp[\int \ln(x_l-a(\alpha))d\alpha]}{\exp[\int 2\ln(x_l-c(\alpha,\beta))d\alpha d\beta]}\eqno(2)$](https://dxdy-04.korotkov.co.uk/f/3/8/5/385aea5cbd761c2b4682e1a1c185988182.png)
как это сделать, см. в прилагаемом файле. Уравнение (2) можно представить в виде

Умножаем (3) на величину

и интегрируем по величине

, получим
![$(\frac{d x_l}{dh_l})^2/2=P(x_l)+\int k(\alpha,\beta)[\ln(x_l-c(\alpha,\beta))+2\pi i p]d\alpha d\beta-\int \frac{b(\alpha,\beta)}{2(x_l-c(\alpha,\beta))}d\alpha d\beta\eqno(3)$ $(\frac{d x_l}{dh_l})^2/2=P(x_l)+\int k(\alpha,\beta)[\ln(x_l-c(\alpha,\beta))+2\pi i p]d\alpha d\beta-\int \frac{b(\alpha,\beta)}{2(x_l-c(\alpha,\beta))}d\alpha d\beta\eqno(3)$](https://dxdy-04.korotkov.co.uk/f/7/6/2/7622d96feaade244c3594390e1af938682.png)
Далее извлекаем корень из левой части и интегрируем это дифференциальное уравнение ., получим
![$\int_{u=x_l^0}^{x_l}\frac{du}{[P(u)+ \int k(\alpha,\beta)[\ln(u-c(\alpha,\beta))+2\pi i p]d\alpha d\beta-\int \frac{b(\alpha,\beta)}{2(u-c(\alpha,\beta))}d\alpha d\beta]^{1/2}}=h_l-h_l^0$ $\int_{u=x_l^0}^{x_l}\frac{du}{[P(u)+ \int k(\alpha,\beta)[\ln(u-c(\alpha,\beta))+2\pi i p]d\alpha d\beta-\int \frac{b(\alpha,\beta)}{2(u-c(\alpha,\beta))}d\alpha d\beta]^{1/2}}=h_l-h_l^0$](https://dxdy-02.korotkov.co.uk/f/d/5/9/d59bcb678d2a53327e3c3cc7aa9c31a382.png)
Получим решение уравнения движения, причем координата зависит от времени и от целого числа p. Считаем интеграл энергии и получаем зависимость энергии от целого числа p. Т.е. получен дискретный спектр энергии.