# データをリセットして点の数を設定
rm(list = ls())
N <- 100000

# 1. Xの乱数を 0 〜 1 で生成（幅は 1）
x <- runif(N, min = 0, max = 1)

# 2. 正規分布の「高さ」を全員分計算する
# ※ 1 / sqrt(2 * pi) * exp(-x^2 / 2) をそのまま書きます
heights <- 1 / sqrt(2 * pi) * exp(-x^2 / 2)

# 3. 幅(1) × 高さの平均値
estimated_area <- 1 * mean(heights)

# 結果の表示
cat("モンテカルロ積分:", estimated_area, "\n")