面白い問題おしえて～な 31問目

面白い問題おしえて～ ..

855:１３２人目の素数さん
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

次ページ

続きを表示

1を表示