создание анимации движения луны и земли вокруг солнца

[VeryStupidPerson]

Я пытался написать анимацию движения луны вокруг земли и земли вокруг солнца, но абсолютно не понимаю, как заставить луну выписывать окружность относительно изменяющихся координат земли.
Пытался я сделать так:

import pygame, sys, time
from pygame.locals import *
import math
pygame.init()
FPS=100
fpsClock=pygame.time.Clock()
width=1500
height=700
x, y = 150, 320
xl, yl = 50, 320
mainSurface=pygame.display.set_mode((width,height),0,32)
pygame.display.set_caption('Keyb moves')
background=pygame.image.load('bg1.jpg')
zemla=pygame.image.load('pict2.gif')
solnze=pygame.image.load('pict1.gif')
luna=pygame.image.load('pict3.gif')
# game loop
a = 8
b = 4
deg = 0
while True:
    fpsClock.tick(FPS) # frame rate
    mainSurface.blit(background,(0,0))
    
    # get all events from the queue
    for event in pygame.event.get():
        # loop events queue
        if event.type == QUIT:
            # window close X pressed
            pygame.quit()
            sys.exit()
    mainSurface.blit(solnze, (650, 280))
    mainSurface.blit(zemla, (x, y))
    mainSurface.blit(luna, (xl, yl))
    rad = math.radians(deg)
    x = a * math.sin(rad)  + x
    y = -b * math.cos(rad) + y
    xl =  2*math.sin(rad) + x + xl # в принципе, вопрос только про эту
    yl = 2* math.cos(rad) + y + xl # и эту строки. (Ну, может быть, ещё про 2 предыдущие.)
    pygame.display.update()
    if deg < 361:
        deg += 1
    else:
        deg = 0

Но луна сразу улетает в неизвестном направлении.
Посоветуйте, пожалуйста, какими формулами описать движение луны относительно осей х и у. (я подозреваю, что это, собственно, не совсем относится к программированию, но на счёт этих формул в физике я информации не нашёл)

[type_Oleg]

Я старый солдат, и Питоном не владею. Поэтому скажу чисто математически.
Сразу скажу, что точное, честное движение Луны - вещь сложная, теория Брауна содержит неск. тысяч слагаемых в формуле. Но вам наверное не это надо, а упрощенное движение. То есть 2 движения по окружности с периодами 365.25 и 27.5 дней ( https://ru.wikipedia.org/wiki/Сидерический_период )
Все просто - сначала считаем гелиоцентрические координаты Земли - R*cos ... , R*sin ..
Потом - геоцентрические координаты Луны - R*cos ... , R*sin .. Конечно, другие R и другой период.
И наконец - складывая их, получаем гелиоцентрические координаты Луны.
Вот и все.
Только не ожидайте, что у вас получатся какие-то завитульки, если все сделать в правильном масштабе то получится вот что - https://youtu.be/SefK_RSrnfM?t=1200
С момента 20.00

[type_Oleg]

Если надо - вот файл Excel, разбирайтесь. Можете поэкспериментировать , поменять период и радиус лунной орбиты. Тогда могут получиться завитульки.

[VeryStupidPerson]

Спасибо. Всё заработало. Просто я очень затупил с программой. К формуле земли надо было прибавлять координаты солнца, а я прибавлял предыдущие координаты земли.