
Mathematica code:
Animate[ DensityPlot[
Sum[.20*Cos[(Cos[n*2*Pi/5] + Sin[n*2*Pi/5])*x + (Cos[n*2*Pi/5] - Sin[n*2*Pi/5])*y + t]^10,
{n, 0, 4, 1}], {x, -25, 25}, {y, -25, 25},
PlotPoints -> 50, Mesh -> False, Frame -> False, ColorFunction -> GrayLevel, ImageSize -> 521],
{t, 0, 2*Pi, Pi/20}]