A = Point(xAxis)
B = Mirror(A, yAxis)
C = Point(yAxis)
D = Mirror(C, xAxis)
F_1: (sqrt(x(A)^(2.0) - y(C)^(2.0)), 0.0)
F_2 = Mirror(F_1, yAxis)
c = Ellipse(F_1, F_2, C)
P = Point(c)
l_1 = Segment(F_1, P)
l_2 = Segment(F_2, P)
a, b = AngularBisector(l_1, l_2)