θφで面素密度に合わせて乱数にしたら良いと思う
dS=cosθdθdφなのでθという値を取る確率をcosθにする
つまりzθの長方形で乱数発生させてz>cosθは除外してθを取る

dV=dxdydz=rdrdSだから
xyzで球の外を除外して正規化しても本質は同じだな