円板上一様分布

> rand1<-runif(1000)
> rand2<-runif(1000)
> rand2<-rand2*2*pi
> X<-rand1*sin(rand2)
> Y<-rand1*cos(rand2)
> X<-sqrt(rand1)*sin(rand2)
> Y<-sqrt(rand1)*cos(rand2)
> plot(X,Y)

参考にした乱数発生に関するMathematicaの記事