A = (2.7664, 1.8184)
B = (1.8664, 3.7784)
C = (5.2064, 2.6184)
ab = Line(B, A)
I = Intersect(d, ab)
A' = Mirror(A, d)
a'b = Line(A', B)
H = Midpoint(A', A)
bk = OrthogonalLine(B, d)
K = Intersect(d, bk)
BK = Segment(K, B)
J = Intersect(d, a'b)
mij = LineBisector(I, J)
O = Intersect(ab, mij)
co = Circle(O, I)
cm = OrthogonalLine(C, d)
M = Intersect(cm, d)
bc = Line(C, B)
I' = Intersect(d, bc)
C' = Mirror(C, d)
bc' = Line(C', B)
J' = Intersect(bc', d)
mij' = LineBisector(I', J')
O' = Intersect(bc, mij')
co' = Circle(O', I')
F, F' = Intersect(co', co)
AH = Segment(A, H)
A'H = Segment(A', H)
CM = Segment(M, C)
C'M = Segment(C', M)
ac = Line(C, A)
I'' = Intersect(d, ac)
a'c = Line(C, A')
J'' = Intersect(a'c, d)
mji'' = LineBisector(I'', J'')
O'' = Intersect(ac, mji'')
co'' = Circle(O'', I'')
FA = Segment(F, A)
FB = Segment(F, B)
FC = Segment(F, C)
e: FA / AH
BC' = Segment(B, C')
c: (x - x(F))^(2) + (y - y(F))^(2) = (e^(2) * x^(2))
F'A = Segment(F', A)
e': F'A / AH
f: (x - x(F'))^(2) + (y - y(F'))^(2) = (e'^(2) * x^(2))