3次元空間において、r = (x, y, z)を1点決めて、
これを中心とした半径nの球を考えます。
次に、この球の内部をグリッドで細かく立方体に分割します。
この立方体の辺の長さはmとします。
このとき、グリッドの交点の座標をすべて計算したいです。

やりたいことのイメージは、3次元の球で図を書くのが難しいので、
2次元の円で図を書きました。以下がそれです。
グリッドの交点の座標（赤色の点）をすべて求めたいです。

このようなことを実現するコードは、Pythonでどのように書けば良いでしょうか？
for文で座標をmずつずらすのをぜんぶ書けば良さそうだと思ったのですが、
球体だと効率的に書くのが思ったより難しい気がしました。

どなたかご教授いただけるとありがたいです。