A = (-1.7801748258528132, 6.768877063792959, -2.569961144618574)
B = (-4.557293992161333, -1.355134905003568, -6.938893903907228E-18)
C = (2.68052451778626, 0.2242050165413687, 0)
poly1, c, a, b = Polygon(A, B, C)
D = (-0.6291514229935149, 0.9264420369595607, 4.065159090857984)
d, faceABD, faceBCD, faceACD, edgeAD, edgeBD, edgeCD = Pyramid(poly1, D)
E = Midpoint(edgeCD)
F = Midpoint(edgeAD)
G = Midpoint(c)
H = Midpoint(a)
I = Midpoint(edgeBD)
J = Midpoint(b)
e = OrthogonalPlane(E, edgeCD)
f = OrthogonalPlane(J, b)
g = OrthogonalPlane(F, edgeAD)
h = IntersectPath(e, f)
i = OrthogonalPlane(G, c)
j = IntersectPath(i, e)
k = OrthogonalPlane(J, b)
l = OrthogonalPlane(G, c)
m = IntersectPath(l, k)
K = Intersect(h, j)
n = Sphere(K, D)
p = Circle(D, C, A)
q = Circle(B, A, D)
r = Circle(B, C, D)
s = Circle(C, B, A)