% ref Mandlebrot, FGoN, p.32 and p.48
\documentclass[border=5mm]{standalone}
\usepackage{luamplib}
\begin{document}
\mplibtextextlabel{enable}
\begin{mplibcode}
path gene; gene = (0,0) -- (1,0) -- (1, -1) -- (2, -1);
gene := gene rotated - angle point length gene of gene;
gene := gene scaled (1 / abs point length gene of gene);
beginfig(1);
for n=5, 6:
path p; p = gene scaled 300;
for i = 1 upto n:
p := for t = 1 upto length p:
subpath (0, 2) of gene
zscaled (point t of p - point t-1 of p)
shifted point t-1 of p -- endfor point length p of p;
endfor
p := p shifted (0, -144n);
for i = 0 upto length p:
fill fullcircle scaled 4
shifted point i of p
withcolor 7/8;
endfor
draw p withcolor 1/2;
endfor
endfig;
\end{mplibcode}
\end{document}