O = Intersect(xAxis, yAxis)
A_{1} = Point(xAxis)
eq1: x^(2) / (x(A_{1}))^(2) + y^(2) = 1
a: x(A_{1})
B_{1}, C_{1} = Root(((-a^(2)) * x^(2)) + a^(4) + (2 * a^(3)) - ((2 * a) * x^(2)) - x^(2))
eq2: x = x(B_{1})
D = Intersect(eq1, eq2, 1)
f = Tangent(D, eq1)
E = Intersect(f, yAxis)
F = Intersect(eq1, yAxis, 1)
g = OrthogonalLine(F, yAxis)
c = Circle(O, E)
G = Intersect(c, f, 2)
h = Segment(E, G)
H = Intersect(c, g, 1)
i = Segment(E, H)
j = Segment(H, G)
A = Point(c)
k, l = Tangent(A, eq1)
B = Intersect(c, k, 2)
m = Segment(A, B)
C = Intersect(c, l, 2)
n = Segment(B, C)
p = Segment(C, A)
P = Intersect(eq1, n, 1)
Q = Intersect(eq1, p, 1)
R = Intersect(eq1, m, 1)
YJT = 1 (slider from 0 to (2 * π))