Bib-bridgslam2

c: 10
x0 = 2.13 (slider from 0 to c / 2, step 0.01)
A: (c / 2, ((-c)) / 2)
B: (c / 2, c / 2)
poly1, f, g, h, i, C, D = Polygon(A, B, 4)
α = 1.0471975511965976 rad (slider from 45° to αmax, step 0.1°)
E: (c / 2 - x0, ((-c)) / 2)
A' = Rotate(A, α, E)
j = Ray(E, A')
k = Ray(A, C)
F = Intersect(j, k)
G: (c / 2, ((-c)) / 2 + x0)
q1, a_0, e, f_1, g_1 = Polygon(A, E, F, G)
A'_{1} = Mirror(A, xAxis)
E1 = Mirror(E, xAxis)
F1 = Mirror(F, xAxis)
G1 = Mirror(G, xAxis)
q1', a'_{1}, e', f', g' = Polygon(A'_{1}, E1, F1, G1)
A'_{2} = Mirror(A, yAxis)
E3 = Mirror(E, yAxis)
F3 = Mirror(F, yAxis)
G3 = Mirror(G, yAxis)
q1'_{1}, a'_{2}, e'_{1}, f'_{1}, g'_{1} = Polygon(A'_{2}, E3, F3, G3)
H = Mirror(A'_{1}, yAxis)
E2 = Mirror(E1, yAxis)
F2 = Mirror(F1, yAxis)
G2 = Mirror(G1, yAxis)
q1'', h_1, e'', f'', g'' = Polygon(H, E2, F2, G2)
q2, f_2, l, m, n = Polygon(F, F1, F2, F3)
q3, f_3, g_2, p, q = Polygon(F, G, G1, F1)
q4, r, s, t, b = Polygon(F1, E1, E2, F2)
q5, d, i_1, j_1, k_1 = Polygon(F2, G2, G3, F3)
q6, l_1, m_1, e_1, f_4 = Polygon(F3, E3, E, F)
n_1 = Segment(G2, G1)
p_1 = Segment(G, G3)
q_1 = Segment(E, E1)
r_1 = Segment(E2, E3)
I = Intersect(p_1, q_1)
K = Intersect(n_1, r_1)
L = Intersect(p_1, r_1)
J = Intersect(n_1, q_1)
q7, i_2, j_2, k_2, l_2 = Polygon(I, J, K, L)
hpyr: (x0 * sqrt((-1 + tan(α)) / (1 + tan(α))))
I0: (c / 2 - x0, ((-c)) / 2 + x0, hpyr)
J0: (c / 2 - x0, c / 2 - x0, hpyr)
L0: (((-c)) / 2 + x0, ((-c)) / 2 + x0, hpyr)
K0: (((-c)) / 2 + x0, c / 2 - x0, hpyr)
q8, i0, j0, k0, l0 = Polygon(I0, J0, K0, L0)
s_1 = Segment(I0, I)
t_1 = Segment(I, E)
M = Intersect(f_4, t_1)
a_1 = Segment(I0, M)
b_1 = Segment(E, M)
c_1 = Segment(M, I)
a: (x0 * tan(α) / (1 + tan(α)))
q9, i_3, j_3, k_3, l_3 = Polygon(I, J, K, L)
β = 2.1862760354652844 rad (slider from 0 to βmax, step 0.1°)
F' = Rotate(F, β, q)
G' = Rotate(G, β, q)
G1' = Rotate(G1, β, q)
F1' = Rotate(F1, β, q)
q3', d_1, h_2, g1', f1' = Polygon(F', G', G1', F1')
F1'_{1} = Rotate(F1, β, b)
E1' = Rotate(E1, β, b)
E2' = Rotate(E2, β, b)
F2' = Rotate(F2, β, b)
q4', f1'_{1}, e1', e2', f2' = Polygon(F1'_{1}, E1', E2', F2')
F3' = Rotate(F3, (-β), f_4)
E3' = Rotate(E3, (-β), f_4)
E' = Rotate(E, (-β), f_4)
F'_{1} = Rotate(F, (-β), f_4)
F2'_{1} = Rotate(F2, (-β), k_1)
G2' = Rotate(G2, (-β), k_1)
G3' = Rotate(G3, (-β), k_1)
F3'_{1} = Rotate(F3, (-β), k_1)
F3'_{2} = Rotate(F3, β, f_4)
E3'_{1} = Rotate(E3, β, f_4)
E'_{1} = Rotate(E, β, f_4)
F'_{2} = Rotate(F, β, f_4)
q6', f3'_{2}, e3'_{1}, m_2, f'_{2} = Polygon(F3'_{2}, E3'_{1}, E'_{1}, F'_{2})
F2'_{2} = Rotate(F2, β, k_1)
G2'_{1} = Rotate(G2, β, k_1)
G3'_{1} = Rotate(G3, β, k_1)
F3'_{3} = Rotate(F3, β, k_1)
q5', f2'_{2}, g2'_{1}, g3'_{1}, f3'_{3} = Polygon(F2'_{2}, G2'_{1}, G3'_{1}, F3'_{3})
γ = Angle(A, E, F)
αmax: 180° - (atan(c / (c - (2 * x0))) * (180°) / π)
βmax: acosd((a - x0) / a)
Vol: ((q2 + q8) * hpyr / 2)
N = (-5, 9)
O = Point(Circle[N, Vol / 10])
Volume = Segment(N, O)