Mathematica code:
Animate[
DensityPlot[
Sum[Cos[(Cos[n*2*Pi/23] + Sin[n*2*Pi/23])*x + (Cos[n*2*Pi/23] - Sin[n*2*Pi/23])*y + t],
{n, 0, 22, 1}], {x, -150, 150}, {y, -150, 150},
PlotPoints -> 150, Mesh -> False, Frame -> False, ColorFunction -> Hue, ImageSize -> 522],
{t, 0, 2*Pi, Pi/10}]