{
	p = [0, 1, 2, 2+I];

	for (n=1, 5,
		e = p[#p];
		p = concat(p, apply(z -> (z-e)*I+e+I, Vecrev(p)));
		e = p[#p];
		p = concat(p, apply(z -> (z-e)/I+e+1, Vecrev(p)));
	);

	n = -1;
	apply (z -> print (n++ " " real(z)), p);
}

quit