dualmoo > euh mais la méthode habituelle pour trouver les orbites c'est de résoudre u en fonction de thêta, pas en fonction de t, si ? (de toute façon c'est ça qui te donne la trajectoire : r en fonction de t ça t'intéresse pas beaucoup pour ça ^^)

parce que si on cherche à résoudre directement en fonction de t, on a r'=-u'/u² et r"=2u'²/u³ - u"/u², donc euh le changement de variable est pas super génial