a = 2.5 (slider from 0 to 5, step 0.1)
R_E = 1 (slider from 1 to 5, step 0.1)
R_P = 0.49 (slider from 0.1 to 0.5, step 0.01)
tps = 2.495 (slider from 0 to 10, step 0.0050)
P = 1.8 (slider from 1 to 6, step 0.1)
θ: 2 * π * tps / P
P_P: (a * cos(θ), 0)
d_{in} = If(abs(x(P_P)) < R_E, true, false)
c_P = Circle(P_P, R_P)
d_{droite} = If(x(P_P) > 0, true, false)
M: (If[d_{droite}, R_E, -R_E], 0)
L: (x(P_P) + If[d_{droite}, -R_P, R_P], 0)
x_P: x(P_P)
d_{ES} = If((abs(x_P) < R_E + R_P) ∧ (abs(x_P) > R_E - R_P), true, false)
O = (0, 0)
c_E = Circle(O, R_E)
J_1, J_2 = Intersect(c_E, c_P)
J: (x(J_1), abs(y(J_1)))
sct_E = CircleSector(O, If[d_{droite}, M, J], If[d_{droite}, J, M])
sct_P = CircleSector(P_P, If[d_{droite}, J, L], If[d_{droite}, L, J])
stri, o, p_p, j = Polygon(P_P, J, O)
sobt: sct_P + sct_E - stri
d_{phase} = If(sin(θ) < 0, true, false)
c_{orb} = Circle(O, a)
P': (a; θ)
scach = If(d_{ES}, 2 * sobt, If[d_{in}, π * R_P^(2), 0])
c_{P'} = Circle(P', R_P)
s_{vue}: π * R_E^(2) - scach * If[d_{phase}, 1, 0]
T: (x(P_P), s_{vue} / (π * R_E^(2)) * 3)
T': T
Trace = true
L': (x(P_P) + If[d_{droite}, R_P, -R_P], 0)