Script

g(x) = x
k = 2.13 (slider from 1 to 10, step 0.01)
f(x) = x^(k) / (x^(k) + (1 - x)^(k))
d = 10 (slider from 0 to 10, step 0.01)
T = 0.59 (slider from 0 to 1, step 0.01)
K = 10 (slider from 1 to 10, step 1)
L = Join({{0}, Sequence[k, k, 1, (6 * K) - 1]})
A = (-4.32, -5)
B: (4.32, -5)
H(x, y) = 1 / 3^((-g(x)) + floor(y / 6))
size = Length(L)
a: 3^(floor(1 / 6))
Fx(x, y, z) = (x * sin((π / 3 * (y + ((-1)^(floor(y / 6)) * f(g(z)))))) / 3^((-g(z)) + floor(y / 6)))
Fy(x, y, z) = (x * cos((π / 3 * (y + ((-1)^(floor(y / 6)) * f(g(z)))))) / 3^((-g(z)) + floor(y / 6)))
Vec = Sequence(Vector[(0, 0), (Fx(d, Element[L, k], T), Fy(d, Element[L, k], T))], k, 1, Length[L])
C = Sequence((x(Element[Vec, k]), y(Element[Vec, k])), k, 1, Length[L])
G(x, y) = ((f(g(x)) * π / 3) * (-1)^(floor(y / 6)))
K_0 = Sequence(Translate[pic1, Vector[Element[Vec, k]]], k, 1, Length[L])
K_1 = Sequence(Dilate[Element[K_0, k], H(T, Element[L, k]), Element[C, k]], k, 1, Length[L])
K_2 = Sequence(Rotate[Element[K_1, k], G(T, Element[L, k]), Element[C, k]], k, 1, Length[L])