Ещё, на всякий случай, озвучу такую мысль:
Можно попробовать оптимизировать. Перед тем как вычислять сумму квадратов разностей, то есть попадание в окружность, можно сразу отсеять то, что не попадает в описанный квадрат. Как бы:
DX:=abs(x1-x2);
IF DX>R THEN EXIT;
DY:=abs(y1-y2);
IF DY>R THEN EXIT;
IF SQR(DX)+SQR(DY)>R2 THEN EXIT;