﻿import Random  # 乱数を扱うためのモジュールをインポート
n = 100000  # 試行回数
Random.seed!(1)  # 乱数のシード値を固定
X = rand(n)  # U(0,1)に従う乱数をn個生成し，Xに代入
Y = (2 .- 2 .* X) .* rand(n)  # U(0,2-2X)に従う乱数をn個生成し，Yに代入
# 確率を表示
println(count(Y .>= 1) / n)  # モンテカルロ法
println((1/2) * (1 - log(2)))  # 参考：理論値