So far I came up with this.
PointsList =
Flatten[Table[{i, j}, {i, -4, 4, 0.5}, {j, -4, 4, 0.5}],
1];(*generate poinst list*)
For[k = 1,
x^2 + y^2 <= 1 /. {x -> PointList[[k, 1]], y -> PointList[[k, 2]]},
k = k + 1]
I still need to go back and check how to take things off the list