>>432
θが[0,Δθ] になる領域は広いが、[π/2-Δθ,π/2] の領域は狭い。
これに同じ個数の点を打てば、θ＝π/2 付近では密集する。
領域の中の点の数と面積が比例すればいい。

-1≦z≦1,0≦φ＜2πの一様乱数を取って
x=√(1-z^2) cosφ, y=x=√(1-z^2) sinφ
とするといいらしい。