O = Intersect(xAxis, yAxis)
Γ = Circle(O, 1)
A' = Point(Γ)
B' = Point(Γ)
C' = Point(Γ)
ta = Tangent(A', Γ)
tc = Tangent(C', Γ)
tb = Tangent(B', Γ)
C_t = Intersect(ta, tb)
ck = Line(C_t, C')
A_t = Intersect(tc, tb)
ak = Line(A_t, A')
K = Intersect(ck, ak)
ab = Line(A', B')
ac = Line(A', C')
bc = Line(B', C')
B_t = Intersect(ta, tc)
AB = Segment(A', B')
BC = Segment(B', C')
AC = Segment(C', A')
ok = Line(K, O)
S, S' = Intersect(Γ, ok)
S_{A} = Midpoint(S, Mirror[S, bc])
S_{B} = Midpoint(S, Mirror[S, ac])
S_{C} = Midpoint(S, Mirror[S, ab])
ss = Line(S_{B}, S_{C})
G': (A' + B' + C') / 3
gy = Line(G', ss)
f = Segment(S, S_{C})
g = Segment(S, S_{A})
h = Segment(S, S_{B})
i = OrthogonalLine(G', gy)
A = Midpoint(B', C')
B = Midpoint(A', C')
C = Midpoint(A', B')
D = Midpoint(A', G')
E = Midpoint(C', G')
c = Conic(C, D, B, E, A)
G_2, G_1 = Focus(c)