牟合方蓋
Author: 藍邦偉　老師

單位向量

座標軸

x

y

z

牟合方蓋

O = (10.799224923904992, 3.561397916649132)
β = 1.3962634015954636 rad (slider from 0 to 6.283185307179586, step 0.1*°)
Z_{unit}: O + (0, sin(β))
w = Vector(O, Z_{unit})
α = 3.633775502652194 rad (slider from 185*° to 265*°, step 0.01*°)
X_{unit}: O + (cos(α), sin(α) * cos(β))
u = Vector(O, X_{unit})
Y_{unit}: O + (-1 * sin(α), cos(α) * cos(β))
v = Vector(O, Y_{unit})
L = 7 (slider from 1 to 10, step 1)
a = Vector(O - L * v, O + L * v)
b = Vector(O - L * w, O + L * w)
z = Vector(O - L * u, O + L * u)
c = false
d = false
Q_2: 0.2 * Corner[2] + 0.8 * Corner[1]
Q_3: (x(Q_2), y(Corner[3]))
g = Segment(Q_2, Q_3)
poly1, h_1, q_1, q_4, i = Polygon(Corner[1], Q_2, Q_3, Corner[4])
r = 1.7 (slider from 1 to 5, step 0.1)
h = 4.96 (slider from 1.8 * r to 5, step 0.1)
A_3: O + h * u + r * w
A_1: O + h * u
A_2: O + h * u + r * v
c_1 = CurveCartesian(x(O + r * cos(t) * v + r * sin(t) * w + h * u), y(O + r * cos(t) * v + r * sin(t) * w + h * u), t, 0, 1.5707963267948966)
c_2 = CurveCartesian(x(O + r * cos(t) * v + r * sin(t) * w), y(O + r * cos(t) * v + r * sin(t) * w), t, 0, 1.5707963267948966)
A_4: O + r * v
A_5: O + r * w
A_6: O
Poly1, a_1, a_2, a_4, a_6 = Polygon(A_1, A_2, A_4, A_6)
Poly2, a_3, a_5, a_7, a_8 = Polygon(A_1, A_3, A_5, A_6)
list_1 = Sequence(CurveCartesian[x(O + r * cos(t) * v + r * sin(t) * w + k * u), y(O + r * cos(t) * v + r * sin(t) * w + k * u), t, 0, 1.5708], k, 0, h, 0.25)
Tri1, a_9, a_{10}, a_{11} = Polygon(A_1, A_2, A_3)
B_1: O + h * v
B_2: O + h * v + r * u
B_3: O + h * v + r * w
c_3 = CurveCartesian(x(O + r * cos(t) * u + r * sin(t) * w + h * v), y(O + r * cos(t) * u + r * sin(t) * w + h * v), t, 0, 1.5707963267948966)
c_4 = CurveCartesian(x(O + r * cos(t) * u + r * sin(t) * w), y(O + r * cos(t) * u + r * sin(t) * w), t, 0, 1.5707963267948966)
B_4: O + r * u
list_2 = Sequence(CurveCartesian[x(O + r * cos(t) * u + r * sin(t) * w + k * v), y(O + r * cos(t) * u + r * sin(t) * w + k * v), t, 0, 1.5707963267948966], k, 0, h, 0.25)
Tri2, b_3, b_1, b_2 = Polygon(B_1, B_2, B_3)
Poly3, b_4, b_5, a_{12}, a_{13} = Polygon(B_1, B_3, A_5, A_6)
Poly4, b_6, b_7, a_{14}, b_8 = Polygon(B_2, B_1, A_6, B_4)
c_5 = CurveCartesian(x(O + r * cos(t) * u + r * sin(t) * w + r * cos(t) * v), y(O + r * cos(t) * u + r * cos(t) * v + r * sin(t) * w), t, 0, 1.5707963267948966)
θ = 0.5235987755982988 rad (slider from 0 to 90*°, step 0.01*°)
H_2: O + r * cos(θ) * u + r * sin(θ) * w
H_4: O + r * cos(θ) * v + r * sin(θ) * w
H_3: O + r * cos(θ) * u + r * sin(θ) * w + r * cos(θ) * v
H_1: O + r * sin(θ) * w
poly2, h_2, h_3, h_4, h_5 = Polygon(H_1, H_2, H_3, H_4)