При построении численной теории физической либрации перед нами встала задача сравнения численной и аналитической теорий орбитального движения Луны. В ходе исследования нам было важно понять, как различия в точности описания движения центра масс скажутся на поведении либрации. На данный момент мы представляем результаты первого этапа сравнения, а именно: описание поведения во времени различий между двумя типами орбитальных теорий. Нами была использована аналитическая теория Гутцвиллера и Шмидта, построенная в рамках главной проблемы, и современная численная теория DE431 JPL NASA, учитывающая большое количество факторов, выходящих за рамки главной проблемы, которые было бы сложно или даже невозможно учесть при получении аналитического решения. Произведя необходимые редукции при приведении обеих теорий к единой системе отсчета, мы сравнили оба решения на интервале в 800 лет. В итоге мы получили, что амплитуда в долготе за этот период не превосходит 80 угловых секунд, а в широте 10 угловых секунд. Основным источником расхождений являются эффекты, не учитываемые в главной проблеме, такие как: планетные возмущения, движение плоскости эклиптики, сжатие Земли, приливные эффекты и эффекты общей теории относительности.

The problem of comparing the numerical and analytical theories of the orbital motion of the Moon came before us at the stage of constructing a numerical theory of physical libration. In the course of the study, it was important for us to understand how the differences in the accuracy of the description of the center of mass motion will affect the behavior of the libration. Now we present the results of the first stage of the comparison, namely: the description of the time behavior of differences between the two types of orbital theories. We work with the analytical theory by Gutzwiller and Schmidt built in the framework of the main problem, and the modern numerical theory DE431 by JPL NASA, which takes into account a large number of factors that go beyond the main problem and which would be difficult or even impossible to take into account when obtaining an analytical solution. Having made the necessary reductions when we brought both theories to a unified reference frame, we compared both solutions on an interval of 800 years. As a result, we have obtained that the amplitude in longitude for this period does not exceed 80 arcseconds, and in latitude 10 arcseconds. The main sources of discrepancies are perturbations from planets and asteroids, motion of the ecliptic plane, the Earth's figure, tidal effects, and general relativity in numerical ephemerides, which were not included in analytical theory.