{
	zz = [0, 1+I];
	e = 3;

	for (k=1, 6,
		zz = concat([
			zz,
			apply(z -> e + real(z)*(I-1) + imag(z)*(I), zz),
			apply(z -> e*I + z, zz),
			apply(z -> e*(1+I) + real(z)*(1-I) + imag(z)*(-I), zz)
		]);
		e *= 2;
	);

	zz = concat(zz, e);
	for (n=1, #zz,
		print (n-1 " " imag(zz[n]));
	);
}

quit