Soleil

per_V = 224.7
a_V = 0.7233
tps = 21.5 (slider from 0 to 4000, step 0.1)
θ_V: tps * 360.0 * ° / per_V
per_T = 365.25
θ_T: tps * 360.0 * ° / per_T
a_T = 1
V: (a_V; θ_V)
T: (a_T; θ_T)
S = (0, 0)
cv = Circle(V, 0.4)
r_0 = 0.2
x_V = 0
y_V = 2
d_{TV} = Distance(T, V)
r_V: r_0 / d_{TV}
c_V = Circle((x_V, y_V), r_V)
ua = 1.4959787E8
R_{Venus} = 6056
φ = Angle(Vector[Vector[V, S]], Vector[Vector[V, T]])
VT = Segment(V, T)
VS = Segment(V, S)
ST = Segment(S, T)
dp_{SV} = OrthogonalLine(V, VS)
I_1, I_2 = Intersect(cv, dp_{SV})
ii = Segment(I_1, I_2)
J_1, J_2 = Intersect(cv, OrthogonalLine[V, VT])
jj = Segment(J_1, J_2)
S' = Intersect(cv, Ray[V, Vector[S, V]])
sct = CircumcircleSector(I_1, S', I_2)
L = If(Distance[I_1, T] < Distance[I_2, T], I_1, I_2)
P = Intersect(jj, OrthogonalLine[L, jj])
sgt_{LP} = Segment(L, P)
M = If(Angle[T, V, J_1] > 180.0 * °, J_2, J_1)
α = Angle(L, V, M)
xsol = If(φ > 180.0 * °, 1.5, -1.5)
C_S: (xsol + x_V, y_V)
c_{sol} = Circle(C_S, 0.15)
psol = If(Angle[V, T, S] > 180.0 * °, 1.5, -1.5)
P': (x_V + r_V * cos(α), y_V)
a_e: r_V
b_e: r_V * cos(α)
elim = CurveCartesian(b_e * cos(t) + x_V, a_e * sin(t) + y_V, t, -1.5707963267948966, 1.5707963267948966)
co_T = Circle(S, a_T)
co_V = Circle(S, a_V)
φ': φ + 180.0 * °
aa: φ' > 0.0 ∧ φ' < 1.01 * 3.141592653589793
texte3: "Phase : " + φ + " " + (If[φ > 89.5 * ° ∧ φ < 90.5 * °, "Plus longue élongation Ouest", If[φ > 269.5 * ° ∧ φ < 270.5 * °, "Plus longue élongation Est", If[φ > 179.5 * ° ∧ φ < 180.5 * °, "Conjonction inférieure", If[φ' > 179.5 * ° ∧ φ' < 180.5 * °, "Conjonction supérieure", ""]]]])
texte1: "diam. Vénus : " + (2.0 * R_{Venus} / (d_{TV} * ua) * 180.0 * 3600.0 / 3.141592653589793) + "''" + "
dist T-V : " + d_{TV} + " u.a."
dc_V = If(xsol > 0.0, Semicircle[(x_V, y_V + r_V), (x_V, y_V - r_V)], Semicircle[(x_V, y_V - r_V), (x_V, y_V + r_V)])