A = (1.3402344343237687, 1.2786458656979036, 0)
B = PointIn(Sphere[A, 1])
a = Segment(A, B)
b, C, D, E, F, G, H, faceABCD, faceADHE, faceABFE, faceBCGF, faceCDHG, faceEFGH, edgeAD, edgeAB, edgeBC, edgeCD, edgeAE, edgeDH, edgeEH, edgeEF, edgeBF, edgeFG, edgeCG, edgeGH = Cube(A, B, xOyPlane)
D' = Mirror(D, H)
D'_1 = Mirror(D, C)
I = Point(edgeAD)
c = Segment(D, D')
d = Segment(D'_1, I)
e = OrthogonalLine(D'_1, xOyPlane)
f = Plane(H, E, A)
g = OrthogonalLine(D', f)
J = Intersect(g, e)
A' = Mirror(A, I)
poly1, d', j, d'_1, i = Polygon(D', J, D'_1, I)
k, O, P, Q, R, faceD'JPO, face4, face3, faceD'IRO, faceOPQR, edgeD'O, edgeJP, edge2, edgeIR, edgeOP, edgePQ, edgeQR, edgeOR = Prism(poly1, 1)
K = Midpoint(O, R)
L = Midpoint(Q, R)
K' = Mirror(K, O)
L' = Mirror(L, Q)
h = Segment(K', A')
l = Segment(R, L')
m = OrthogonalLine(L', xOyPlane)
n = OrthogonalLine(K', f)
M = Intersect(n, m)
poly2, m_1, k', r, l' = Polygon(M, K', R, L')
p, V, W, Z, A_1, faceK'MVW, faceK'RZW, face1, face2, face5, edgeMV, edgeK'W, edgeRZ, edge1, edgeVW, edgeWZ, edge4, edge3 = Prism(poly2, -1)
q = Line(D'_1, Q)
s = Line(D', O)
t = Line(A_1, Z)
a_1 = Line(Z, W)
N = Intersect(q, t)
S = Intersect(s, a_1)
N' = Mirror(N, A_1)
S' = Mirror(S, W)
b_1 = OrthogonalLine(N', xOyPlane)
c_1 = OrthogonalLine(S', f)
T = Intersect(c_1, b_1)
poly3, s', t_1, n', z_1 = Polygon(S', T, N', Z)
u, E_1, F_1, G_1, H_1, face14, face13, face11, face12, face15, edge14, edge15, edge13, edge12, edge16, edge18, edge19, edge17 = Prism(poly3, 1)
U = Mirror(A_1, N')