Je réserve ici
Spécialité Physique-Chimie Terminale

Ressources · Terminale · Spécialité Physique-Chimie

Satellites et planètes : mouvement circulaire et lois de Kepler

Champ de gravitation d’un astre et choix du référentiel ; mouvement circulaire uniforme d’un satellite, vecteurs vitesse et accélération ; période de révolution et troisième loi de Kepler ; satellite géostationnaire ; les trois lois de Kepler sur une orbite elliptique ; données astronomiques et satellitaires testées en Python.

À propos de cette page
Dix-huit exercices corrigés de terminale, trois par cas, du plus simple au plus difficile : force du Soleil sur Mercure et point où les attractions de la Terre et de la Lune se compensent, accélération de la Lune et d’un module en orbite lunaire, satellites de navigation et satellites de Jupiter, erreur des 24 h et délai des communications par satellite, orbites de Pluton, d’une comète et de type Molniya, programmes Python sur des planètes et des satellites.

18 exercices, classés par cas comme le cours. Écris ta réponse, puis « Vérifier ».

Cas 1Champ de gravitation et référentiel

Réflexe — Je convertis l’altitude en distance au centre, r=R+h, en mètres ; je choisis le référentiel lié au centre de l’astre attracteur ; seul soumis à son attraction, un corps a l’accélération a→=G→. Revoir la méthode →

Exercice 1★★★Facile

Le Soleil a pour masse MS=1,99×1030 kg ; Mercure a pour masse m=3,30×1023 kg, et la distance moyenne entre leurs centres vaut d=5,79×1010 m. On prend G=6,67×10−11 m³·kg⁻¹·s⁻² et l’on ne tient compte que de l’attraction du Soleil.

  1. Référentiel adapté à l’étude du mouvement de Mercure : référentiel
    Mercure tourne autour du Soleil : on se place dans le référentiel héliocentrique (centre du Soleil, axes vers des étoiles lointaines), supposé galiléen.
  2. Valeur de la force exercée par le Soleil sur Mercure : F= × 10²² N
    F=GMSmd2=6,67×10−11×1,99×1030×3,30×1023(5,79×1010)2≈ 1,31 × 10²² N, dirigée de Mercure vers le centre du Soleil.
  3. Valeur du champ de gravitation créé par le Soleil à cette distance : G= × 10⁻² N·kg⁻¹
    G=GMSd2=1,33×10203,35×1021≈ 3,96 × 10⁻² N·kg⁻¹ ; on retrouve aussi Fm.
  4. Valeur de l’accélération de Mercure : a= × 10⁻² m·s⁻²
    Deuxième loi de Newton : ma→=mG→, donc a→=G→ et a= 3,96 × 10⁻² m·s⁻². La masse de Mercure s’est simplifiée : c’est la même valeur que le champ, et elle ne dépend pas de la planète.
Exercice 2★★★Moyen

Un satellite d’observation tourne à l’altitude h=700 km. Données : RT=6380 km ; GMT=3,98×1014 m³·s⁻².

  1. Distance du satellite au centre de la Terre : r= km
    r=RT+h=6380+700= 7 080 km, soit 7,08×106 m.
  2. Champ de gravitation terrestre à cette altitude : G= N·kg⁻¹
    G=GMTr2=3,98×1014(7,08×106)2≈ 7,94 N·kg⁻¹.
  3. Un élève prend r=700 km : il trouve G= N·kg⁻¹
    3,98×1014(7,00×105)2≈ 812 N·kg⁻¹, plus de 80 fois le champ au sol : un résultat absurde, qui signale l’oubli du rayon terrestre.
  4. Ce champ vaut % de sa valeur au sol, calculée avec le même modèle
    Au sol r=RT : le rapport vaut (RTr)2=(63807080)2≈0,812, soit 81,2 %. Ni G ni MT ne sont utiles pour un rapport.
  5. Altitude à laquelle le champ terrestre vaut le neuvième de sa valeur au sol : h= km
    Diviser G par 9 revient à multiplier r par 9=3 : r=3RT, donc h=r−RT=2RT= 12 760 km, et non 3RT.
Exercice 3★★★Difficile

Une sonde, moteurs coupés, traverse le segment qui joint les centres de la Terre et de la Lune. Données : MT=5,97×1024 kg ; ML=7,35×1022 kg ; distance entre les centres d=3,84×108 m ; GMT=3,98×1014 m³·s⁻² ; G=6,67×10−11 SI. On néglige l’attraction du Soleil.

  1. Rapport des masses MTML=
    5,97×10247,35×1022≈ 81,2.
  2. Distance au centre de la Terre du point où les deux champs ont la même valeur : × 10⁸ m
    À la distance x de la Terre : GMTx2=GML(d−x)2, donc xd−x=MTML≈9,01 et x=9,0110,01d≈ 3,46 × 10⁸ m, à environ 38 400 km du centre de la Lune.
  3. Valeur de chacun des deux champs en ce point : × 10⁻³ N·kg⁻¹
    GT=3,98×1014(3,456×108)2≈ 3,33 × 10⁻³ N·kg⁻¹, et le champ lunaire a la même valeur.
  4. Accélération de la sonde en ce point : a= m·s⁻²
    Les deux champs ont même valeur et des sens opposés (l’un vers la Terre, l’autre vers la Lune) : leur somme est nulle, donc a= 0. Ce n’est pas un point où la gravitation disparaît : deux attractions s’y compensent.
  5. Au-delà de ce point, en allant vers la Lune, l’accélération de la sonde est dirigée vers
    Plus près de la Lune, le champ lunaire l’emporte : l’accélération est dirigée vers la Lune, et une sonde qui a franchi ce point « tombe » vers elle.

Cas 2Vitesse et accélération d’un satellite en orbite circulaire

Réflexe — Système, référentiel, une seule force centripète ; la deuxième loi de Newton et le repère de Frenet donnent dvdt=0 et v2r=GMr2, d’où v=GMr. Revoir la méthode →

Exercice 4★★★Facile

Un satellite décrit une orbite circulaire de rayon r=8500 km autour du centre de la Terre, dans le référentiel géocentrique supposé galiléen. Donnée : GMT=3,98×1014 m³·s⁻².

  1. Le vecteur accélération du satellite est dirigé vers le de la Terre
    Seule force : l’attraction terrestre, dirigée vers le centre de la Terre ; d’après la deuxième loi de Newton, a→ a la même direction et le même sens : il est centripète.
  2. Coordonnée tangentielle de l’accélération : dvdt= m·s⁻²
    La force n’a pas de coordonnée selon u→t : dvdt= 0, la valeur de la vitesse est constante.
  3. Valeur de l’accélération : a= m·s⁻²
    a=GMTr2=3,98×1014(8,50×106)2≈ 5,51 m·s⁻².
  4. Valeur de la vitesse : v= km·s⁻¹
    v=GMTr=3,98×10148,50×106≈6,84×103 m·s⁻¹, soit 6,84 km·s⁻¹.
  5. Angle entre le vecteur vitesse et le rayon de l’orbite : °
    Le vecteur vitesse est tangent au cercle, donc perpendiculaire au rayon : 90°. Il est perpendiculaire au vecteur accélération.
Exercice 5★★★Moyen

La Lune décrit autour de la Terre une orbite considérée comme circulaire, de rayon r=3,84×108 m, en T=27,3 jours. Données : GMT=3,98×1014 m³·s⁻² ; 1 jour = 86 400 s.

  1. Valeur de la vitesse de la Lune, calculée à partir de r et de T : v= km·s⁻¹
    T=27,3×86400≈2,359×106 s et v=2πrT=2π×3,84×1082,359×106≈1,02×103 m·s⁻¹, soit 1,02 km·s⁻¹.
  2. Accélération déduite du mouvement, a=v2r=4π2rT2 : a= × 10⁻³ m·s⁻²
    a=4π2×3,84×108(2,359×106)2≈ 2,72 × 10⁻³ m·s⁻².
  3. Champ de gravitation terrestre à la distance de la Lune : G= × 10⁻³ N·kg⁻¹
    G=3,98×1014(3,84×108)2≈ 2,70 × 10⁻³ N·kg⁻¹.
  4. Écart relatif entre ces deux valeurs, arrondi à l’unité : %
    2,725−2,6992,699≈0,01, soit 1 % environ. La deuxième loi de Newton prévoit a=G : la mesure du mouvement la confirme ; le petit écart vient surtout de ce que la Terre est, elle aussi, attirée par la Lune : le centre de la Terre, origine du référentiel géocentrique, n’est pas tout à fait immobile, et la Lune tourne en réalité autour du centre de masse du système Terre-Lune ; l’attraction du Soleil réduit un peu cet effet.
  5. Un élève affirme : « la Lune tourne à vitesse constante, donc aucune force ne s’exerce sur elle ». Ce qui change en permanence, c’est la du vecteur vitesse
    La valeur de la vitesse est constante, mais sa direction tourne : le vecteur vitesse varie, l’accélération n’est pas nulle, et une force s’exerce, l’attraction terrestre, qui courbe la trajectoire sans changer la valeur de la vitesse.
Exercice 6★★★Difficile

Un module habité tourne autour de la Lune sur une orbite circulaire, à l’altitude h=100 km. Données : ML=7,35×1022 kg ; RL=1737 km ; G=6,67×10−11 m³·kg⁻¹·s⁻².

  1. Produit GML= × 10¹² m³·s⁻²
    GML=6,67×10−11×7,35×1022≈ 4,90 × 10¹² m³·s⁻², 81 fois moins que pour la Terre.
  2. Vitesse du module : v= km·s⁻¹
    r=1737+100=1837 km et v=4,90×10121,837×106≈1,63×103 m·s⁻¹, soit 1,63 km·s⁻¹.
  3. Un élève remplace r par l’altitude : il trouve v= km·s⁻¹
    4,90×10121,00×105≈ 7,00 km·s⁻¹ : plus de quatre fois trop, alors que l’orbite n’est qu’à 100 km du sol lunaire.
  4. Accélération d’un astronaute qui flotte dans le module : a= m·s⁻²
    L’astronaute est, comme le module, soumis à la seule attraction lunaire : a=G=GMLr2=4,90×1012(1,837×106)2≈ 1,45 m·s⁻². Il flotte parce que le module a la même accélération, pas parce que le champ serait nul.
  5. Altitude à laquelle un module en orbite circulaire autour de la Lune aurait une vitesse de 1,00 km·s⁻¹ : km
    De v2=GMLr : r=GMLv2=4,90×1012(1,00×103)2≈4902 km, donc h=r−RL≈ 3 165 km. Pour aller moins vite, il faut monter.

Cas 3Période de révolution et troisième loi de Kepler

Réflexe — J’écris T=2πrv ou T2r3=4π2GM en unités SI ; pour deux corps du même astre, T2=T1(r2r1)3/2. Revoir la méthode →

Exercice 7★★★Facile

Un satellite de navigation décrit une orbite circulaire de rayon r=26580 km autour du centre de la Terre. Données : GMT=3,98×1014 m³·s⁻² ; jour sidéral : 86 164 s.

  1. Vitesse du satellite : v= km·s⁻¹
    v=3,98×10142,658×107≈3,87×103 m·s⁻¹, soit 3,87 km·s⁻¹.
  2. Période de révolution : T= × 10⁴ s
    T=2πrv=2π×2,658×1073,87×103≈ 4,32 × 10⁴ s.
  3. Cette période vaut 11 h min
    431593600≈11,99 h ; 0,99×60≈59 : 11 h 59 min.
  4. Nombre de révolutions effectuées pendant un jour sidéral :
    8616443159≈ 2,00 : le satellite repasse au-dessus du même point du sol après deux tours, une fois par jour.
Exercice 8★★★Moyen

Io, satellite de Jupiter, décrit une orbite quasi circulaire de rayon r=4,22×108 m en T=1,77 jour. Ganymède, autre satellite de Jupiter, tourne à r′=1,07×109 m. Données : G=6,67×10−11 SI ; MT=5,97×1024 kg ; 1 jour = 86 400 s.

  1. Période d’Io en secondes : T= × 10⁵ s
    T=1,77×86400≈ 1,53 × 10⁵ s.
  2. Masse de Jupiter : MJ= × 10²⁷ kg
    MJ=4π2r3GT2=4π2×(4,22×108)36,67×10−11×(1,529×105)2≈ 1,90 × 10²⁷ kg.
  3. Jupiter est fois plus massive que la Terre (à l’unité près)
    1,902×10275,97×1024≈ 319 (318 avec la masse arrondie à 1,90 × 10²⁷ kg).
  4. Période de Ganymède : T′= jours
    Même astre central : T′=T(r′r)3/2=1,77×(1,07×1094,22×108)1,5≈ 7,15 jours, sans convertir en secondes : les unités se simplifient.
Exercice 9★★★Difficile

Un satellite d’observation de la Terre, sur une orbite circulaire, effectue exactement 14,5 révolutions par jour de 86 400 s. Données : GMT=3,98×1014 m³·s⁻² ; RT=6380 km.

  1. Période : T= s
    T=8640014,5≈ 5 959 s, soit 99,3 min.
  2. Rayon de l’orbite : r= km
    r=(GMTT24π2)1/3=(3,98×1014×595924π2)1/3≈7,10×106 m, soit 7 100 km. Attention à la racine cubique.
  3. Altitude : h= km
    h=r−RT=7100−6380= 720 km.
  4. Vitesse : v= km·s⁻¹
    v=2πrT=2π×7,100×1065959≈ 7,49 km·s⁻¹.
  5. Un second satellite effectue 16 révolutions par jour : son altitude vaut h′= km
    T′=5400 s : r′=r×(54005959)2/3≈6649 km, donc h′≈ 269 km. Plus de tours par jour, c’est une orbite plus basse et plus rapide.

Cas 4Le satellite géostationnaire

Réflexe — Orbite circulaire, dans le plan de l’équateur, d’ouest en est, de période T=86164 s ; le rayon vient de la troisième loi de Kepler. Revoir la méthode →

Exercice 10★★★Facile

Un satellite de télévision est géostationnaire. Le rayon de son orbite vaut r=42140 km.

  1. Durée d’un tour de la Terre sur elle-même dans le référentiel géocentrique : 86 164 s, soit 23 h min 4 s
    86164−23×3600=3364 s =56×60+4 : 23 h 56 min 4 s, le jour sidéral.
  2. Plan de l’orbite : plan de l’
    La force passe par le centre de la Terre, et un point du sol décrit un cercle centré sur l’axe des pôles : seule l’orbite contenue dans le plan de l’équateur permet de rester au-dessus du même point.
  3. Sens de parcours : celui de la rotation de la Terre, d’ouest en
    D’ouest en est, comme la rotation de la Terre : en sens inverse, le satellite défilerait dans le ciel.
  4. Angle dont tourne le satellite autour du centre de la Terre en 6,00 h : °
    Un tour, 360°, en 86 164 s : en 21 600 s, 360°×2160086164≈ 90,2°, exactement l’angle dont a tourné la Terre : le satellite reste au-dessus du même point.
Exercice 11★★★Moyen

Un élève calcule le rayon de l’orbite géostationnaire en prenant T=86400 s (24 h) au lieu du jour sidéral. Données : GMT=3,98×1014 m³·s⁻² ; le bon rayon, calculé avec 86 164 s, vaut 4,2143×107 m.

  1. Rayon obtenu par l’élève : r′= km
    r′=(3,98×1014×8640024π2)1/3≈4,2220×107 m, soit 42 220 km.
  2. Écart avec le bon rayon : km
    42219,9−42143,0≈ 77 km : la période n’a changé que de 0,27 %, mais l’erreur atteint plusieurs dizaines de kilomètres.
  3. Sur cette orbite, le satellite fait un tour en 86 400 s. Pendant que la Terre fait un tour (86 164 s), il prend un retard angulaire de °
    En 86 164 s, le satellite ne parcourt que 360°×8616486400≈359,02° : un retard de 0,98° par jour sidéral, environ 1°.
  4. Vu du sol, ce satellite dérive lentement vers l’
    Il tourne d’ouest en est moins vite que la Terre : par rapport au sol, il glisse vers l’ouest et sort en quelques semaines du faisceau des antennes fixes.
Exercice 12★★★Difficile

Un satellite géostationnaire relaie des communications. Son altitude vaut h=35760 km ; les signaux se propagent à c=3,00×108 m·s⁻¹. On considère des stations situées à la verticale du satellite.

  1. Durée du trajet aller-retour d’un signal station → satellite → station : s
    Δt=2hc=2×3,576×1073,00×108≈ 0,238 s.
  2. Lors d’une conversation relayée, délai minimal entre la fin d’une question et le début de la réponse entendue : s
    La question monte et redescend (2h), puis la réponse aussi (2h) : 4hc≈ 0,477 s, un silence que l’on remarque au téléphone.
  3. Un satellite décrit une orbite circulaire équatoriale de même rayon, mais d’est en ouest. Il repasse au-dessus du même point de l’équateur toutes les h (à 0,01 près)
    Même rayon, donc même période, 86 164 s, mais en sens inverse : 1S=1T+1Tsid=286164, donc S=43082 s, soit 11,97 h.
  4. Un satellite décrit d’ouest en est une orbite circulaire équatoriale de période 6,00 h. Il repasse au-dessus du même point de l’équateur toutes les h (à 0,01 près)
    1S=121600−186164, donc S≈28830 s, soit 8,01 h : plus que sa période, car le sol tourne dans le même sens que lui.

Cas 5Les lois de Kepler sur une orbite elliptique

Réflexe — Je calcule a=rp+ra2 ; j’applique la loi des aires aux deux sommets, rpvp=rava ; je calcule une période par T=a3/2 (années, ua) ou par T2a3=4π2GM (SI). Revoir la méthode →

Exercice 13★★★Facile

Pluton passe au périhélie à 29,7 ua du Soleil et à l’aphélie à 49,3 ua. Pour le système solaire, T2a3=1 en années et en unités astronomiques.

  1. Le centre du Soleil occupe un de l’ellipse décrite par Pluton
    Première loi de Kepler : le Soleil occupe l’un des deux foyers, pas le centre de l’ellipse.
  2. Demi-grand axe : a= ua
    a=rp+ra2=29,7+49,32= 39,5 ua.
  3. Période de révolution : T= ans
    T=a3/2=39,51,5≈ 248 ans.
  4. Rapport des vitesses vpva=
    Loi des aires aux deux sommets : vpva=rarp=49,329,7≈ 1,66.
Exercice 14★★★Moyen

Une comète décrit autour du Soleil une orbite elliptique de demi-grand axe a=17,8 ua ; elle passe au périhélie à rp=0,586 ua du Soleil. On travaille en années et en ua.

  1. Période de révolution : T= ans
    T=a3/2=17,81,5=17,8×17,8≈ 75,1 ans.
  2. Distance au Soleil à l’aphélie : ra= ua
    ra=2a−rp=35,6−0,586≈ 35,0 ua, au-delà de l’orbite de Neptune.
  3. Rapport des vitesses vpva=
    vpva=rarp=35,010,586≈ 59,8 : près de 60 fois plus vite au périhélie.
  4. Loi qui explique qu’on ne voie la comète que quelques mois tous les 75 ans : loi des
    La comète ne brille que près du Soleil, là où elle va le plus vite : la loi des aires montre qu’elle passe l’essentiel de sa période dans la partie lointaine et lente de l’orbite.
Exercice 15★★★Difficile

Un satellite de télécommunications suit une orbite très allongée de type « Molniya », de période T=43082 s (un demi-jour sidéral), dont le périgée est à 600 km d’altitude. Données : GMT=3,98×1014 m³·s⁻² ; RT=6380 km.

  1. Demi-grand axe : a= km
    Troisième loi avec le demi-grand axe : a=(GMTT24π2)1/3=(3,98×1014×4308224π2)1/3≈ 26 550 km.
  2. Distance du centre de la Terre à l’apogée : ra= km
    rp=6380+600=6980 km (distance au centre, pas altitude) ; ra=2a−rp=53097−6980≈ 46 120 km.
  3. Altitude de l’apogée : km
    46117−6380≈ 39 740 km, plus haut qu’un satellite géostationnaire.
  4. Rapport des vitesses au périgée et à l’apogée :
    vpva=rarp=461176980≈ 6,61 : le satellite passe très lentement à l’apogée, où il reste longtemps visible.
  5. Durée du trajet du périgée à l’apogée : h
    Une demi-orbite dure une demi-période : 430822=21541 s, soit 5,98 h.

Cas 6Tester les lois de Kepler avec Python

Réflexe — Troisième loi : je calcule T2a3 pour chaque corps ou je trace T2 en fonction de a3 ; deuxième loi : je compare les aires 12|x1y2−x2y1| balayées pendant des durées égales. Revoir la méthode →

Exercice 16★★★Facile

On exécute le programme suivant.

a = 1.524   # demi-grand axe de Mars (ua)
T = 1.881   # période de Mars (années)
k = T**2 / a**3
print(round(k, 2))

  1. Valeur affichée :
    1,88121,5243=3,5383,540≈0,9996, arrondi à 1.0 : Mars vérifie la troisième loi, avec la constante 1 du système solaire en années et en ua.
  2. En Python, a3 s’écrit a3
    La puissance s’écrit ** : a**3. Le symbole ^ désigne une autre opération et provoque ici une erreur.
  3. Pour un astéroïde fictif de demi-grand axe 4.00 ua, un élève écrit T = a**3/2 : valeur obtenue
    Python calcule d’abord a**3 =64, puis divise par 2 : 32, ce qui n’est pas la période.
  4. Période correcte de cet astéroïde : ans
    T=a3/2=4,001,5=43= 8,00 ans, ce qui s’écrit a**1.5 ou a**(3/2).
Exercice 17★★★Moyen

Une simulation donne les positions d’un astéroïde fictif autour du Soleil, placé à l’origine, en ua, relevées à 0,05 an d’intervalle. Près du périhélie : P1(1,000;0,000) puis P2(0,952;0,379). Près de l’aphélie : Q1(−3,000;0,000) puis Q2(−2,995;−0,128). Le programme calcule A = 0.5 * abs(x1*y2 - x2*y1).

  1. Aire du triangle Soleil-P1-P2 : ua²
    A1=12|1,000×0,379−0,952×0|≈ 0,190 ua².
  2. Aire du triangle Soleil-Q1-Q2 : ua²
    A2=12|(−3,000)×(−0,128)−(−2,995)×0|= 0,192 ua² : les deux aires sont égales à 1 % près (le triangle remplace l’arc), la loi des aires est vérifiée.
  3. Distance P1P2 parcourue en 0,05 an : ua
    P1P2=(0,952−1,000)2+0,3792≈ 0,382 ua.
  4. Distance Q1Q2 : ua
    Q1Q2=0,0052+0,1282≈ 0,128 ua.
  5. Rapport des vitesses moyennes près du périhélie et près de l’aphélie :
    0,3820,128≈ 2,98, proche de rarp=3,0001,000=3 : la relation rpvp=rava est retrouvée.
Exercice 18★★★Difficile

Un programme exploite les orbites, quasi circulaires, de trois satellites de la Terre : un satellite Galileo, un satellite GLONASS et un satellite géostationnaire. Données : G=6,67×10−11 SI.

# Trois satellites de la Terre
# r en km, T en s
r = [29600, 25510, 42164]
T = [50685, 40548, 86164]
for i in range(3):
    k = T[i]**2 / (r[i]*1e3)**3
    print(k)

  1. Valeur affichée pour le satellite Galileo : × 10⁻¹⁴ s²·m⁻³
    506852(2,96×107)3=2,569×1092,593×1022≈ 9,91 × 10⁻¹⁴ s²·m⁻³.
  2. Valeur affichée pour le satellite GLONASS : × 10⁻¹⁴ s²·m⁻³
    405482(2,551×107)3≈ 9,90 × 10⁻¹⁴ s²·m⁻³ : la même valeur à 0,02 % près, et de même pour le satellite géostationnaire.
  3. Masse de la Terre déduite de la valeur commune k des trois satellites, arrondie à trois chiffres : × 10²⁴ kg
    Les trois valeurs valent k≈9,90×10−14 s²·m⁻³, d’où MT=4π2Gk=39,486,67×10−11×9,90×10−14≈ 5,98 × 10²⁴ kg, la masse de la Terre à 0,2 % près.
  4. Un élève oublie le facteur 1e3 : la masse de la Terre qu’il déduit de la valeur affichée vaut × 10¹⁵ kg
    Avec r en kilomètres, r3 est 109 fois plus petit : k est multiplié par 109, donc M=4π2Gk est divisée par 109 : 5,98 × 10¹⁵ kg, une masse absurde pour une planète. Les trois valeurs affichées restent pourtant égales entre elles : la loi semble vérifiée, seule la masse trahit l’erreur d’unité.
Bloqué sur ce chapitre ?

Moi c'est Ben, j'ai créé ce site pour aider le plus d'élèves possible, partout en France. Je peux aussi reprendre avec toi ce qui coince en cours particulier : en visio, ou à domicile si tu es à Marseille.

Prof de physique-chimie à Marseille · Cours particuliers au lycée · Aide aux devoirs