Quasicrystals

wave2 = Compile[{{x, _Real}, {y, _Real}, {t, _Real}},
   Sum[Cos[Sin[\[Theta]] x + Cos[\[Theta]] y + t],
    {\[Theta], 0, Pi (1 - 1/n), Pi/n}]
   , CompilationTarget -> "C"];
n = 7;
width = 200;
X2[t_] := Table[wave2[x, y, t]
   , {x, -width, width}, {y, -width, width}];
f[x_] := (1 + Tanh[10 (x - 0.5)])/2;
Manipulate[
 ArrayPlot[X2[t],
  ColorFunction -> Function[x, GrayLevel[f[x]]]
  ]
 , {t, 0, 2 Pi}
 ]