\documentclass[border=5mm]{standalone}

\documentclass[border=5mm]{standalone}
\usepackage{luamplib}
\begin{document}
\begin{mplibcode}
beginfig(1);
path s; s = unitsquare shifted -(1/2, 1/2) scaled 21;
for n = 1 upto 4:
numeric N; N = 9n + 8;
numeric a; a = 30 normaldeviate;
for t=0 upto N-1:
draw s rotated if odd n: - fi 13 shifted (5N * right) rotated (360 / N * t + a)
withpen pencircle scaled 1 withcolor t mod 2;
endfor
endfor
picture P;
P = currentpicture; clearit;
fill bbox P withcolor 0.52; draw P;
endfig;
\end{mplibcode}
\end{document}