xxxxxxxxxx
from itertools import product
def points_in_circle(radius):
for x, y in product(range(int(radius) + 1), repeat=2):
if x**2 + y**2 <= radius**2:
yield from set(((x, y), (x, -y), (-x, y), (-x, -y),))
list(points_in_circle(2))