> 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の記事