Rotate around 
x-axis

Rotate around 
y-axis

vertices

lateral surfaces

light rays

projection

I_1 = (0, 6)
k_2 = Circle(I_1, 0.6)
K_1 = Point(k_2)
J_1 = Point(k_2)
b_1 = Angle(J_1, I_1, K_1)
I = (-2, 6)
k_1 = Circle(I, 0.6)
K = Point(k_1)
J = Point(k_1)
a_1 = Angle(J, I, K)
q = -0.5 (slider from -5 to 5, step 0.1)
a: a_1
b: b_1
C: (q * (cos(a) * (cos(b) - sin(b)) + sin(a)) + cos(b) + sin(b), q * (cos(a) * (cos(b) - sin(b)) + sin(a)) - sin(a) * (cos(b) - sin(b)) + cos(a))
A: (q * (cos(a) * (cos(b) + sin(b)) - sin(a)) - cos(b) + sin(b), q * (cos(a) * (cos(b) + sin(b)) - sin(a)) - sin(a) * (cos(b) + sin(b)) - cos(a))
B: (q * (cos(a) * (cos(b) - sin(b)) - sin(a)) + cos(b) + sin(b), q * (cos(a) * (cos(b) - sin(b)) - sin(a)) - sin(a) * (cos(b) - sin(b)) - cos(a))
D: (q * (cos(a) * (cos(b) + sin(b)) + sin(a)) - cos(b) + sin(b), q * (cos(a) * (cos(b) + sin(b)) + sin(a)) - sin(a) * (cos(b) + sin(b)) + cos(a))
E: (q * (cos(a) * (-cos(b) + sin(b)) - sin(a)) - cos(b) - sin(b), q * (cos(a) * (-cos(b) + sin(b)) - sin(a)) - sin(a) * (-cos(b) + sin(b)) - cos(a))
F: (q * (cos(a) * (-cos(b) - sin(b)) - sin(a)) + cos(b) - sin(b), q * (cos(a) * (-cos(b) - sin(b)) - sin(a)) - sin(a) * (-cos(b) - sin(b)) - cos(a))
G: (q * (cos(a) * (-cos(b) - sin(b)) + sin(a)) + cos(b) - sin(b), q * (cos(a) * (-cos(b) - sin(b)) + sin(a)) - sin(a) * (-cos(b) - sin(b)) + cos(a))
H: (q * (cos(a) * (-cos(b) + sin(b)) + sin(a)) - cos(b) - sin(b), q * (cos(a) * (-cos(b) + sin(b)) + sin(a)) - sin(a) * (-cos(b) + sin(b)) + cos(a))
A': (1 / 2 * q * ((-1 + cos(a)) * cos(b) - sin(a) + (1 + cos(a)) * sin(b) - 8) + 1 / 2 * ((-1 + cos(a)) * cos(b) - sin(a) + (1 + cos(a)) * sin(b) + 8), -cos(a) + sin(a) * (-sin(b) - cos(b)) + 1 / 2 * q * ((-1 + cos(a)) * cos(b) - sin(a) + (1 + cos(a)) * sin(b) - 8))
B': (1 / 2 * q * ((1 + cos(a)) * cos(b) - sin(a) + (1 - cos(a)) * sin(b) - 8) + 1 / 2 * ((1 + cos(a)) * cos(b) - sin(a) + (1 - cos(a)) * sin(b) + 8), -cos(a) + sin(a) * (sin(b) - cos(b)) + 1 / 2 * q * ((1 + cos(a)) * cos(b) - sin(a) + (1 - cos(a)) * sin(b) - 8))
C': (1 / 2 * q * ((1 + cos(a)) * cos(b) + sin(a) + (1 - cos(a)) * sin(b) - 8) + 1 / 2 * ((1 + cos(a)) * cos(b) + sin(a) + (1 - cos(a)) * sin(b) + 8), cos(a) + sin(a) * (sin(b) - cos(b)) + 1 / 2 * q * ((1 + cos(a)) * cos(b) + sin(a) + (1 - cos(a)) * sin(b) - 8))
D': (1 / 2 * q * ((-1 + cos(a)) * cos(b) + sin(a) + (1 + cos(a)) * sin(b) - 8) + 1 / 2 * ((-1 + cos(a)) * cos(b) + sin(a) + (1 + cos(a)) * sin(b) + 8), cos(a) + sin(a) * (-sin(b) - cos(b)) + 1 / 2 * q * ((-1 + cos(a)) * cos(b) + sin(a) + (1 + cos(a)) * sin(b) - 8))
i = Segment(A', B')
j = Segment(B', C')
k = Segment(C', D')
l = Segment(D', A')
E': (1 / 2 * q * ((-1 - cos(a)) * cos(b) - sin(a) + (-1 + cos(a)) * sin(b) - 8) + 1 / 2 * ((-1 - cos(a)) * cos(b) - sin(a) + (-1 + cos(a)) * sin(b) + 8), -cos(a) + sin(a) * (-sin(b) + cos(b)) + 1 / 2 * q * ((-1 - cos(a)) * cos(b) - sin(a) + (-1 + cos(a)) * sin(b) - 8))
F': (1 / 2 * q * ((1 - cos(a)) * cos(b) - sin(a) + (-1 - cos(a)) * sin(b) - 8) + 1 / 2 * ((1 - cos(a)) * cos(b) - sin(a) + (-1 - cos(a)) * sin(b) + 8), -cos(a) + sin(a) * (sin(b) + cos(b)) + 1 / 2 * q * ((1 - cos(a)) * cos(b) - sin(a) + (-1 - cos(a)) * sin(b) - 8))
m = Segment(B', F')
n = Segment(F', E')
p = Segment(A', E')
G': (1 / 2 * q * ((1 - cos(a)) * cos(b) + sin(a) + (-1 - cos(a)) * sin(b) - 8) + 1 / 2 * ((1 - cos(a)) * cos(b) + sin(a) + (-1 - cos(a)) * sin(b) + 8), cos(a) + sin(a) * (sin(b) + cos(b)) + 1 / 2 * q * ((1 - cos(a)) * cos(b) + sin(a) + (-1 - cos(a)) * sin(b) - 8))
H': (1 / 2 * q * ((-1 - cos(a)) * cos(b) + sin(a) + (-1 + cos(a)) * sin(b) - 8) + 1 / 2 * ((-1 - cos(a)) * cos(b) + sin(a) + (-1 + cos(a)) * sin(b) + 8), cos(a) + sin(a) * (-sin(b) + cos(b)) + 1 / 2 * q * ((-1 - cos(a)) * cos(b) + sin(a) + (-1 + cos(a)) * sin(b) - 8))
r = Segment(E', H')
s = Segment(G', H')
t = Segment(G', F')
i_1 = Segment(G', C')
j_1 = Segment(H', D')
l_1 = Ray(A', A)
m_1 = Ray(B', B)
n_1 = Ray(F', F)
p_1 = Ray(E', E)
q_1 = Ray(G', G)
r_1 = Ray(H', H)
s_1 = Ray(D', D)
t_1 = Ray(C', C)
L = (6.5615613973377025, -0.27645510002626955)
M = (5.5831648468670325, 4.395043922643557)
N = (7.195452120177855, 2.3969101223865517)
p_2 = Ellipse(L, M, N)
c = Segment(A, B)
d = Segment(A, E)
e = Segment(E, F)
f = Segment(B, F)
g = Segment(F, G)
h = Segment(H, E)
c_1 = Segment(D, A)
d_1 = Segment(C, B)
e_1 = Segment(D, C)
f_1 = Segment(H, G)
g_1 = Segment(G, C)
h_1 = Segment(D, H)
o = false
Vieleck1, a_2, b_2, c_2, d_2 = Polygon(A, B, C, D)
Vieleck2, a', b', c', d' = Polygon(A', B', C', D')
u = true
Vieleck3, c_3, b_3, f_2, g_2 = Polygon(C, B, F, G)
Vieleck4, c'_1, b'_1, f', g' = Polygon(C', B', F', G')
Vieleck5, c_4, d_3, h_2, g_3 = Polygon(C, D, H, G)
Vieleck6, c'_2, g'_1, h', d'_1 = Polygon(C', G', H', D')
Vieleck7, e_2, f_3, g_4, h_3 = Polygon(E, F, G, H)
Vieleck8, e', f'_1, g'_2, h'_1 = Polygon(E', F', G', H')
Vieleck9, a_3, b_4, f_4, e_3 = Polygon(A, B, F, E)
Vieleck10, a'_1, b'_2, f'_2, e'_1 = Polygon(A', B', F', E')
Vieleck11, a_4, d_4, h_4, e_4 = Polygon(A, D, H, E)
Vieleck12, a'_2, d'_2, h'_2, e'_2 = Polygon(A', D', H', E')
v = true
w = true