---

r[0] := sphere([3, 0, 0], 1, grid = [25, 25]):
a := (1/4) * Pi:
c := 1:
while evalf(a - 2*Pi) < 0 do r[c] := rotate(r[0], a, [[0, 0, 0], [0, 0, 1]]); a := a + (1/4) * Pi; c := c + 1 end do:
r := display([seq(r[i], i = 0..c-1)], scaling = constrained, style = hidden, lightmodel = light4, orientation = [10, -125], shading = zhue):

for i from 1 to 10 do q[i] := homothety(r, 0.5^i): w[i] := translate(q[i], 0, 0, i/2): end do:
qq := display([seq(q[i], i = 1..10)]):
qa := display(qq, r);
display(qa);