---

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): 
display(r, r);