## Modules
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt

## Constantes
# Masse de la Lune
m = 7.36e22  # kg
# Constante de gravitation
G = 6.67e-11  # N m²/kg²
# Masse de la Terre
MT = 5.97e24  # kg

# Conditions initiales
theta0 = 0  # rad
dtheta0 = ...  # rad/s
r0 = 384e6  # m
dr0 = ...  # m/s

## Fonctions

def Fc(r: float) -> float:
    """Force centrale."""

    return ...

def F(X, t):
    ... 
    return


## Résolution numérique

# Durée de la simulation
T = ...  # s
# Nombre de points
N = ...
# Liste des instants
t = ...



## Tracé de la trajectoire


plt.figure()
plt.axis('equal')
plt.plot(...)
plt.plot(0, 0, 'ko', markersize = 10)
plt.xlabel('$x$ (m)')
plt.ylabel('$y$ (m)')
plt.grid()
plt.show()