tps = 4000 (slider from 0 to 4000, step 0.041666666666666664)
per_V = 583.92
Θ_V: tps * 360 * ° / per_V
a_V = 0.7233
S = (0, 0)
r_0 = 0.2
x_V = 0
y_V = 2
co_V = Circle(S, a_V)
V: (a_V; Θ_V)
T = (0, -1)
c_V = Circle(V, 0.4)
Vap: (x_V, y_V)
d_{VT} = Distance(V, T)
r_{Vap}: r_0 / d_{VT}
c_{Vap} = Circle(Vap, r_{Vap})
dp_{VS} = OrthogonalLine(V, Segment[V, S])
I_1, I_2 = Intersect(c_V, dp_{VS})
ii = Segment(I_1, I_2)
S' = Intersect(c_V, Ray[V, Vector[S, V]])
scv = CircumcircleArc(I_1, S', I_2)
dp_{VT} = OrthogonalLine(V, Segment[V, T])
J_1, J_2 = Intersect(c_V, dp_{VT})
jj = Segment(J_1, J_2)
L = If(Distance[I_1, T] > Distance[I_2, T], I_2, I_1)
P = Intersect(jj, OrthogonalLine[L, jj])
sgt_{LP} = Segment(L, P)
M = If(Angle[T, V, J_1] > 180°, J_2, J_1)
α = Angle(L, V, M)
P': (x_V + r_{Vap} * cos(α), y_V)
a_e: r_{Vap}
b_e: r_{Vap} * cos(α)
elim = CurveCartesian(b_e * cos(t) + x_V, a_e * sin(t) + y_V, t, -1.5707963267948966, 1.5707963267948966)