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)
M = Point(c)
a = OrthogonalLine(M, b)
e = Segment(M, O)
P' = Intersect(t, e, 1)
f = OrthogonalLine(P', d)
P = Intersect(f, a)
g = OrthogonalLine(M, e)
I = Intersect(b, g)
j = Line(P, I)
H = Intersect(b, a)
l = Segment(M, H)