détails

d: y = 0
F = (2.242848221255715, 1.6541124873013304)
K = Midpoint(F, Mirror[F, d])
p = Vector(F, K)
c = Parabola(F, d)
fy = OrthogonalLine(F, d)
P = (10.396357972236412, 10.671034552349159)
ax = OrthogonalLine(P, fy)
H = Intersect(fy, ax)
Ω = Translate(H, p)
HΩ = Vector(H, Ω)
f = OrthogonalLine(Ω, fy)
V = Mirror(Ω, H)
U = Mirror(V, P)
g = Segment(V, U)
h: y = (sgn(x(P) - x(H)) * (x - x(Ω))) + y(Ω)
M = Intersect(h, g)
i = OrthogonalLine(P, g)
N = Intersect(h, i)
I = Midpoint(M, N)
e = Circle(I, M)
k = Semicircle(Ω, I)
T = Intersect(e, k, 1)
q = Circle(Ω, T)
F'_1, F_1 = Intersect(q, h)
r = Hyperbola(F_1, F'_1, P)
P_2, P_3, P_1 = Intersect(c, r)
a = false
n_1 = OrthogonalLine(P_1, Tangent[P_1, c])
n_2 = OrthogonalLine(P_2, Tangent[P_2, c])
n_3 = OrthogonalLine(P_3, Tangent[P_3, c])
P_{4} = Point(c)
n_p = OrthogonalLine(P_{4}, Tangent[P_{4}, c])
j = Segment(P_{4}, F)
l = OrthogonalLine(F, j)
J = Intersect(n_p, l)
C = Mirror(P_{4}, J)
lieu1 = Locus(C, P_{4})
J' = Dilate(J, 1.5, F)
m = Segment(F, J')