﻿rm(list = ls())  # ワークスペースをクリア
z <- qnorm(0.95, 0, 1)
z  # z_0.05
mu <- 1.0
pnorm(z, sqrt(33)*mu/2, 1, lower.tail=FALSE)  # 検出力(n=33)
pnorm(z, sqrt(34)*mu/2, 1, lower.tail=FALSE)  # 検出力(n=34)
pnorm(z, sqrt(35)*mu/2, 1, lower.tail=FALSE)  # 検出力(n=35)
pnorm(z, sqrt(36)*mu/2, 1, lower.tail=FALSE)  # 検出力(n=36)