How to numerically create a Delone set, $Del(r, R)$, in a square box in Euclidean space?

Question

I want to numerically create a Delone set, $Del\,(r, R)$, in a square box in Euclidean space. I was successful in creating a set of points which are at least $2*r$ distance away from each other. But the next challenge is to make sure that there's no empty sphere of radius $R$ in the box. I am unable to find an algorithm to successfully satisfy this condition.