C = Point(xAxis)
B = Point(yAxis)
C' = Mirror(C, yAxis)
B' = Mirror(B, xAxis)
e = Ellipse(C, C', B)
A', A = Intersect(e, xAxis)
P = Point(e)
P' = Point(e)
r = Line(P, P')
S = Intersect(yAxis, r)
c = Circle(S, C, C')
G, D = Intersect(c, r)
luogo1 = Locus(D, P')