A = Point(xAxis)
B = Point(yAxis)
c: y(B)^(2.0) * x^(2.0) + x(A)^(2.0) * y^(2.0) = y(B)^(2.0) * x(A)^(2.0)
F_1: (sqrt(x(A)^(2.0) - y(B)^(2.0)), 0.0)
F_2: (-sqrt(x(A)^(2.0) - y(B)^(2.0)), 0.0)
P = Point(c)
a = Tangent(P, c)
b = OrthogonalLine(P, a)
u = Vector(F_2, P)
C = Point(b)
α = Angle(F_2, P, C)
C' = Rotate(C, α, P)
β = Angle(C, P, C')
d = Ray(P, F_1)
D, E = Intersect(c, d)
v = Vector(P, E)