20/03/17 07:53:05.59 Ze9EuNOD.net
>>820
四面体で１００万回シミュレーションして平均値をだしてみた。
vertices <- function(r=1){
a=runif(2,-pi,pi) # 角度Θ,Φを一様分布で選らんで
theta=a[1]
phi=a[2]
x=r*sin(theta)*cos(phi)
y=r*sin(theta)*sin(phi)
z=r*cos(theta)
c(x,y,z) # 直交座標を返す
}
sim <- function(r=1){
vectors=replicate(4,vertices(r)) # ４点の直交座標
abs(det(vectors[,2:4]-vectors[,1]))/6 # 四面体の体積
}
k=1e6
re=replicate(k,sim())
mean(re)
> mean(re)
[1] 0.1069067