{
	squ = [];

	mx = 200;
	mmx = 2000;

	for (x=0, oo,
		if (norm(x^2) > mmx^2,
			break;
		);
		for (y=0, oo,
			z = (x+I*y)^2;
			if (norm(z) > mmx^2,
				break;
			);
			squ = concat(squ, [ (+x+I*y)^2, (+x-I*y)^2, (-x+I*y)^2, (-x-I*y)^2 ]); 
		);
	);

	squ = Set(squ);

	i = 0;
	for (d=0, 200,
		for (y=0, d,
			x = d-y;
			my (d=vecmin(apply(z -> norm( (x+I*y) - z ), squ)));
			print (i++ " " d);
		);
	);
}

quit