Mathematica code: Tile[k_, rx_, ry_, x_, y_, r_] := Table[ Translate[ Rotate[ {AbsoluteThickness[k], Circle[{i, i}, {rx, ry}, {i*Pi, Pi/2 + i*Pi}]}, r, {.5, .5}], {x, y}], {i, 0, 1, 1}]rr[Q_] := (SeedRandom[Q]; RandomReal[])GCOS[k_, rx1_,