20/03/17 21:18:37 k85T9ON2.net
こっちの方がx,y,zともに一様分布になっている。
vertex <- function(r=1){
x=runif(1,-1,1) # x ~ 一様分布[-1,1]
phi=runif(1,-pi,pi) # φ ~ 一様分布[-π,π]
y=sqrt(1-x^2)*cos(phi) # √(1-x^2)*cos(φ)
z=sqrt(1-x^2)*sin(phi) # √(1-x^2)*sin(φ)
r*c(x,y,z)
}
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
これでやってみると、四面体の場合
> mean(tetra)
[1] 0.1201118
> summary(tetra)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0000003 0.0372266 0.0922805 0.1201118 0.1794858 0.5104545