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)
a = Segment(F_1, P)
b = Segment(F_2, P)
n, t = AngularBisector(a, b)
E = Point(n)
α = Angle(F_2, P, E)
β = Angle(E, P, F_1)
F = Point(n)
G = Point(t)
γ = Angle(G, P, F)
δ = Angle(F_2, P, F_1)