A = Point(Circle[(1, 0, 0), Distance[(0, 0, 0), (1, 0, 0)], Segment[(0, 0, 0), (1, 0, 0)]])
a, B, G, H, E, F, face1, face3, face2, face4, faceABFE, faceCDEF, edge2, edge1, edge4, edgeAB, edge3, edgeBF, edgeCF, edgeCD, edge5, edgeDE, edgeAE, edgeEF = Cube((0, 0, 0), (1, 0, 0), A)
f = Line((0.5, 0.5, 0), Vector[(0, 0, 1)])
g = Line((0.5, 0, 0.5), Vector[(0, 1, 0)])
h = Line((0, 0.5, 0.5), Vector[(1, 0, 0)])
G_1 = Point(h)
b = Sphere(G_1, H)
H_1 = Point(h)
c = Sphere(H_1, G)
I = Point(g)
d = Sphere(I, G)
J = Point(g)
e = Sphere(J, E)
K = Point(f)
i = Sphere(K, H)
L = Point(f)
j = Sphere(L, A)
k = IntersectConic(d, i)
C: (0, 0)
D: (1, 0)
M, H' = Intersect(k, b)
G', O_1 = Intersect(k, c)
p = IntersectConic(d, b)
N, D' = Intersect(p, j)
q = IntersectConic(j, c)
B', Q = Intersect(q, e)
R, C' = Intersect(q, d)
r = IntersectConic(j, e)
A', S = Intersect(r, b)
s = IntersectConic(i, e)
F', T = Intersect(s, c)
P, E' = Intersect(s, b)
t = Circle(G', C', D')
c_1 = Circle(D', A', E')
d_1 = Circle(F', E', A')
e_1 = Circle(H', E', F')
f_1 = Circle(B', C', G')
g_1 = Circle(A', D', C')
U = Center(e_1)
p_1 = Plane(H', F', E')
l = OrthogonalLine(U, p_1)
q_1 = Plane(C', D', A')
V = Center(g_1)
m = OrthogonalLine(V, q_1)
O = Intersect(m, l)
n = Sphere(O, C')
o = IsInRegion(A', n)
u = IsInRegion(B', n)
v = IsInRegion(C', n)
w = IsInRegion(D', n)
a_1 = IsInRegion(E', n)
b_1 = IsInRegion(F', n)
h_1 = IsInRegion(G', n)
i_1 = IsInRegion(H', n)