O = (2, 1)
V = Point(yAxis)
c = Circle(O, V)
B = Point(c)
C = Point(c)
A = Point(c)
poly1, c_1, a, b = Polygon(A, B, C)
d = OrthogonalLine(A, xAxis)
A1 = Point(xAxis)
A2 = Intersect(c, d, 2.0)
e = OrthogonalLine(B, b)
B1 = Point(b)
B2 = Intersect(c, e, 2.0)
f = OrthogonalLine(C, c_1)
C1 = Point(c_1)
C2 = Intersect(c, f, 2.0)
H = Point(d)
C'' = Midpoint(C, H)
B'' = Midpoint(B, H)
A'' = Midpoint(A, H)
A' = Midpoint(B, C)
B' = Midpoint(A, C)
C' = Midpoint(A, B)
Ω = Midpoint(O, H)
g = Circle(Ω, B')
h = Line(A1, B'')
i = Line(B1, A'')
R1 = Intersect(h, i)
j = Line(C, R1)
k = Line(C'', B1)
l = Line(B'', C1)
P1 = Intersect(k, l)
m = Line(P1, R1)
n = Line(P1, Ω)
p = OrthogonalLine(Ω, m)
S = Intersect(m, p)
q = Line(B2, A)
r = Line(A2, B)
R2 = Intersect(q, r)
s = Line(B, C2)
t = Line(C, B2)
P2 = Intersect(s, t)
a_1 = Line(P2, R2)
I' = Intersect(l, t)
K = Intersect(h, q)
b_1 = Line(K, I')