{ 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