# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add("Distributions")
import Random  # 乱数を扱うためのモジュールをインポート
using Distributions
function prob(k)
    Random.seed!(1)  # 乱数のシード値を固定
    # 二項分布B(10000, 0.5)に従う乱数をk個生成し，Snに代入
    Sn = rand(Binomial(10000, 0.5), k)
    target_prob = count(4900 .<= Sn .<= 5100) / k
    return target_prob
end
# モンテカルロ法で求めた確率を表示 (参考：理論値は0.9544)
println(prob(100))
println(prob(1000))
println(prob(10000))
println(prob(100000))