Bien le bonjour à tous.
Dans le cadre de nos cours nous devons réaliser la "Mise en équation et résolution numérique des équations différentielles régissant le mouvement des planètes et de leur satellite. Calcul de la trajectoire optimale d'un engin spatial."
Nous avons "simplifié" le problème en ne prenant que le système soleil-terre-lune.
Nous avons alors trois équations différentielles du second ordre, avec trois variables, positions du soleil, de la terre et de la lune :
Soleil: msr''s(t) = ( msmt (rs-rt) / |rs-rt|3 ) + ( msml (rs-rl) / |rs-rl|3 )
Terre: mtr''t(t) = ( mtms (rt-rs) / |rt-rs|3 ) + ( mtml (rt-rl) / |rt-rl|3 )
Lune: mlr''l(t) = ( mtms (rt-rs) / |rt-rs|3 ) + ( mtml (rt-rl) / |rt-rl|3 )
Hélas, une fois dans mathematica, impossible de les résoudre...
Alors mes questions sont les suivantes :
1) Les équations sont-elles justes ?
2) Est-il impossible de les résoudre ou ont-elles été mal entrées dans mathematica ?
3) Y a-t-il une autre approche plus facile ?
Merci d'avance !
-----