f: y = 0
g: y = 1
A = (-1.8182510397961338, 0.7867007412246129)
B = Point(g)
h = OrthogonalLine(B, f)
C = Intersect(f, h)
D = Midpoint(B, C)
d = Segment(B, D)
j = Segment(D, C)
k = Line(D, g)
c = Circle(A, Segment[D, B])
M, M' = Intersect(c, k)
d_{1} = Circle(M, A)
e = Circle(M', A)