import numpy as np
import matplotlib.pyplot as plt

#Liste des demi-grand axe en unité astronomique
a_ua=[0.387,0.723,1,1.524,5.203,9.555,19.22,30.10]
#Liste des demi-grand axe en m à compléter
a=[]
#Liste des périodes en jours
T_jours=[87.97,224.7,365.26,365.26+322,11*365.26+315,29*365.26+167,84*365.26+27,164*365.26+280]
#Liste des périodes en secondes à compléter
T=[]

#Liste des T^2 à compléter
T2=[]

#Liste des a^3 à compléter
a3=[]

# Courbe a^3 en fonction de T^2 
plt.plot() # à compléter
plt.show()
# Régression linéaire
pente,ordao = np.polyfit() # à compléter
# Obtention de la masse du Soleil d'après la 3e loi de Kepler
print("masse Soleil=",,"kg") # à compléter et à comparer à la valeur tabulée : 1.9891e30 kg