O = Intersect(xAxis, yAxis)
B = Point(yAxis)
A = Point(yAxis)
c = Circle(O, A)
h = OrthogonalLine(B, d)
C = Intersect(c, h, 1)
i = OrthogonalLine(C, h)
F' = Intersect(b, i)
F = Mirror(F', O)
k = Ellipse(F', F, B)
t = Circle(O, B)
D = (1.0094, 1.4122)
I = Point(b)
a = Line(B, D)
J = Intersect(b, a)
e = Line(A, J)
Δ = Line(D, I)
g = OrthogonalLine(D, b)
K = Intersect(b, g)
D' = Intersect(g, e)
Δ' = Line(I, D')
M', M = Intersect(c, Δ')
l = OrthogonalLine(M', b)
m = OrthogonalLine(M, b)
P = Intersect(Δ, m)
P' = Intersect(Δ, l)
H = Intersect(b, m)
H' = Intersect(b, l)
n = Segment(B, D)
p = Segment(B, J)
q = Segment(A, D')
r = Segment(A, J)
f = Segment(M, H)
j = Segment(M', H')
s = Segment(D', K)