﻿rm(list = ls())  # ワークスペースをクリア
n <- 100000  # 試行回数
total <- numeric(n)  # 3個のさいころの目の合計を格納するベクトル
set.seed(1)  # 乱数のシード値を固定
for (i in 1:n) {
  total[i] <- sum(sample(1:6, 3, replace=TRUE))
}
# 合計が9になる確率を表示 (※モンテカルロ法で求めた値，理論値の順)
sum(total == 9) / n ; 25 / 216
# 合計が10になる確率を表示 (※モンテカルロ法で求めた値，理論値の順)
sum(total == 10) / n ; 1 / 8