(D)

P

T

e = 0.7 (slider from 0 to 5, step 0.1)
c: x^(2) + y^(2) = (e^(2) * (x + 1)^(2))
F = Intersect(xAxis, yAxis)
f = Circle(F, 1)
C = Intersect(d, a)
H = Point(d)
b = Line(H, F)
I = Dilate(H, e / (1 + e), F)
J = Dilate(H, e / (e - 1), F)
m = Midpoint(I, J)
g = Circle(m, I)
h = OrthogonalLine(H, d)
i = OrthogonalLine(m, h)
P_1, P_2 = Intersect(g, h)
K = Intersect(h, i)
P = Dilate(P_1, sgn(1 - e), K)
P' = Mirror(P, K)
j = Segment(P, m)
k = OrthogonalLine(P, j)
T = Intersect(d, k)
p = Circle(H, P, F)
l = Segment(P, F)
n = Segment(F, T)
α = Angle(P, F, T)
q = OrthogonalLine(I, b)
P_3 = Intersect(q, h)
T_3 = Intersect(q, d)
r = Circle(H, P_3, F)
s = Segment(P_3, F)
t = Segment(F, T_3)