Erbites des exoplanètes
Author: phm

Excentricité

Masse planète (en masses Terre)

Période (en jours)

Masse étoile (en masses solaires)

Excentricité

P_{SYS} = 56.9 (slider from 1 to 400, step 0.1)
e = 0.37 (slider from 0 to 1, step 0.0050)
M_E = 1 (slider from 0.1 to 20, step 0.1)
M_P = 11550.1 (slider from 0.1 to 33000, step 0.1)
a: cbrt((P_{SYS} / 365.25)^(2))
a_E: a * M_P / (M_E * 330000)
a_P: a - a_E
F2_P: -2 * a_P * e
F2_E: 2 * a_E * e
el_P = Ellipse((0, 0), (F2_P, 0), a_P)
el_S = Ellipse((0, 0), (F2_E, 0), a_E)
tps = 17 (slider from 0 to 500, step 0.1)
M = Mod(tps * 360 / P_{SYS}, 360)
u_0: 3.141592653589793 / 180 * M
u_1: u_0 + (u_0 + e * sin(u_0) - u_0) / (1 - e * cos(u_0))
u_2: u_1 + (u_0 + e * sin(u_1) - u_1) / (1 - e * cos(u_1))
u_3: u_2 + (u_0 + e * sin(u_2) - u_2) / (1 - e * cos(u_2))
u_4: u_3 + (u_0 + e * sin(u_3) - u_3) / (1 - e * cos(u_3))
u: u_4 + (u_0 + e * sin(u_4) - u_4) / (1 - e * cos(u_4))
v: 2 * atan(tan(u / 2) * sqrt((1 + e) / (1 - e))) * 180 * ° / 3.14159
ρ: a_P * (1 - e^(2)) / (1 + e * cos(v))
P: (ρ; v)
d1 = Line(P, (0, 0))
S_1, S_2 = Intersect(el_S, d1)
S: S_2
sp = Segment(S, P)