F_1 = Point(xAxis)
B = Point(xAxis)
c = Circle(F_1, B)
F_2 = Point(xAxis)
D = Point(xAxis)
d = Circle(F_2, D)
Q = Point(c)
e = Circle(Q, Segment[F_2, D])
a = Line(F_1, Q)
F, G = Intersect(e, a)
b = Line(G, F_2)
M = Midpoint(G, F_2)
f = OrthogonalLine(M, b)
P = Intersect(a, f)
luogo1 = Locus(P, Q)
g = Circle(P, Q)
h = Segment(P, F_2)