C = Point(xAxis)
C' = Mirror(C, yAxis)
B = Point(yAxis)
e = Ellipse(C, C', B)
B' = Mirror(B, xAxis)
A', A = Intersect(e, xAxis)
P = Point(e)
d = Circle(P, 0.001)
P', F = Intersect(d, e)
r = Line(P, P')
T = Intersect(yAxis, r)
c = Circle(T, C, C')
a = OrthogonalLine(P, r)
N_1, N = Intersect(c, yAxis)