Ω: (0, 0, 0)
z0: z(Ω)
y0: y(Ω)
x0: x(Ω)
c = 10 (slider from 0 to 10)
yA: c / 2
xA: yA
h = 7.0710678118654755 (slider from 0 to ωa)
A: (xA + x0, yA + y0, z0)
B: ((-xA) + x0, yA + y0, z0)
lc = 20 (slider from 0 to 50)
poly1, f, g, i, j, C, D = Polygon(A, B, 4, xOyPlane)
E: (x0, y0, z0 + h)
F: (x0, y0, z0 - h)
t1, e, a, b = Polygon(A, B, E)
t2, e_1, b_1, c_1 = Polygon(B, C, E)
t3, e_2, c_2, d = Polygon(C, D, E)
t4, e_3, d_1, a_1 = Polygon(D, A, E)
t5, f_1, a_2, b_2 = Polygon(A, B, F)
t6, f_2, b_3, c_3 = Polygon(B, C, F)
t7, f_3, c_4, d_2 = Polygon(C, D, F)
t8, f_4, d_3, a_3 = Polygon(D, A, F)
r = 4.083 (slider from 0 to (2 * h), step 0.001)
k = Sphere(Ω, r)
I: (A + B + E) / 3
J: (B + C + E) / 3
K: (C + D + E) / 3
L: (D + A + E) / 3
I': (A + B + F) / 3
J': (B + C + F) / 3
K': (C + D + F) / 3
L': (D + A + F) / 3
l = Line(I, K')
m = Line(J, L')
n = Line(K, I')
p = Line(L, J')
ωi: nroot((x(I) - x0)^(2) + (y(I) - y0)^(2) + (z(I) - z0)^(2),2)
ωa: nroot((x(A) - x0)^(2) + (y(A) - y0)^(2) + (z(A) - z0)^(2),2)
r': nroot(r^(2) - ωi^(2),2)
q = Circle(I, r', l)
s = Circle(J, r', m)
t = Circle(K, r', n)
g_1 = Circle(L, r', p)
h_1 = Circle(I', r', n)
k_1 = Circle(J', r', p)
p_1 = Circle(K', r', l)
q_1 = Circle(L', r', m)
o, r_1, u = Cylinder(q, (-lc))
v, v_2, v_3 = Cylinder(s, (-lc))
w, w_2, w_3 = Cylinder(t, (-lc))
i_1, s_1, j_1 = Cylinder(g_1, (-lc))
l_1, t_1, m_1 = Cylinder(h_1, lc)
n_1, n_2, n_3 = Cylinder(k_1, lc)
o_1, o_2, o_3 = Cylinder(p_1, lc)
u_1, g_2, v_1 = Cylinder(q_1, lc)
to = 1 (slider from 0 to 1)
tc = 1 (slider from 0 to 1)
tb = 1 (slider from 0 to 1, step 1)