# モンテカルロ法のヒストグラム

# 試行回数
trial <- 100

# 点の個数
n <- 1000

# 面積をn回求める
r <- replicate(n, sum(runif(n, -1, 1) ^ 2 + runif(n, -1, 1) ^ 2 < 1) / n * 4)
print(r)

# 面積のヒストグラム（相対度数）
hist(r, freq = FALSE)