rotation

O = Intersect(xAxis, yAxis)
B = Point(xAxis)
i = Vector(O, B)
C = Point(yAxis)
j = Vector(O, C)
theta = 2.705260340591211 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
phi = 4.014257279586958 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
u: cos(phi) * i - sin(theta) * sin(phi) * j
w: -1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * sin(phi) * i - sin(theta) * cos(phi) * j
v': -1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * sin(theta) * j
v: cos(theta) * j
A = placer(theta, phi, 5.0, 5.0, 5.0)
D = placer(theta, phi, -5.0, 5.0, 5.0)
E = placer(theta, phi, 5.0, -5.0, 5.0)
F = placer(theta, phi, 5.0, 5.0, -5.0)
G = placer(theta, phi, -5.0, -5.0, 5.0)
H = placer(theta, phi, -5.0, 5.0, -5.0)
I = placer(theta, phi, -5.0, -5.0, -5.0)
J = placer(theta, phi, 5.0, -5.0, -5.0)
a = Segment(D, G)
b = Segment(G, E)
c = Segment(E, A)
d = Segment(A, D)
e = Segment(H, D)
f = Segment(G, I)
g = Segment(E, J)
h = Segment(A, F)
k = Segment(F, J)
l = Segment(J, I)
m = Segment(I, H)
n = Segment(H, F)
K = Midpoint(E, G)
L = Midpoint(E, J)
M = Midpoint(J, I)
N = Midpoint(G, D)
P = Midpoint(D, A)
Q = Midpoint(A, E)
R = Midpoint(F, A)
S = Midpoint(H, D)
T = Midpoint(I, H)
U = Midpoint(J, F)
V = Midpoint(H, F)
W = Midpoint(G, I)
poly1, k_1, l_1, m_1, w_1 = Polygon(K, L, M, W)
poly2, l_2, k_2, q = Polygon(K, Q, L)
poly3, w_2, k_3, n_1 = Polygon(K, N, W)
poly4, m_2, l_3, r = Polygon(L, R, M)
poly5, t, w_3, m_3 = Polygon(W, M, T)
poly6, u_1, l_4, m_4 = Polygon(L, M, U)
poly7, m_5, w_4, t_1 = Polygon(W, T, M)
poly8, k_4, q_1, p, n_2 = Polygon(K, Q, P, N)
poly9, l_5, q_2, r_1, u_2 = Polygon(L, Q, R, U)
poly10, n_3, w_5, t_2, s = Polygon(N, W, T, S)
poly11, m_6, u_3, v_1, t_3 = Polygon(M, U, V, T)
poly12, r_2, p_1, q_3 = Polygon(P, Q, R)
poly13, s_1, p_2, n_4 = Polygon(P, N, S)
poly14, v_2, s_2, t_4 = Polygon(S, T, V)
poly15, u_4, v_3, r_3 = Polygon(V, R, U)
poly16, v_4, s_3, p_3, r_4 = Polygon(V, S, P, R)

Tool: point3D (command placer)
placer[theta,phi,x,y,z]
Input: theta, phi, a, b, c
Output: D

D: (-1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * a * sin(phi) + b * cos(phi), -1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * 1.0 * a * sin(theta) * cos(phi) - b * sin(theta) * sin(phi) + c * cos(theta))