A = Point(xAxis)
B = Point(yAxis)
C = (0.9099446884904197, 2.2539144352419527, 0)
D = Point(zAxis)
a, faceABC, faceABD, faceBCD, faceACD, edgeAB, edgeBC, edgeAC, edgeAD, edgeBD, edgeCD = Pyramid(A, B, C, D)
f = AngularBisector(B, D, C)
g = OrthogonalLine(C, f, space)
i = AngularBisector(C, D, A)
j = OrthogonalLine(C, i, space)
E = Intersect(f, g)
F = Intersect(i, j)
k = Line(F, E)
G = Intersect(g, edgeBD)
H = Intersect(j, edgeAD)
I = Intersect(g, edgeBD)
l = Line(H, G)
poly1, g_1, c, h_1 = Polygon(C, H, G)
h = Segment(H, I)
m = Segment(C, H)
n = Segment(H, C)
p = Segment(I, C)